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ABSTRACT 


This  thesis  examines  the  pool  boiling  heat  transfer 
characteristics  of  HCFC-124  and  HCFC- 124/oil  mixtures  with  up 
to  10%  (by  weight)  miscible  alkylbenzene  oil.  One  smooth  and 
4  enhanced  tubes  were  tested:  a  19  and  26  low  integral -fin 
tube  (GEWA-K) ;  a  modified  finned  tube  (TURBO-B) ;  and  a  porous 
coated  tube  (HIGH  FLUX) .  The  tests  were  carried  out  using  the 
procedure  used  for  CFC-114  at  the  same  saturation  temperature 
of  2.2  ’C.  This  allowed  for  direct  comparison  of  the  pool 
boiling  heat  transfer  characteristics  between  the  two 
refrigerants . 

The  smooth  and  GEWA-K  19  fin  per  inch  tube  performance  in 
pure  HCFC-124  and  HCFC- 124/oil  mixtures  ranged  between  10  to 
50%  better  than  in  pure  CFC-114  and  CFC- 114/oil  mixtures  for 
all  heat  fluxes.  The  HIGH  FLUX  and  TURBO-B  tubes  were  similar 
in  performance.  With  pure  HCFC-124,  the  finned  tubes  typically 
provided  enhancements  in  the  heat  transfer  coefficient  between 
2  and  3  times  that  of  a  smooth  tube.  The  HIGH  FLUX  and  TURBO-B 
surfaces  typically  provided  additional  enhancements  2  times 
that  of  the  finned  tubes.  With  the  addition  of  oil,  the  heat 
transfer  increased  from  the  smooth  and  finned  tubes  but 
decreased  from  the  HIGH  FLUX  and  TURBO-B  tubes.  The  HIGH  FLUX 
and  TURBO-B  tubes  therefore  exhibited  enhancements  less  than 
the  finned  tubes  at  high  oil  concentrations  and  high  heat 
fluxes . 
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I.  INTRODUCTION 


In  an  ongoing  effort  to  eliminate  the  use  of 
chlorof luorocarbons  (CFCs) ,  which  have  been  linked  to  the 
depletion  of  the  earth's  protective  ozone  layer,  the  United 
States  Navy  is  looking  into  using  recently  developed 
alternative  refrigerants.  There  are  2  principal  alternative 
refrigerants:  hydrofluorocarbons  (HFCs),  and 
hydrochlorof luorocarbons  (HCFCs) .  HFCs  contain  no  chlorine  and 
have  zero  ozone  depletion  potential  (ODP) .  HCFCs  do  contain 
chlorine,  but  with  the  attached  hydrogen  atom  to  the  CFC 
structure,  the  molecule  is  less  stable  and  the  ozone  damaging 
chlorine  is  released  in  the  lower  atmosphere  prior  to  reaching 
the  ozone  layer.  The  ODP  of  HCFC-124  is  about  2%, 
significantly  less  than  the  70%  ODP  for  CFC- 114.  Figure  1.1, 
provided  by  Baehr  [Ref.  1],  displays  the  advantages  of  HCFC- 
124  over  CFC- 114  with  respect  to  both  ozone  depletion 
potential  (ODP)  and  global  warming  potential  (GWP) . 

The  US  Navy  is  currently  using  CFC-114  (R-114)  in  all  of 
its  larger  capacity  centrifugal  air  conditioning  units  (>100 
tons).  The  search  for  a  'drop-in'  replacement  has  identified 
HCFC-124  (R-124)  as  a  possible  alternative  to  R-114.  'Drop-in' 
here  means  operating  at  the  same  evaporating  temperature  as 
existing  systems.  Therefore  the  operating  parameters  of  R-124 
under  similar  conditions  as  present  R-114  chiller  systems 


alocarbon  Ozone  Depletion  Potentials  and  Global  Warming  Potentials 


must  be  taken  into  account.  For  example,  the  shell  must  be 
capable  of  withstanding  higher  pressures  (11  psig  for  R-124  at 
2.2  °C  compared  to  0  psig  for  R-114  at  2.2  °C) .  The  shutdown 
system  pressure  is  also  a  significant  factor  when  the  system 
is  subjected  to  higher  temperatures  from  a  surrounding  engine 
room  or  perhaps  a  casualty  condition.  A  temperature  of  100  °F 
would  result  in  a  system  vapor  pressure  of  31  psig  for  R-114 
compared  to  68  psig  for  R-124.  If  the  heat  exchanger  shell 
can  handle  the  increased  system  pressure,  hardware 
modifications  should  be  limited  to  alterations  such  as 
compressor  impeller  size  and  changes  to  sealing  materials  to 
maintain  leak- tight  integrity. 

If  hardware  modifications  can  be  made,  then  the  other 
major  consideration  is  the  heat  transfer  characteristics  of  R- 
124  compared  to  those  of  R-114.  This  is  not  such  a  simple 
comparison  as  the  two  fluids  may  exhibit  varying  boiling 
characteristics  for  each  type  of  enhanced  surface.  In 
addition,  the  oil  used  for  each  fluid  is  different  since  each 
must  be  miscible  in  its  particular  refrigerant  and  this  may 
significantly  affect  heat  transfer  performance. 

To  evaluate  and  compare  the  heat  transfer  performance  of 
R-124  with  that  of  R-114,  the  following  thesis  objectives  were 
established: 

l.  Install  a  new  fiberglass -strengthened  evaporator  and 
condenser  section  into  existing  apparatus  to  withstand 
the  higher  saturation  pressure  associated  with  R-124. 
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2.  Collect  pool  boiling  heat  transfer  data  of  pure  R-114  and 

R- 114/oil  mixtures  over  a  range  of  increasing  and 

decreasing  heat  fluxes  for  the  newly  manufactured  tubes 
and  check  their  repeatability  with  existing  R-114  data 
using  the  newly  modified  apparatus. 

3.  Modify  existing  data  reduction  program  to  account  for  the 
thermophysical  properties  of  R-124. 

4.  Collect  pool  boiling  heat  transfer  data  of  pure  R-124  and 

R- 124/oil  mixtures  over  a  range  of  increasing  and 

decreasing  heat  fluxes  using  a  smooth  tube  and  four 
enhanced  surface  tubes. 

5.  Study  effect  of  oil  on  boiling  hysteresis  for  a  smooth 
and  porous  coated  tube. 

6.  Compare  single  tube  heat  transfer  characteristics  of  R- 
124  with  previously  obtained  R-114  data  taken  by  Sugiyama 
[Ref.  2] . 
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II.  POOL  BOILING  OF  REFRIGERANTS 

A.  POOL  BOILING  CURVE 

The  pool  boiling  curve,  as  provided  by  Bar-Cohen  [Ref. 3], 
is  shown  in  Figure  2.1.  The  applicable  regions  within  the 
curve  where  single  tube  pool  boiling  heat  transfer  tests  of  R- 
114  and  R- 124  refrigerants  were  conducted  are  reviewed. 


Figure  2.1  Pool  Boiling  Curve 
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1.  Natural  Convection 


This  region  (a-b)  occurs  during  the  increasing  heat 
flux  phase  where  there  are  no  vapor  bubbles  being  produced  at 
the  tube  surface.  In  this  region,  heat  transfer  is  due  solely 
to  single-phase  free  convection  effects.  Several  correlations 
are  available  to  predict  heat  transfer  performance  in  this 
region,  including  the  Churchill  and  Chu  [Ref.  4]  correlation 
shown  below  in  equation  (1) . 

h=k/D0  [ .  6+  .  387  (g  •  j8  *©*D03/ v  •  ot) 1/6 /  [1+  ( 0 . 559/Pr)  9/16]  8/27]  2  (1) 


2.  Incipient  Boiling 

The  transition  from  the  natural  convection  region  to 
the  nucleate  boiling  region  is  the  subject  of  numerous 
studies.  The  superheat  (Twall-Tsat)  at  which  bubble  inception 
(point  b)  occurs  is  referred  to  as  the  onset  of  nucleate 
boiling  (ONB)  or  incipient  boiling.  Considerable  reduction  in 
tube  wall  temperature  occurs  at  this  point,  resulting  in  a 
considerable  reduction  in  superheat  noted  as  ATex  in  Figure 
2.1.  This  'jump'  in  performance  only  occurs  during 
increasing  heat  flux  and  is  referred  to  as  the  thermal 
overshoot.  A  hysteresis  effect  occurs  when  decreasing  heat 
flux  data  (from  point  c  to  point  a)  are  compared  to  the 
initial  ascending  path  (a-b-c) .  Several  factors  believed  to 
influence  the  thermal  overshoot  at  the  ONB  for  highly-wetting 
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liquids  were  studied  by  Bar-Cohen  [Ref.  3].  Probable 
influential  factors  included:  vibration,  surface  tension, 
enhanced  commercial  surfaces,  pressure,  time  between  test 
runs,  oil  mixture,  temperature  attained  at  end  of  decreasing 
heat  flux  phase  of  the  cycle,  and  variation  with  cool-down 
procedure.  Because  of  these  numerous  influences  on  the  boiling 
incipience  superheat,  past  studies,  including  Bar- Cohen  [Ref. 
3] ,  reveal  highly  variable  and  widely  distributed  superheats 
at  ONB.  Furthermore,  the  process  of  first  nucleation  is 
inherently  stochastic,  leading  to  variations  in  behavior  from 
one  day  to  the  next . 

3.  Nucleate  Bolling  Region 

Once  nucleate  boiling  commences,  during  increasing 
heat  flux  (point  c)  ,  sustained  improvement  in  heat  transfer  is 
realized  throughout  the  entire  fully- developed  nucleate 
boiling  region  (c-d) .  As  vapor  bubbles  form  at  nucleation 
sites  and  separate  from  the  surface,  dramatic  increases  in 
heat  transfer  are  attained.  There  are  several  factors 
affecting  where  this  region  is  located  on  the  pool  boiling 
curve,  the  most  important  of  which  is  the  number  of  active 
sites  on  the  surface,  as  may  occur  on  different  boiling 
surfaces.  Thome  [Ref.  5]  provides  detailed  insight  into 
enhanced  boiling  heat  transfer  from  various  commercially 
available  surfaces.  For  plain  tubes,  several  correlations  have 
been  proposed  to  predict  boiling  performance.  Two  such 
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correlations  were  provided  by  Stephan  and  Abdelsalam  [Ref.  6] 
and  Chongrungreong  and  Sauer  [Ref.  7]  . 

B.  SINGLE  TUBE  BOILING  EXPERIMENTS  AT 

NAVAL  POSTGRADUATE  SCHOOL 

Heat  transfer  measurements  were  collected  by  McManus 
[Ref. 8]  for  single  tube  pool  boiling  of  R-114  at  a  saturation 
temperature  of  13.8°C  with  up  to  6%  oil  concentration.  The 
boiling  tubes  in  these  tests  were  heated  internally  by  warm 
water.  HIGH  FLUX,  TURBO-B  and  smooth  tubes  were  tested  in  this 
manner.  McManus  [Ref.  8]  reported  enhancements  (compared  to 
the  smooth  tube)  in  pure  R-114  as  high  as  14.6  and  6.4  for  the 
HIGH  FLUX  and  TURBO-B  tubes  respectively,  while  performance 
decreased  to  enhancements  of  7.0  and  4.9  respectively  with  a 
6%  oil  concentration. 

Pool  boiling  heat  transfer  coefficients  in  pure  R-114  and 
R- 114/oil  mixtures  up  to  10%  (by  weight)  were  measured  by 
Wanniarachchi  et  al  [Ref.  9]  for  electrically-heated  smooth 
and  enhanced  tubes.  These  tests  were  run  at  a  saturation 
temperature  of  2.2°C.  At  3  0  kW/m2  Wanniarachchi  et  al  [Ref.  9] 
reported  enhancement  results  for  the  HIGH  FLUX,  Thermoexcel - E, 
Thermoexcel -HE,  and  GEWA-T  tubes  of  9.1,  8.2,  6.8,  and  4.4 
respectively.  Also  reported  were  the  enhancements  for  these 
tubes  with  3%,  and  10%  oil  concentration. 

Sugiyama  [Ref.  2]  conducted  extensive  pool  boiling  tests 
with  R-114  and  R- 114/oil  mixtures  using  a  large  variety  of 
enhanced  surface  tubes.  For  the  re-entrant  cavity  tubes  with 
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pure  R-114,  Sugiyama  [Ref.  2]  reported  enhancements  of  6.4  and 
6.2  for  the  HIGH  FLUX  and  TURBO-B  tubes  respectively  at  35 
kW/m:.  With  10%  oil  these  enhancements  were  reduced  to  3.8  and 
5.6  respectively.  This  degradation  in  performance  was 
partially  attributed  to  the  clogging  of  surface  pores  with 
oil.  For  GEWA-K  26  and  40  fpi  tubes,  enhancements  of  2.4  and 
3.1  were  reported  with  pure  R-114  at  35  kW/m;.  With  10%  oil 
these  same  tubes  were  surprisingly  reported  to  have 
enhancements  of  3.2  and  4.6  respectively.  This  indicated  an 
improved  heat  transfer  performance  with  addition  of  oil  in  the 
case  of  finned  tubes,  believed  to  be  caused  by  foaming  action 
created  within  the  pool  that  was  most  prevalent  at  higher  oil 
concentrations  (3-10%)  and  heat  fluxes  <>30  kW/m2)  . 

The  effect  of  oil  on  the  incipience  of  nucleate  boiling 
for  pure  R-114  and  R- 114/oil  mixtures  was  studied  by  Memory 
and  Marto  [Ref.  10]  with  smooth,  finned  and  re-entrant  cavity 
tubes.  They  observed  a  significantly  lower  heat  flux  at 
incipience  for  the  HIGH  FLUX  and  TURBO-B  tubes  compared  with 
the  smooth  and  finned  tubes  in  pure  R-114.  However,  consistent 
delay  in  the  boiling  incipience  of  the  HIGH  FLUX  and  TURBO-B 
tubes  was  exhibited  in  the  presence  of  3  and  10%  oil 
concentrations.  At  these  oil  concentrations,  boiling 
incipience  for  the  smooth  and  finned  tubes  showed  no 
dependence  with  oil.  Also  emphasized  was  the  inherent 
randomness  and  variation  connected  with  temperature  overshoot 
associated  with  boiling  incipience. 
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C.  COMPARISON  OF  R-114  TO  R-124 

1.  Property  Comparison  of  R-114/R-124 

Table  1  is  provided  to  show  direct  comparison  of  the 
properties  of  R-114  and  R-124  at  Tsat=2.2°C 


Table  1.  R-114/R-124  PROPERTY  COMPARISON 


Property 

R-114 

R-124 

Ps*/  (kPa) 

96 

177 

P  (liq) / (kg/m2) 

1526 

1429 

p  (vapor)  /  (kg/m3) 

am 

11.3 

cp (liq) / ( J/kg- K) 

932 

1060 

hfg/ (kJ/kg) 

134.7 

158.9 

k ( liq) / (W/m- K) 

0.0701 

0.0744 

M  (liq)  /  (g/m-  s) 

0.4449 

0.3375 

a/ (N/m) 

0.0136 

0.0129 

2.  Need  for  Oil  Mixtures 

The  oil  mixtures  used  are  to  simulate  the  actual 
environment  in  refrigeration  systems  which  can  collect  lube 
oil  within  a  unit's  flooded  evaporator  to  concentrations  as 
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high  as  10%.  The  oils  used  for  R-124  and  R-114  are  both 
miscible  within  their  respective  refrigerants.  However,  with 
R-114,  a  York-C  mineral  oil  was  used,  while  with  R-124,  an 
alkylbenzene  oil  was  used. 

3.  Need  for  R-124  Data 

There  is  presently  no  data  on  single  tube  pool  boiling 
heat  transfer  with  R-124  or  R- 124/oil  mixtures.  Tests  must 
therefore  be  conducted  with  R-124  to  parallel  R-114  results 
obtained  to  date. 
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III.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 


A.  SYSTEM  DESCRIPTION 

Details  of  the  relocation  and  early  modifications  to  the 
original  apparatus  are  given  by  Sugiyama  [Ref.  2]  and  Lake 
[Ref.  11].  To  accommodate  the  increased  operating  pressures 
associated  with  R-124,  stronger  evaporator  and  condenser  units 
were  installed.  The  data  reduction  program  DRP71  was  modified 
(DRPGB)  to  include  the  thermophysical  properties  of  R-124. 

A  schematic  diagram  representing  the  experimental 
apparatus  is  shown  in  Figure  3.1  and  comprises  the  following 
components:  1)  an  evaporator  tee  used  to  boil  the  refrigerant; 
(2)  a  condenser  tee  used  to  condense  the  refrigerant;  (3)  a 
reservoir  for  liquid  refrigerant;  (4)  a  ref rigerant/oil 
subsystem; (5)  a  refrigeration  cooling  support  system;  (6)  a 
vacuum  pump;  (7)  a  data  acquisition  and  instrumentation 
system;  (8)  an  aluminum/plexiglas  framework  housing  items  (1) 
to  (4)  . 

The  general  operation  of  the  system  was  closed  loop  with 
vapor  generated  in  the  evaporator  by  the  heated  tube  located 
within  the  refrigerant  pool  of  the  evaporator.  The  vapor 
passed  through  an  aluminum  ' L- shaped'  tube  to  the  condenser 
where  the  vapor  was  condensed  by  a  chilled  water/ethylene 


12 


Figure  3.1  Schematic  of  Apparatus 
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glycol  mixture  circulated  through  the  condenser  coil  by  an  8 
gpm  turbine  type  positive  displacement  pump.  The  R-114  or  R- 
124  condensate  was  then  gravity  fed  back  to  the  evaporator 
through  a  return  line.  R-502  and  R-12  refrigeration  systems 
maintained  the  water/ethylene  glycol  mixture  in  the  sump 
between  -12  and  -18  degrees  Celsius.  Sump  cooling  was 
accomplished  using  a  counter  flow  heat  exchanger  '  / 
recirculating  the  sump  coolant  (using  an  identical  positive 
displacement  pump)  over  the  R-502  side  of  the  heat  exchanger. 
To  further  maintain  this  coolant  temperature,  the  R-12  system 
was  located  directly  within  the  water/ethylene  glycol  sump. 

The  addition  of  oil  to  the  refrigerant  within  the 
evaporator  was  accomplished  by  using  a  graduated  cylinder  to 
measure  the  precise  quantity.  This  simply  drained  into  the 
evaporator  by  gravity. 

Accurate  horizontal  positioning  of  the  boiling  tube  was 
critical  and  accomplished  by  supporting  each  end  of  the  tube 
by  Teflon  inserts.  The  tube  was  sealed  using  0- rings.  For 
safety,  a  relief  valve  on  the  aluminum  pipe  between  the 
evaporator  and  condenser  was  set  at  a  gage  pressure  of  49  psi. 

B.  BOILING  TEST  SECTION 

1 .  Evaporator 

The  evaporator  was  a  Pyrex  glass  'T'  section  coated 
with  a  continuous  winding  of  fiberglass  filaments  impregnated 
with  a  modified  epoxy  resin.  This  upgraded  evaporator  shown  in 
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Figure  3.2  was  designed  for  a  working  pressure  of  50  psi  gage 
and  initially  hydrostatically  safety  tested  to  75  psi  gage. 
The  original  glass  was  designed  for  a  working  pressure  of  30 
psi  gage.  The  inside  surface  was  the  same  as  before  and 
maintained  the  advantages  of  minimizing  nucleate  boiling  at 
the  surface  as  well  as  being  corrosion  resistant.  The  two  ends 
of  the  evaporator  were  fitted  with  a  cast-iron  flange  and 
gasket  assembly  as  shown  in  Figure  3.2.  They  mated  with 
aluminum  flanges  at  each  end  of  the  evaporator,  which  in  turn 
contained  the  thermocouple  housings,  oil  entry  connection  and 
Teflon  inserts  containing  the  evaporator  tube. 

2 .  Evaporator  tubes 

Figure  3.3  shows  a  schematic  of  a  boiling  tube.  These  were 
accurately  held  in  place  using  Teflon  inserts  and  sealed  by  0- 
rings.  A  4 -stud  strong-back  secured  the  Teflon  inserts  inside 
the  aluminum  flanges. 

Five  different  tubes  were  tested  within  the  apparatus. 
Dimensions  of  each  tube  are  provided  in  Table  5  of  Appendix  C 
while  an  expanded  view  of  the  enhanced  surfaces  are  shown  in 
Figure  3.4.  The  tubes  used  were: 

1 .  Smooth  tube 

2.  GEWA-K  (19  fpi)  finned  tube 

3.  GEWA-K  (26  fpi)  finned  tube 

4.  HIGH  FLUX  (porous  coated)  enhanced  tube 

5.  TURBO-B  (deformed  surface)  enhanced  tube 


15 


a)  Corning  Pyrex  Glass  Evaporator  CD  x  Ii=  402x51  mm 
L  =  178  mm,  L*  *  127  mm) 


INSERT 

b)  Cast  Iron  Flange  and  Gasket  (d**  190mm,  do=  210mm 
1^  =  14mm,  A*  21*  ) 


Figure  3.2  Sketch  of  Pyrex  Glass  Vessels  and  Flanges 


Figure  3.3  Schematic  of  boiling  tube 
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HIGH  FLUX  SURFACE 


Figure  3.4  Surface  Profile#  of  Enhanced  Tube# 
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The  heater  used  within  each  boiling  tube  was  nominally  a 
1  kW,  240-Volt  stainless  -  steel  cartridge  heater,  measuring 
6.35  mm  in  outside  diameter  and  203.2  mm  in  length.  The  actual 
heated  length  of  the  heater  was  190  mm.  The  HIGH  FLUX  and 
GEWA-K  26  fpi  tubes  were  previously  manufactured  while  the 
TJRBO-B,  GEWA-K  19  fpi  and  smooth  tubes  required  assembly 
prior  to  testing  with  R-124  or  R-114.  The  manufacturing 
process  of  these  tubes  is  described  below. 

The  heater  was  placed  inside  a  copper  sleeve  (Figure  3.5 
provides  dimensions  of  smooth  tube  sleeve)  with  a  tight 
mechanical  fit  of  .001  inch  clearance.  The  copper  sleeve  was 
machined  with  6  grooves  (1.27mm  x  1.27mm)  on  its  outer 
diameter  of  various  lengths.  Thermocouples  were  laid  in  each 
of  these  grooves  and  the  edges  of  the  grooves  were  peened  over 
with  a  punch.  Any  resulting  imperfections  were  smoothed  out 
with  250  grit  emery  paper.  The  thermocouple  wires  and  heater 
leads  all  exited  the  copper  sleeve  at  the  same  end.  The  copper 
sleeve  containing  the  heater  and  thermocouples  was  then  placed 
on  2  aluminum  V-blocks  and  the  outer  surface  of  the  sleeve  was 
covered  with  a  flux  and  solder  paste  (50/50  tin/lead) ;  great 
care  was  taken  to  ensure  complete  coverage  of  the  thermocouple 
wires  lying  in  the  grooves.  The  thermocouple  and  heater  leads 
were  then  passed  through  the  enhanced  boiling  tube  (the  outer 
diameter  of  the  copper  sleeve  was  machined  to  fit  tightly  into 
each  respective  tube) .  To  prevent  shorting  between  the  heater 
and  outer  tube  shell,  a  high  temperature  RTV  silicone  rubber 
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(DIMENSIONS  IN  mm) 


Figure  3 . 5  Schematic  of  copper  sleeve 
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was  applied  to  the  heater  and  thermocouple  leads  at  the  point 
where  they  leave  the  sleeve.  The  copper  sleeve  was  then 
tapped  gently  into  the  tube  and  positioned  accurately  so  as  to 
line  up  precisely  with  the  8  inch  active  length  of  the 
enhanced  surface.  This  enhanced  surface  length  was  different 
from  the  actual  heated  length  which  was  only  7.5  inches. 

The  entire  assembly  was  then  placed  horizontally  in  an 
oven  and  heated  to  between  400°F  and  500°F.  This  was  the 
temperature  'window'  at  which  the  flux  melted,  but  the 
thermocouple  leads  did  not.  The  temperature  of  the  tube  was 
monitored  by  attaching  the  thermocouples  to  a  simple 
temperature  indicator.  After  stabilizing  the  oven  temperature 
between  these  temperatures  and  ensuring  that  the  solder  had 
melted,  the  tube  was  taken  out  and  quenched  with  a  damp  cloth 
at  one  end  (as  close  to  the  end  as  possible)  where  the 
thermocouple  and  heater  leads  protrude.  The  oven  was  then 
placed  vertically  and  the  tube  reinserted  from  below  keeping 
the  lower  end  quenched  while  holding  both  tube  and  cloth  with 
a  pair  of  channel  locks .  When  the  temperature  increased  to 
>400°F,  the  solder  again  began  to  melt.  Additional  40/60  tin- 
lead  solder  was  then  applied  by  hand  from  the  end  which 
protrudes  from  the  top  of  the  oven.  Remelting  and  addition  of 
extra  solder  in  this  way  ensured  a  uniform  circumferential 
layer  of  solder  and  helped  to  prevent  formation  of  air 
pockets.  The  tube  was  then  removed  from  the  oven  and  quenched 
over  its  whole  length  with  a  damp  cloth  while  in  a  vertical 
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position.  It  was  then  left  to  cool  in  air  for  24  hours.  This 
procedure  differed  from  earlier  manufactured  tubes  in  that  the 
previous  copper  sleeves  had  4  thermocouples  leaving  each  end 
of  the  sleeve.  Trying  to  repeat  this  previous  procedure 
accurately  led  to  large  discrepancies  in  the  thermocouple 
readings  (>20°C)  as  heater  power  was  increased.  It  was  felt 
this  was  due  to  air  pockets  created  in  the  fabrication 
procedure  (the  tube  could  never  be  placed  vertically  in  the 
oven  with  leads  protruding  from  both  ends)  .  The  above 
procedure  was  much  simpler  and  led  to  excellent  uniformity  in 
the  thermocouple  values  at  all  heater  powers. 

C .  CONDENSER  SECTION 

The  condenser  was  also  a  fiberglass  coated  Corning  Pyrex 
glass  Tee  identical  in  size  and  shape  to  the  evaporator  shown 
in  Figure  3.2.  The  refrigerant  was  condensed  on  a  3/8  inch 
outside  diameter  copper  tube  that  was  coiled  within  the 
condenser.  The  condenser  was  mounted  vertically.  An  O-ring 
sealed  aluminum  flange  mounted  at  the  top  of  the  condenser 
contained  not  only  the  water  ethylene/glycol  coolant 
inlet/outlet,  but  also  a  3/8  inch  outside  diameter  copper  tube 
vent,  which  was  connected  either  to  a  vacuum  pump  via  valve 
Vll  or  to  a  charging  cylinder  via  valve  V10.  The  bottom  of  the 
condenser  was  also  capped  with  an  O-ring  sealed  aluminum 
flange  and  contained  the  gravity  feed  return  line  to  the 
evaporator  via  valve  V5 .  The  refrigerant  vapor  from  the 
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evaporator  entered  the  condenser  through  an  aluminum  L- shaped 
tube,  connected  to  the  condenser  by  O-ring  and  gasketed 
flange.  A  pressure  gage  ranging  from  30  inches  mercury  vacuum 
to  150  psig  was  mounted  on  the  L- shaped  tube;  this  was  used 
only  as  a  rough  visual  check  of  pr'^sure  within  the  facility. 

D.  OIL  ADDITION  SECTION 

To  provide  oil  addition  into  the  refrigerant,  a 
cylindrical  aluminum  reservoir,  6  inches  in  diameter  and  6 
inches  high,  and  a  graduated  cylinder  355  mm  in  length  and 
25.4  mm  in  diameter  is  used  as  shown  in  Figure  3.1.  The  oil 
used  with  R-124  was  a  miscible  alkylbenzene  oil  while  the  oil 
used  with  R-114  was  a  York-C  mineral  oil,  also  miscible.  The 
graduated  cylinder  was  filled  from  the  reservoir  via  valve  V2 , 
while  the  reservoir  was  replenished  through  valve  V3  at  the 
top.  Oil  entered  the  evaporator  from  the  graduated  cylinder  by 
gravity  feed  via  valve  VI.  Oil  concentrations  (by  weight)  of 
0,1, 2, 3, 6,  and  10%  were  used  for  each  of  the  tubes. 

E.  COOLING  SYSTEM 

1.  Water /Ethylene- Glycol  Coolant  Tank 

At  the  base  of  the  apparatus,  a  5.35  cubic  foot  tank 
stored  30  gallons  of  ethylene-glycol/water  mixture. 
Recirculation  and  discharge  ports  were  connected  to  the  tank 
together  with  a  thermocouple  to  monitor  sump  temperature. 
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2 .  Refrigeration  Plants 

To  provide  the  cooling  of  the  water/ethylene  glycol 
mixture,  a  1/2  ton  R-502  and  1/4  ton  R-12  refrigeration  system 
were  installed  near  the  sump  as  seen  in  Figure  3.1.  Both 
plants  contained  a  hermetically  sealed  compressor  assembly, 
an-air  cooled  condenser,  receiver,  filter-dryer,  pressure 
regulator,  temperature  control  switch,  and  thermostatic 
expansion  valve.  The  R-502  plant  (Figure  3.6)  used  a  counter 
current  heat  exchanger  while  the  R-12  plant  used  a  coiled 
evaporator  placed  within  the  sump.  The  sump  temperature  was 
maintained  at  approximately  -17  °C  when  both  these  systems 
were  used.  Each  plant  was  controlled  by  a  temperature  control 
switch  and  thermostatic  expansion  valve. 

3 .  Pump  and  Control  Valve 

An  8  gpm,  positive  displacement  pump  was  installed 
with  a  1  inch  diameter  PVC  pipe  connected  to  the 
water/ethylene  glycol  sump.  The  discharge  of  the  pump  was 
piped  to  the  condenser  via  valve  VC  as  shown  in  Figure  3.1.  A 
by-pass  valve  V9  was  placed  upstream  of  the  control  valve  VC 
to  avoid  overloading  the  pump  when  control  valve  VC  was  closed 
and  provides  for  mixing  of  the  return  coolant  from  the 
condenser.  At  high  heat  fluxes,  the  by-pass  valve  needed  to  be 
closed  down  to  increase  flow  of  coolant  through  the  condenser. 
An  identical  pump  discharged  the  ethylene  glycol/water  coolant 
through  the  R-502  counter- current  heat  exchanger. 
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Figure  3.6  R-502  Refrigerant  plant 
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F.  REFRIGERANT  RESERVOIR 


The  reservoir  used  to  store  the  R-114  or  R-124  was  made 
from  an  aluminum  cylinder,  9  inches  in  diameter  and  10  inches 
in  height.  A  sight  glass  attached  on  the  reservoir  allowed  for 
the  liquid  level  of  refrigerant  to  be  monitored.  The  reservoir 
was  connected  to  the  'L'  shaped  tube  through  valve  V7  and  to 
the  liquid  line  into  the  evaporator  through  valve  V6  as  shown 
in  Figure  3.1. 

G.  FRAME 

To  house  the  entire  assembly,  an  aluminum  frame  was 
constructed  measuring  1070.1  mm  x  509.8  mm  x  610.1  mm.  With  a 
view  to  keeping  the  apparatus  in  an  isothermal  atmosphere,  as 
well  as  for  safety  considerations,  the  sides  were  covered  with 
12.7  mm  thick  plexiglas  with  hinges  mounted  on  the  2  sides  to 
provide  easy  access  to  the  assembly.  The  top  of  the  frame  was 
covered  by  plywood  while  the  bottom  was  covered  by  an  aluminum 
plate.  The  entire  frame  rested  on  the  coolant  sump  tank.  The 
front  side  of  plexiglas  had  small  holes  to  accommodate  the 
valve  stems  of  valves  VI  through  V8,  so  they  can  be  operated 
while  maintaining  the  integrity  of  the  enclosed  frame. 

H.  INSTRUMENTATION 

1.  Power  Measurement 

The  boiling  tube  heater  was  powered  by  a  240  volt  AC 
source,  adjusted  to  0-220  volts  and  0-5  amps  by  a  variac 
control.  Power  input  to  the  heater  was  measured  with  an  AC 


26 


Figure  3.7  Data  acquisition/control  unit 
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current  inductive  sensor  and  a  voltage  sensor.  The  output  was 
a  proportional  signal  in  volts  having  been  converted  through 
an  AC-DC  true  R.M.S.  converter.  Both  the  AC  current  sensor  and 
the  AC-DC  true  R.M.S.  converter  were  connected  to  the  data 
acquisition/control  unit  as  shown  in  Figure  3.7. 

2 .  Temperature  Measurement 

The  primary  means  of  obtaining  temperatures  was  the 
use  of  30  gage  copper- constantan  thermocouples.  Specific 
temperatures  recorded  were  from  6  thermocouples  in  the  boiling 
tube  (8  in  the  earlier  manufactured  tubes) ,  3  thermocouples  in 
the  evaporator  pool  (housed  within  separate  sealed 
thermocouple  wells  within  the  evaporator  pool) ,  and  one 
thermocouple  to  monitor  the  water/ethylene  glycol  sump  tank. 
A  Hewlett-Packard  3497A  data  acquisition  system  controlled  by 
a  Hewlett-Packard  9826  computer  read  these  thermocouple 
measurements . 

I.  DATA  ACQUISITION  AND  REDUCTION  PROGRAM 

Data  collection  and  reduction  program  DRP71,  used  by 
Sugiyama  [Ref. 2]  for  use  with  R-114  was  modified  for  use  with 
R-124.  The  thermophysical  properties  of  saturated  R-124  were 
obtained  from  Ref  Prop  software  [Ref.  12]  and  were  used  to 
generate  equations  for  use  in  the  calculations  of  the 
reduction  program.  Appendix  A  contains  plots  of  each  property 
versus  temperature  and  the  respective  equation.  The  modified 
program  was  renamed  DRPGB . 
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The  Hewlett-Packard  9826  computer  reduced  all  the 
thermocouple,  current,  and  R.M.S.  voltage  readings  using  the 
DRPGB  program.  The  program  was  run  through  the  keyboard 
allowing  regulation  of  the  data  acquisition  unit.  A  printout 
of  the  reduced  data  was  provided  by  a  Hewlett-Packard  Inkjet 
printer.  Table  2  lists  the  channel  assignments  for  the  sensor 
inputs  to  the  data  acquisition  system. 


J.  PROCEDURE  PROM  KEYBOARD 

The  following  procedure  outlines  the  procedure  used  to 
collect  and  process  data  from  the  Hewlett  Packard  data 
acquisition  system: 


1.  Selected  the  'taking  data'  option  after  initial 
loading  of  program  DRPGB. 

2.  Selected  the  refrigerant  being  tested,  tube  heating 

mode,  (ie  electrically  or  water  heated)  and 

thermocouple  type  from  DRPGB. 

3.  Assigned  name  to  file  of  data  to  be  processed  and 
stored. 

4.  Identified  and  entered  any  defective  thermocouples  by 

channel  #  1-8,  noting  that  the  new  tubes  have  only  6 
thermocouples  and  channels  7  &  8  were  always 

' defective' . 

5.  Selected  test  tube  type. 

6.  Set  desired  saturation  temperature  (2.22°C)  for 

present  tests. 

7.  Set  desired  heat  flux  desired  using  variac  adjustment. 

8.  Attained  desired  saturation  temperature  by  adjusting 
flow  of  coolant  through  condenser  coils  with  control 
valve  VC. 
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9.  Once  saturation  temperature  was  achieved,  waited  5 
minutes  for  steady  state  conditions  prior  to  taking 
data . 

10.  Prompted  data  acquisition  unit  to  scan  all  the 
channels  listed  in  Table  2.  All  channel  readings  were 
made  in  volts  and  stored  in  user  specified  fields. 

11.  Calculated  temperature  and  power  from  these  voltages. 

12.  Calculated  heat  transfer  rate  from  cartridge  heater. 

13.  Calculated  the  average  wall  temperature  of  the  test 
tube  and  calculated  wall  superheat  (Twall-  Tsat) . 

14.  Calculated  physical  properties  of  R-124  using  property 
correlations  from  Ref Prop  [Ref .12]  at  film  temperature 
(Twall+Tsat) / 2 . 

15.  Calculated  the  natural  convection  heat  -  transfer 
coefficient  of  R-124  from  the  unheated  ends  of  the 
test  tube. 

16 .  Computed  heat  losses  from  unheated  ends  of  the  test 
tube . 

17.  Calculated  heat  flux  from  heated  length  of  the  test 
tube . 

18.  Calculated  the  heat  -  transfer  coefficient  from  the 
heated  length  of  the  test  tube. 

19.  Stored  heat  flux  and  wall  superheat  for  each  data  set 
in  user  specified  fields. 

20.  Plotted  data  on  available  software. 


Sample  calculations  of  the  above  procedure  are  given  in 
Appendix  C. 
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TABLE  2.  HP  3497A  CHANNEL  ASSIGNMENTS 


Channel 

Channel  Assignments 

00-07 

tube  wall  temperatures 

08-10 

pool  liquid  temperature 

11 

sump  temperature 

20 

RMS  voltage 

21 

current  sensor 
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IV.  EXPERIMENTAL  PROCEDURE 


A.  ASSEMBLY  PREPARATION 

1.  Vacuum  test  of  the  apparatus 

After  placing  the  desired  boiler  tube  in  position 
within  the  evaporator,  the  system  was  evacuated  down  to  a 
pressure  of  29  inch  Hg.  vacuum  using  a  portable  vacuum  pump 
via  valves  Vll  and  V8  in  Figure  3.1.  The  system  was  then 
isolated  for  approximately  30  minutes  to  check  for  leak- tight 
integrity  (indicated  by  any  significant  vacuum  loss) .  If  a 
leak  was  detected,  a  positive  air  pressure  (20  psig)  was 
placed  on  the  system  via  valves  V10  and  V8.  A  soap/water 
solution  was  applied  on  the  sealing  areas  to  identify  any 
leaks  (if  there  is  a  leak,  bubbles  are  seen) .  If  identified, 
the  leak  source,  was  corrected  and  the  system  integrity 
retested. 

2.  Charging  apparatus  with  refrigerant 

After  system  in  egrity  was  ensured,  the  evaporator  was 
filled  with  refrigerant  to  a  level  of  20  mm  above  the  test 
tube.  The  refrigerant  was  drained  by  gravity  from  the 
reservoir  into  the  evaporator  through  valves  V6  and  V5 .  Valve 
V7  was  also  open  to  equalize  pressure  in  the  system.  The 
reservoir  was  then  isolated  from  the  system  with  the  closure 
of  valves  V5  through  V7. 
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3 .  Degassing  and  Data  Acquisition  Channel  Check 

Prior  to  connecting  all  thermocouple  leads,  the 
resistance  between  the  heater  leads  and  outer  tube  wall  was 
checked  with  a  simple  multimeter.  Readings  were  usually 
greater  than  1  .  Each  channel's  output  was  checked  using 

program  SETUP71.  Any  problem  thermocouples  were  either  fixed 
or  excluded  from  the  calculations  in  program  DRPGB.  Before 
energizing  the  heater  and  starting  the  tests,  the  evaporator 
pool  was  slowly  cooled  to  a  saturation  pressure  of  2.22  °C  by 
allowing  a  low  coolant  flow  through  the  condenser.  Once 
completed,  the  heating  surface  flux  was  set  to  approximately 
90  kW/m2  to  degas  the  refrigerant;  any  non- condensible  gases 
were  collected  in  the  condenser  and  vented  off  using  the 
vacuum  pump  for  approximately  one  minute.  The  apparatus  was 
then  secured  for  a  period  of  approximately  8  hours  prior  to 
operation  to  allow  the  surface  to  cool  and  become  fully  wetted 
in  the  refrigerant. 

B .  OPERATIONAL  PROCEDURE 

The  following  procedure  was  used  to  operate  the  facility 
(this  follows  closely  the  procedure  used  by  Sugiyama  [Ref .2]  !  : 

1.  The  R-502  and  R-12  refrigeration  units  were  utilized 
until  the  temperature  of  the  water- ethylene  glycol  sump 
reached  a  minimum  of  -8  °C  before  starting  tests  (with 
continued  cooling  to  -17  °C) . 

2.  The  data  acquisition  unit,  computer  and  variac  control 
panel  were  switched  on. 
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3.  The  computer  program  SETUP  was  loaded  and  run.  All  data 
acquisition  channels  were  rechecked.  The  refrigerant 
pool  temperature  was  then  cooled  down  slowly  to  2.2  °C 
by  circulating  coolant  through  the  condenser. 

4.  The  data  acquisition/reduction  program  DRPGB  was  loaded 
and  run. 

5.  The  desired  heat  flux  setting  for  the  test  tube  was 
input  to  the  program  DRPGB  and  the  variac  was  adjusted 
to  attain  that  heat  flux.  Starting  heat  flux  was 
between  500  and  600  W/m2 . 

6.  Control  valve  VC  was  adjusted  to  control  coolant  flow 
through  the  condenser  to  maintain  a  constant  saturation 
temperature  of  2.22  °C  +  0.1  °C. 

7.  At  each  heat  flux  setting,  the  thermocouple  values  were 
allowed  to  stabilize  for  about  5  minutes  prior  to 
taking  data. 

8.  The  following  data  were  then  taken  and  used  in  the 
acquisition/reduction  program  DRPGB:  pool  temperature ; 
test  tube  wall  temperatures;  sump  temperature;  current 
and  voltage  sensor  readings. 

9.  Data  were  taken  at  each  desired  heat  flux  setting  and 
always  at  the  same  saturation  temperature  (the  coolant 
flow  in  the  condenser  was  adjusted  to  ensure  this) .  At 
each  heat  flux  setting,  program  DRPGB  generated  a 
permanent  printout  listing:  wall  temperatures  of  the 
test  tubes;  refrigerant  pool  temperatures;  sump 
temperature;  wall  superheat;  heat  -  transfer  coefficient 
and  the  corrected  heat  flux  (ie.  to  account  for  heat- 
transfer  lost  from  each  end) . 

10.  For  each  heat  flux  setting,  steps  6  through  9  were 
repeated.  The  superheat  (Twall-Tsat)  and  the  respective 
heat  flux  were  plotted  on  a  log -log  scale. 

11.  Table  3  provides  a  summarized  listing  of  all  data  runs. 
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TABLE  3.  LISTING  OF  DATA  RUNS 


Data  File 

Refrigerant 

Tube  Type 

Purpose 

D0821SM0 

R-  114 

Smooth 

repeatability 

D0824SM0 

R-  114 

Smooth 

repeatability 

D1015SM0 

R-  124 

Smooth 

Data 

D1017SM0 

R- 124 

Smooth 

Data 

D1223SM0 

R-  124 

Smooth 

Data 

D1222SM0 

R-  124 

Smooth 

Data 

D1221SM0 

R- 124 

Smooth 

Data/samp  calc 

D1220SM0 

R- 124 

Smooth 

uncert  analys 

D1219SM0 

R-  124 

Smooth 

Data 

D1217SM0 

R-  124 

Smooth 

Data 

D1016SM4 .44* 

R-  124 

Smooth 

Data  Tsat4.4°C 

D1114SM0 

R- 124 

Smooth 

Data 

D1114SM1 

R- 124 

Smooth 

Data 

D1114SM2 

R  '124 

Smooth 

Data 

D0108SM3 

R- 124 

Smooth 

Data 

D0107SM3 

R- 124 

Smooth 

Data 

D0106ASM3 

R- 124 

Smooth 

Data 

D0106SM3 

R- 124 

Smooth 

Data 

D0105SM3 

R- 124 

Smooth 

Data 

D0104SM3 

R- 124 

Smooth 

Data 

D1227SM3 

R- 124 

Smooth 

Data 

D1116SM3 

R- 124 

Smooth 

Data 

D1116SM6 

R-  124 

Smooth 

Data 

D0112ASM10 

R-  124 

Smooth 

Data 

DO 112SM10 

R-  124 

Smooth 

Data 

D0111ASM10 

R- 124 

Smooth 

Data 

D0111SM10 

R-  124 

Smooth 

Data 

Continuation  of  Table  3 


Data  File 

Refrigerant 

tube  type 

Purpose 

D0110SM10 

R- 124 

Smooth 

Data 

D0109ASM10 

R-  124 

Smooth 

Data 

D0109SM10 

R- 124 

Smooth 

Data 

D1117SM10 

R- 124 

Smooth 

Data 

D1011190 

R- 114 

GK-19 

Data 

D1010190 

R-  114 

GK-19 

Data 

D0901190 

R-  114 

GK-19 

Data 

D0907193 

R-  114 

GK-19 

Data 

D09091910 

R-  114 

GK-19 

Data 

D1105190 

R-  124 

GK-19 

Data 

D1118190 

R-  124 

GK-19 

Data 

D1118191 

R-  124 

GK-19 

Data 

D1118192 

R-  124 

GK-19 

Data 

D1119193 

R- 124 

GK- 1*3 

Data 

D1120196 

R- 124 

GK-19 

Data 

D11211910 

R- 124 

GK-19 

Data 

D0803260 

R- 114 

GK-26 

Repeatability 

D1021260 

R-  124 

GK-26 

Data 

D1122260 

R-  124 

GK-26 

Data 

D1122261 

R- 124 

GK-26 

Data 

D1122262 

R-  124 

GK-26 

Data 

D1123263 

R-  124 

GK-26 

Data 

D1123266 

R- 124 

GK-26 

Data 

D11252610 

R- 124 

GK-26 

Data 

D11242610 

R-  124 

GK-26 

Data 

D1004TB0 

R- 114 

TURBO -B 

Data 

D1006TB0 

R- 114 

TURBO -B 

Repeatability 

Continuation  of  Table  3 . 


Data  File 

Refrigerant 

Tube  Type 

Purpose 

D1002GYX0 

R- 114 

GEWA- YX 

Data 

D1110TB0 

R- 124 

TURBO -B 

Data 

D1110TB0 

R-  124 

TURBO  -  B 

Data 

D1111TB1 

R-  124 

TURBO -B 

Data 

D1111TB2 

R-  124 

TURBO -B 

Data 

D1112TB3 

R- 124 

TURBO -B 

Data 

D1112TB6 

R- 124 

TURBO -B 

Data 

D1113TB10 

R- 124 

TURBO -B 

Data 

DO 119AHF0 

R-  124 

HIGH  FLUX 

Data 

D0118BHF0 

R-  124 

HIGH  FLUX 

Uncert  Analys 

D0118HF0 

R-  124 

HIGH  FLUX 

Data 

D0114HF0 

R-  124 

HIGH  FLUX 

Data 

D0113AHF0 

R-  124 

HIGH  FLUX 

Data 

D0113HF0 

R-  124 

HIGH  FLUX 

Data 

D1108HF0 

R-  124 

HIGH  FLUX 

Data 

D1126HF0 

R- 124 

HIGH  FLUX 

Data 

D1126HF1 

R- 124 

HIGH  FLUX 

Data 

D1126HF2 

R-  124 

HIGH  FLUX 

Data 

D0125HF3 

R- 124 

HIGH  FLUX 

Data 

DO 122HF3 

R-  124 

HIGH  FLUX 

Data 

D0121CHF3 

R- 124 

HIGH  FLUX 

Data 

D0121HF3 

R- 124 

HIGH  FLUX 

Data 

D0120AHF3 

R- 124 

HIGH  FLUX 

Data 

D0120HF3 

R-  124 

HIGH  FLUX 

Data 

D0119AHF3 

R- 124 

HIGH  FLUX 

Data 

D1128HF3 

R- 124 

HIGH  FLUX 

Data 

D1128HF6 

R- 124 

HIGH  FLUX 

Data 
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Continuation  o£  Table  3 . 


Data  File 

Refrigerant 

Tube  Type 

Purpose 

D1130HF10 

R-124 

HIGH  FLUX 

Data 

D0128HF10 

R-124 

HIGH  FLUX 

Data 

DO 127CHF10 

R-124 

HIGH  FLUX 

Data 

D0127BHF10 

R-124 

HIGH  FLUX 

Data 

D0127HF10 

R-124 

HIGH  FLUX 

Data 

D0126BHF10 

R-124 

HIGH  FLUX 

Data 

D0126AHF10 

R-124 

HIGH  FLUX 

Data 

D0126HF10 

R-124 

HIGH  FLUX 

Data 

The  tabulated  data  was  filed  using  the  following  file  name 
system: 


Example  (D0126AHF10) : 

1.  First  letter  D  simply  refers  to  data 

2.  The  following  4  characters  are  always  the  date;  in  this 
case  Jan  26. 

3 .  The  next  letter  represents  order  of  the  day  test  was 
completed  on  the  respective  date; 

no  letter-  1st  run 
A-  2nd  run 
B-  3rd  run 
etc . 

4.  The  following  2  letters  represented  the  tube  used; 

SM=  smooth 
19=  19  fpi  GEWA-K 
26=  26  fpi  GEWA-K 
HF=  HIGH  FLUX 
TB=  TURBO -B 

5.  The  last  1  or  2  numbers  represent  the  percentage  of 
oil;  in  this  case  10%. 

*  this  one  run  was  pure  R-124  tested  at  4.4  °C  vice  2.2  °C 
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V.  RESULTS  AND  DISCUSSION 


A.  REPRODUCTION  OP  EXPERIMENTAL  DATA 

The  first  tests  conducted  were  to  verify  the 
reproducibility  of  the  single  tube  data  taken  by  Sugiyama 
[Ref.  2]  with  R-114.  A  GEWA-K  26  fin  per  inch  (fpi)  tube 
comparison  with  R-114  is  given  in  Figure  5.1  and  shows  very 
similar  results  with  previous  data.  Uncertainty  bands  for  both 
wall  superheat  and  heat  flux  are  given  on  the  figure. 
Consequently,  at  low  heat  flux  (<  1000  W/m2)  ,  there  is  more 
scatter  due  to  the  lower  wall  superheats  and  hence  higher 
uncertainty.  Overall,  the  agreement  between  the  two  data  sets 
is  satisfactory. 

In  Figure  5.2,  an  increasing  heat  flux  data  run  using 
a  smooth  tube  is  compared  with  Sugiyama' s  [Ref.  2]  data.  It 
can  be  seen  that  in  the  low  heat  flux,  nature  1  convection 
region  and  the  high  heat  flux,  nucleate  boiling  region  the 
data  agreement  is  good.  In  the  mid-range  of  heat  flux,  (3 
kW/m2<  q  <  15  kW/m2)  ,  the  data  varies  somewhat,  depending  upon 
how  the  tube  nucleates.  Curve  #1  nucleates  at  a  relatively  low 
wall  superheat  (12.5  °C)  while  curve  #2  (both  taken  by 
Sugiyama  [Ref.l])  nucleates  at  a  high  wall  superheat  of  about 
35  °C  (point  C) .  The  present  data  (curve  #3)  exhibits  partial 
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Figure  5.1  Comparison  of  R-114  Data  with  0%  Oil  for  GEWA-K 
26  fpi  tube  (decreasing  flux) 


40 


iiiim 


CM 


Figure  5.2  Comparison  of 
Tube  (increasing 


% 

♦ 

5 

x 

♦ 

□ 

♦ 


□ 


♦  □ 

X  □ 


* 


<3 


♦ 


fl.R-IU.Smooth.SufllyamaJHaf  2] 
#2, R-114, Smooth, Sufllyama{R#f  2] 
#3,R-114,  Smooth 


T  I  Mil! 


1  10  100 

Wall  Suparhaat  (K) 


R-114  Data  with  0%  Oil  for  Smooth 
flux) 


41 


nucleation  (see  point  A)  where  3  of  the  six  thermocouples  had 
a  low  wall  superheat  (<  14  °C)  and  the  remaining  3  had  a  wall 
superheat  (>  18  °C) .  Observing  the  tube  surface,  it  could  be 
clearly  seen  that  half  the  tube  was  boiling  while  the  other 
half  was  in  the  natural  convection  mode.  At  point  B,  all  six 
thermocouples  were  within  0.9  °C  of  each  other,  indicating 
complete  nucleation  over  the  entire  surface.  At  this  point 
(and  above),  agreement  with  the  previous  data  is  excellent. 
The  discrepancies  in  this  mid  heat  flux  region  demonstrate  the 
random  behavior  of  nucleate  boiling  incipience;  this  is 
further  investigated  later  in  this  chapter.  After  all  sites 
nucleated  (Point  B)  ,  the  standard  deviation  was  0.9  °C  with  an 
average  thermocouple  temperature  of  12.9°C.  The  point  of 
incipience  and  its  random  occurrence  at  different  levels  of 
superheat  will  be  examined  later  in  this  chapter. 

Sugiyama  [Ref  .2]  did  not  conduct  tests  on  a  19  fpi  GEWA- 
K  tube.  This  is  one  of  the  most  common  tubes  used  in 
commercial  evaporators.  Therefore,  base-line  data  for  the  19 
fpi  GEWA-K  tube  with  R-114  and  R-114/oil  mixtures  were  taken 
during  this  thesis  before  evacuating  and  cleaning  the 
apparatus,  and  filling  with  R-124.  The  R-114  data  obtained  was 
incorporated  in  the  R-124/R-114  comparison  of  the  19  fpi  tube 
( Figure  5.12)  . 

Figure  5.3  shows  the  repeatability  of  decreasing  heat 
flux  data  for  the  newly  constructed  19  fpi  GEWA-K  tube  using 
R-124.  The  good  agreement  is  an  indication  of  the  consistency 
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Figure  5.3  Repeatability  Comparison  of  R*124  Data  with  0%  Oil 
for  19  fpi  GEWA-K  Tube  (decreasing  flux) 
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of  the  data  reduction  program  DRPGB,  the  apparatus  and  the 
operator.  Decreasing  heat  flux  runs  provide  better  comparison 
as  temperature  overshoot  and  hysteresis  effects  are  avoided. 

B.  POOL  BOILING  HEAT -TRANSFER  COEFFICIENT  CORRELATIONS 

Among  the  correlations  available  to  predict  pool  boiling 
heat  transfer  coefficients  for  a  smooth  tube,  most  can  be 
expressed  as  heat  flux  as  a  function  of  wall  superheat.  The 
pure  R-124  smooth  tube  data  is  fitted  very  closely  by  the 
equation: 

q-17.0  aT3-0  (1) 

A  correlation  in  terms  of  fluid  physical  properties  and 
tube  geometry  is  given  by  Chongrungreong  and  Sauer  [Ref.  7] . 
For  R-124  at  177  kPa,  this  reduces  to: 

q=45 . 8  aT2-32  (2) 

Another  simpler  correlation  based  only  on  system  pressure 
is  also  given  by  Chongrungreong  and  Sauer  [Ref.  7] : 

q=76 . 9  aT2-22  (3) 

Stephan  and  Abdelsalem  [Ref.  6]  present  a  statistical 
correlation  for  four  different  groups  of  fluids.  For 
refrigerants,  their  correlation  is  presented  in  the  form: 

q=  ( CIaT)  3 ' 92  (4) 

where  Cl  is  given  as  a  function  of  refrigerant  type  and  system 
pressure.  Unfortunately,  as  R-124  is  a  relatively  new 
refrigerant,  their  graphical  representation  of  Cl  does  not 
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include  R-124.  However,  Che  present  data  could  be  fitted  to  an 
equation  of  the  form  given  in  equation  (4) : 

q= ( 1 . 34  AT)3-92  (5) 

Hence,  at  a  pressure  of  1  atm,  Cl=  1.34. 

Figure  5.4  shows  the  present  R-124  smooth  tube  data  for 
decreasing  heat  flux  along  with  the  above  correlations.  At 
lower  heat  fluxes  the  prediction  of  the  Chongrungreong  and 
Sauer  [Ref.  7]  correlation  (equation  (2)  above)  exhibits 
considerably  closer  prediction  than  the  other  correlations.  In 
the  more  complex  nucleate  boiling  region,  the  Stephan  and 
Abdelsalem  [Ref.  6]  correlation  and  the  Chongrungreong  and 
Sauer  [Ref.  7]  correlation  (equation  (3))  display  a  fair 
prediction  of  the  slope. 

C.  SMOOTH  TUBE  DATA  IN  REFRIGERANT/ OIL  MIXTURES 

Although  tests  were  conducted  with  oil  concentrations  of 
0,1, 2, 3, 6,  &  10%,  only  data  with  0,3  and  10%  have  been  plotted 
in  most  of  the  figures  for  clarity.  Figures  5.5  and  5.6  show 
the  heat  transfer  performance  of  pure  R-124,  R- 124/3%  oil,  and 
R- 124/10%  oil  for  a  smooth  tube  during  increasing  and 
decreasing  heat  flux  respectively.  It  can  be  seen  that  at  the 
highest  heat  fluxes,  as  oil  concentration  is  increased  to  3%, 
the  heat  transfer  is  enhanced  by  a  maximum  of  17%  over  pure  R- 
124.  For  further  increases  in  oil  concentrations  to  10%, 
performance  drops  off.  This  drop-off  in  thermal  performance 
has  also  been  observed  by  Jensen  and  Jackman  [Ref.  13]  who 
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Figure  5.4  Comparison  of  Smooth  Tube  Performance  for 
Pure  R- 124  with  Prediction 
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Wall  Suparhaat  (K) 


Figure  5.5  Performance  Comparison  for  Boiling  R-124/  0%,  3%, 
6  10%  Mixtures  for  Smooth  Tube  (increasing  flux) 


47 


Figure  5.6  Performance  Comparison  for  Boiling  R-124/  0%,  3%, 
6  10%  Mixtures  for  Smooth  Tube  (decreasing  flux) 
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attributed  the  decline  to  a  steeper  oil  concentration 
gradient,  inhibiting  bubble  growth. 

When  oil  is  mixed  with  R-124,  two  competing  mechanisms 
occur:  (1)  Oil  accumulates  in  the  vicinity  of  the  vapor  bubble 
surface  as  the  more  volatile  refrigerant  is  evaporated; 
(2)  diffusion  of  this  accumulated  oil  back  into  the  bulk 
mixture  establishes  an  equilibrium  concentration  gradient 
near  the  departing  bubble.  Mechanisms  (1)  and  (2)  should  act 
to  inhibit  bubble  growth  and  lead  to  an  overall  decrease 
in  the  heat  transfer  with  oil  addition  compared  with  pure 
refrigerant.  It  therefore  appears  that  there  is  an  additional 
mechanism  that  has  a  significant  impact  on  the  thermal 
performance  of  pool  boiling  of  refrigerants  with  small 
quantities  of  oil.  For  R- 124/oil  mixtures  of  3%  and  greater, 
significant  foaming  in  the  evaporator  was  observed  at  mid  to 
high  heat  fluxes  (>  10000  kW/m2) .  It  is  possible  that  a  third 
mechanism  which  could  be  affecting  heat  transfer  is  this 
foaming  that  occurred  within  the  bulk  R- 124/oil  mixture. 
Udomboresuwan  and  Mesler  [Ref  .14]  reported  enhancement  in  pool 
boiling  in  the  presence  of  foam.  They  attributed  this 
enhancement  to  the  foam  causing  the  liquid/vapor  interface  to 
be  closer  to  the  tube  wall.  This  would  create  a  thin  liquid 
film  between  a  vapor  bubble  and  the  heated  tube  surface, 
thereby  providing  large  heat  -  transfer  coefficients.  Higher 
heat  fluxes  and  greater  oil  concentrations  tend  to  increase 
the  foaming,  which  could  explain  why  this  enhancement  is  seen 


49 


more  predominately  at  high  heat  fluxes.  However,  as  oil 
concentration  increases,  eventually  the  oil  concentration 
gradient  becomes  so  steep  that  mechanisms  (1)  and  (2)  tend  to 
dominate  and  heat  transfer  performance  drops  off  (ie.  the 
foaming  enhancement  effect  is  offset) . 

At  low  heat  fluxes,  the  scatter  in  the  data  is  attributed 
to  uncertainty.  For  the  increasing  heat  flux  data  (Fig  5.5), 
it  appears  that  the  addition  of  oil  tends  to  delay  the  point 
at  which  boiling  commences.  Due  to  the  statistical  nature  of 
incipience,  however,  one  should  be  wary  of  drawing  much 
conclusion  from  this  figure.  The  effect  of  oil  on  incipience 
is  discussed  in  greater  detail  later  in  this  chapter. 

Figure  5.7  compares  smooth  tube  decreasing  heat  flux  data 
from  this  thesis  for  R-124  with  Sugiyama  [Ref.  2]  R-114  data 
at  oil  concentrations  of  0%,  3%,  and  10%.  At  all  levels  of 
heat  flux,  the  R-124  data  consistently  demonstrate  improvement 
in  heat  transfer  over  R-114.  This  ranges  from  50%  at  low  heat 
flux  to  15%  at  high  heat  flux.  The  heat  flux  from  a  boiling 
surface  is  proportional  to  bubb1 3  diameter  (db)  ,  bubble 
departure  frequency  (fb) ,  number  of  active  nucleation  sites 
per  unit  area,  vapor  density  (p)  ,  and  latent  heat  of 
vaporization  (hfg)  .  Barthau  [Ref.  15]  concluded,  after 
determination  of  db,  fb,  and  active  nucleation  site  density 
for  a  smooth  copper  tube  in  R-114  using  optical  techniques, 
that  as  saturation  pressure  was  increased,  db  decreased,  while 
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rigoM  3.7  Comparison  of  HCFC-124  and  HCFC-124/Oil  Mixtures 
to  CFC-114  and  CFC-114/OiI  Mixtures  for  Smooth  Tube 


fb>  nucleation  site  density  and  heat  transfer  increased.  While 
conducting  the  tests  at  2.2  °C  the  corresponding  saturation 
pressure  of  11  psig  for  R-124  is  greater  than  the  pressure  of 
0  psig  for  R-114.  This  suggests  that  the  nucleation  site 
density  was  greater  for  R-124  than  for  R-114  at  the  same 
saturation  temperature.  Added  to  this  is  the  higher  latent 
heat  of  vaporization,  vapor  density,  and  thermal  conductivity 
for  R-124  at  this  saturation  temperature.  The  R-114/R-124 
property  comparison  at  2.2  °C  saturation  temperature  is 
summarized  in  Table  1.  It  is  therefore  reasonable  to  expect  a 
higher  heat  transfer  coefficient  for  R-124  at  the  tested 
saturation  temperature  of  2.2  °C. 

D.  BOILING  PERFORMANCE  OF  19  AND  26  FPI  GEWA-K  TUBES 

Figures  5.8  and  5.9  show  the  heat  transfer  performance  of 
a  19  fpi  GEWA-K  tube  in  pure  R-124,  R- 124/3%  oil,  and  R- 
124/10%  oil  during  increasing  and  decreasing  heat  flux 
respectively.  There  is  a  consistent  increase  in  thermal 
performance  for  oil  concentrations  of  3%  and  10%  over  pure  R- 
124  at  high  heat  fluxes.  This  increase  is  substantially 
greater  than  that  obtained  with  the  GEWA-K  26  fpi  tube  for  the 
same  oil  mixtures,  as  shown  in  Figures  5.10  and  5.11  for 
increasing  and  decreasing  flux  respectively.  Surprisingly  for 
both  tubes,  there  seems  to  be  a  significant  influence  of  oil 
in  the  natural  convection  region  (Figs  5.8  &  5.10).  There 
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Figure  5.8  Performance  Comparison  for  Boiling  R-124/  0%,3%,& 
10%  Mixtures  for  19  fpi  GEWA-K  Tube  (increasing  flux) 
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.9  Performance  Comparison  for  Boiling  R- 124/0%,  3%,  10% 
Mixtures  for  19  fpi  GEWA-K  Tube  (decreasing  flux) 
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&  10%  Oil  Mixtures  for  26  fpl  GEWA-K  Tube  (Increasing  flux) 
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Will  Suptrhett  (K) 


Figure  5.11  Performance  Comparison  for  Boiling  R-124/  0%,3%, 
&  10%  Oil  Mixtures  for  26  fpi  GEWA-K  Tube  (decreasing  flux) 


tends  to  be  a  decrease  in  performance  for  increasing  oil 
concentration.  This  was  not  so  obvious  with  Sugiyama  [Ref.  2] 
R- 114  data.  Once  again,  the  point  of  incipience  seems  to  be 
slightly  delayed  with  increasing  oil  percentage.  The  thermal 
performance  of  the  finned  tube  is  better  than  the  smooth  tube, 
evidenced  by  thermal  enhancements  of  2.1  to  2.9  with  the  19 
fpi  tube  and  1.3  to  2.1  with  the  26  fpi  tube.  However,  it 
should  be  noted  that  the  heat  flux  calculation  does  not 
account  for  the  increased  surface  area  of  the  finned  tube,  but 
rather  uses  a  root  diameter.  Some  of  the  enhancement  ic 
therefore  due  to  the  increased  surface  area  of  the  fins  (Figs. 
5.24-5.26  compare  heat  transfer  coefficients  with  smooth 
tube)  ,  with  each  fin  behaving  essentially  as  the  copper 
surface  of  a  smooth  tube  in  terms  of  thermal  performance. 

Figures  5.12  and  5.13  compare  R-124  and  R-114  decreasing 
heat  flux  data  for  the  19  and  26  fpi  GEWA-K  tubes  respectively 
for  0%,  3%,  and  10%  oil  concentrations.  The  19  fpi  tube  shows 
considerable  improvement  over  the  entire  heat  flux  range  while 
the  26  fpi  tube  shows  improvement  only  in  the  low  heat  flux 
region.  The  improvement  of  the  lower  fin  density  tube  can  be 
speculated  to  be  due  in  large  part  to  better  circulation  and 
mixing  between  the  fins.  For  both  tubes  at  low  oil 
concentrations  and  high  heat  flux,  the  foaming  effect  provides 
improvement  in  heat  transfer  similar  to  that  seen  with  the 
smooth  tube. 
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Figure  5.12  Comparison  of  IICIC-124  and  HCFC-124/Oil  Mixtures 
to  CFC-114  and  CFC-114/Oil  Mixtures  for  19  fpi  Tube 


Figure  5.13  Comparison  of  HCFC-124  and  HCFC-124/OU  Mixtures 
to  CFC-114  and  CFC-114/Oil  Mixtures  for  26  fpi  Tube 


E.  BOILING  PERFORMANCE  OF  HIGH  FLUX  AND  TURBO-B  TUBES 

Figures  5.14  and  5.15  show  the  heat  transfer  performance 
of  a  HIGH  FLUX  tube  in  pure  R-124,  R-124/3%  oil,  and  R-124/10% 
oil  during  increasing  and  decreasing  heat  flux  respectively. 
Uncertainty  analysis  (Appendix  D)  was  performed  on  the  HIGH 
FLUX  tube  and  typical  uncertainty  bands  are  also  given.  The 
uncertainty  is  largest  at  low  wall  superheat  and  low  heat  flux 
since  the  uncertainty  is  of  the  order  of  the  measured  values. 
It  can  be  seen  that  the  heat  transfer  performance  decreases 
for  all  heat  fluxes  as  oil  concentration  increases.  This 
decrease  in  performance  is  due  most  likely  to  the  oil  becoming 
trapped  within  the  porous  coating.  As  the  refrigerant  is 
evaporated,  a  very  steep  concentration  gradient  is  developed 
within  the  pores;  oil  diffusion  is  also  retarded  because  of 
the  presence  of  the  interconnected  pores;  this  is  particularly 
severe  at  high  heat  fluxes  with  10%  oil  mixture,  where  a  very 
large  drop-off  in  performance  is  seen.  In  the  natural 
convection  region,  the  effect  of  oil  is  not  as  pronounced  as 
at  high  flux,  exhibiting  only  minor  reduction  in  heat  transfer 
compared  to  the  pure  refrigerant.  Any  advantageous  effect  of 
foaming  is  negated  due  to  the  pores  remaining  impenetrable  to 
the  foam.  Figure  5.16  compares  R-124  and  R-114  decreasing  heat 
flux  data  for  the  HIGH  FLUX  tube  for  0%,3%,  and  10%  oil 
concentrations.  The  trend  of  the  data  for  both  fluids  is  very 
similar.  With  such  a  high  density  of  nucleation  sites  already 
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Figure  5.15  Performance  Comparison  for  Boiling  R- 124/0%, 3% 
&  10%  Oil  Mixtures  for  HIGH  FLUX  Tube  (decreasing  flux) 
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Figure  5.16  Comparison  of  HCFC424  and  HCFC*124/Oil  Mixtures 
to  CFC-114  and  CFC-114/Oil  Mixtures  for  HIGH  FLUX  Tube 


present  on  the  surface,  there  does  not  appear  to  be  a 
significant  thermal  enhancement  at  higher  operating  pressures 
as  seen  with  the  finned  and  smooth  tubes. 

Figures  5.17  and  5.18  show  the  heat  transfer  performance 
of  a  TURBO-B  tube  in  pure  R-124,  R-124/3%  oil,  and  R-124/10% 
oil  mixtures  during  increasing  and  decreasing  heat  flux 
respectively.  The  results  of  oil  addition  are  similar  to  the 
HIGH  FLUX  tube.  The  notable  exception  is  that  at  high  oil 
concentrations  and  high  heat  fluxes,  the  heat  transfer 
performance  of  the  TURBO-B  tube  does  not  drop  off  as  rapidly 
as  the  HIGH  FLUX  tube.  Again,  only  a  small  effect  of  oil  in 
the  natural  convection  region  was  observed. 

Figure  5.19  compares  R-124  and  R-114  decreasing  heat  flux 
data  for  the  TURBO  B  tube  for  0%,  3%,  and  10%  oil 

concentrations.  Again,  this  comparison  is  similar  to  the 
results  already  discussed  with  the  HIGH  FLUX  tube. 

F.  EFFECT  OF  R- 124/OIL  MIXTURES  ON  INCIPIENCE  OF  NUCLEATE 

BOILING 

Separate  tests  were  conducted  to  analyze  the  incipience  of 
nucleate  boiling  or  the  onset  of  nucleate  boiling  (ONB)  in 
pure  R-124  and  how  this  was  affected  by  the  presence  of  oil. 
These  tests  were  taken  in  accordance  with  the  general 
procedures  outlined  in  the  preceeding  chapter.  However, 
smaller  heat  flux  increments  were  applied  to  the  test  tube 
during  the  natural  convection  region  in  order  to  obtain  the 
most  accurate  wall  superheat  at  the  ONB.  Also,  during  these 
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.17  Performance  Comparison  for  Boiling  R-124/0%,3% 
&  10%  Oil  Mixtures  for  TURBO-B  Tube  (increasing  flux) 
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Figure  5.18  Performance  Comparison  for  Boiling  R-124/0%,3%, 
&  10%  Oil  Mixtures  for  TURBO-B  Tube  (decreasing  flux) 
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ngur*  5.19  Comparison  of  HCFC-124  and  HCFC-124/Oil  Mixtures 
to  CFC-114  and  CFC-114/Oil  Mixtures  for  TURBO-B  Tube 


incipience  tests  the  refrigerant  pool  was  cooled  down  at 
approximately  half  the  rate  of  the  prior  increasing  and 
decreasing  heat  flux  runs.  Each  incipience  test  was  conducted 
with  a  minimum  of  8  hours  between  between  re- testing  to  allow 
adequate  time  for  the  R-124  pool  to  settle  and  for  the 
nucleation  sites  to  de-activate.  The  maximum  time  allowed  for 
the  aging  of  the  tube  surface  was  48  hours.  Similar  tests  were 
repeated  seven  times  on  the  smooth  and  HIGH  FLUX  tubes.  Figure 
5.20  summarizes  the  data  on  a  probability  chart.  Due  to  the 
statistical  nature  of  ONB,  You  et  al  [Ref.  16]  proposed  this 
chart  as  the  best  means  of  presenting  incipient  data.  The 
seven  incipient  wall  superheats  for  each  condition  are  given 
for  the  smooth  and  HIGH  FLUX  tube  at  0%,  3%,  and  10%  oil 
concentrations  as  the  percent  probability  of  nucleation 
occurring.  For  example,  the  smooth  tube  in  pure  refrigerant 
never  nucleated  at  a  wall  superheat  less  than  12°C.  Also  the 
data  show  that  for  a  smooth  tube  in  pure  R-124,  nucleation 
occurred  four  times  at  a  wall  superheat  greater  than  20  °C, 
and  three  times  at  a  lower  wall  superheat;  this  corresponds  to 
a  probability  of  nucleation  occurring  at  a  superheat  of  20  °C 
between  3/7  (43)%  and  4/7  (57)%.  As  the  superheat  further 
increases,  the  probability  increases  also  until  at  a  superheat 
of  28°C  there  is  a  100%  probability  of  nucleation.  Bar-Cohen 
[Ref.  3]  also  reported  widely  varying  R-113  superheats  when 
'departure  from  convection'  or  nucleation  occurred.  The 
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Figure  5.2  0  Probability  of  ITucleation  Occurring  with  Smooth 
and  HIGH  FLUX  Tube  for  0%,3%,  &  10%  Oil  Mixture 
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probability  chart  of  Figure  5.20  can  be  made  more  accurate  by 
plotting  a  greater  number  of  incipient  superheats  for  each 
tube . 

For  the  smooth  tube,  there  is  a  distinct  shift  to  higher 
superheats  with  increases  in  oil  concentration.  This 
represents  a  delay  in  incipience  caused  by  the  addition  of 
oil.  Memory  and  Marto  [Ref.  10]  conducted  similar  tests  with 
R-114  and  R- 114/oil  mixtures  and  reported  the  oil 
concentration  to  have  no  systematic  effect  on  the  point  of 
incipient  boiling  for  the  smooth  tube.  Wanniarachchi  et  al 
[Ref.  9]  reported  an  increase  in  the  wall  superheat  needed  to 
initiate  boiling  for  the  smooth  tube  in  R-114  with  10%  oil, 
noting  that  3%  oil  caused  no  increase.  Memory  and  Marto  [Ref. 
10]  precicted  the  effect  of  oil  on  superheat  needed  to 
nucleate  from  a  vapor- filled  cavity  by  incorporating  the 
physical  properties  of  the  refrigerant/oil  mixture  into  a 
nucleation  parameter  (N)  given  r  '  Marto  and  Lepere  [Ref.  17] : 

"  a  '  Tsat/ ^vapor  '  g 

If  the  parameter  were  to  increase,  a  delay  in  the  incipient 
point  with  increasing  oil  concentrations  would  appear  likely. 
However,  it  is  not  clear  how  the  oil  effects  this  parameter; 
with  the  smooth  tube  in  R- 124/oil  mixtures  it  resulted  in  a 
delay  in  incipience,  while  the  HIGH  FLUX  tube  exhibited  no 
distinct  delay  in  incipience.  Also  the  HIGH  FLUX  tube  has  a 
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whole  range  of  sites  that  can  be  easily  activated  while  the 
smooth  tube  only  tends  to  have  very  small  sites  which  may  move 
outside  the  range  of  sites  that  can  be  activated  when  a 
refrigerant/oil  mixture  is  used. 

In  the  case  of  pure  refrigerant,  it  does  appear  that  with 
the  relatively  few  incipient  tests  conducted,  the  wall 
superheats  at  incipience  for  the  smooth  and  HIGH  FLUX  tubes 
are  larger  on  average  for  R-114  than  for  R-124.  This  initial 
result  would  support  the  approximately  50%  greater  nucleation 
parameter  of  R-114  over  R-124.  This  subject  requires  further 
study. 

Figure  5.21  shows  an  expanded  view  of  the  smooth  tube's 
heater  sleeve  rolled  out  to  show  the  exact  locations  of  each 
thermocouple  (longitudinal  and  circumferential)  imbedded  on 
the  surface  of  the  sleeve.  Several  data  sets  were  taken  for 
the  smooth  tube  and  then  repeated  after  reversing  the  position 
of  the  tube  within  the  pool.  This  procedure  was  conducted  to 
see  if  there  were  any  noticeable  trends  in  the  thermocouple 
readings  which  may  be  an  indication  of  local  hot  or  cold 
regions  within  the  pool.  However,  this  reversal  of  the  tube 
had  no  consistent  effect  on  the  thermocouple  readings  or  the 
point  at  which  nucleation  occurred.  This  appears  to  underline 
the  uniformity  of  thermocouple  output  within  the  pool.  A  more 
detailed  survey  of  this  data  needs  to  be  done. 
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Left  Side 
of  Evaporator 


Right  Side 
of  Bvaporator 


5.21  Expanded  View  of  Thermocouple  Locations 
the  Heater  Sleeve  of  Smooth  Tube  (#1  at  12  0' 


within 

clock) 


G.  SUMMARY  OF  THE  R-124  DATA 


Figures  5.22  and  5.23  summarize  the  effect  of  all  oil 
concentrations  ( 0 , 1 , 2 , 3 , 6 , &10% )  in  R-124  for  the  five  tubes  at 
heat  fluxes  of  25  and  90  KW/m2  respectively;  2  5  kW/m2 
represents  a  typical  operating  heat  flux.  For  each  tube,  the 
ratio  of  the  heat- transfer  coefficient  at  a  given  oil 
concentration  divided  by  that  for  the  same  tube  in  pure  R-124 
is  given.  At  low  heat  flux,  the  smooth  tube  and  the  19  fpi 
tube  give  an  increase  in  performance  as  oil  is  added,  up  to 
concentrations  of  6%.  The  smooth  tube  does  the  best  in  this 
respect,  giving  an  increase  in  heat  transfer  of  about  23% 
(compared  to  a  smooth  tube  in  pure  R-124)  at  6%  oil 
concentration.  The  HIGH  FLUX  and  TURBO -B  tubes  show  a  decrease 
in  performance  as  oil  is  added,  while  the  26  fpi  tube  shows  no 
change.  At  the  highest  oil  concentration  (10%) ,  performance 
for  all  tubes  is  worse  than  that  at  6%,  although  the  smooth 
tube  still  indicates  better  heat  transfer  than  when  no  oil  was 
added,  which  is  rather  surprising.  Similar  trends  were 
obtained  for  R-114,  although  at  10%  oil,  all  tubes  had  poorer 
heat  transfer  than  with  no  oil.  At  higher  heat  flux,  general 
effects  are  similar,  although  now  the  19  fpi  tube  shows  the 
best  performance  of  about  35%  at  6%  oil  concentration. 
Interestingly,  between  3  and  6%  oil  concentrations,  all  tubes 
except  the  HIGH  FLUX  tube  show  improved  heat  transfer  over  the 
same  tube  in  pure  refrigerant.  Although  below  3%  oil  there  is 
an  unexplained  'dip'  in  the  heat  transfer  coefficient  ratio 
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with  the  TURBO- B  tube.  The  HIGH  FLUX  surface  has  very  poor 
performance  at  10%  oil  concentration,  being  only  35%  of  that 
in  the  pure  fluid. 

Figures  5.24,  5.25  and  5.26  compare  the  five  tubes 
directly  for  oil  concentrations  of  0,  3,  and  10%  respectively. 
From  these  figures,  enhancements  (ie.  improvement  in  heat 
transfer  compared  to  the  smooth  tube)  can  be  calculated  at  any 
desired  heat  flux.  These  are  provided  in  Table  4  at  heat 
fluxes  25  and  90  kW/m2 .  The  advantages  of  enhanced  surfaces 
are  clearly  apparent.  The  smooth  tube  always  has  the  poorest 
performance  (as  expected)  and  is  used  as  a  base-line  with 
which  the  other  tubes  are  compared.  Between  the  finned  tubes, 
the  19  fpi  tube  consistently  outperforms  the  26  fpi  tube  at 
all  heat  fluxes  and  oil  concentrations:  enhancements  range 
between  2. 1-2. 9  for  the  19  fpi  tube  and  1.3 -2.1  for  the  26  fpi 
tube.  Indeed,  at  the  highest  heat  fluxes  at  10%  oil 
concentration,  the  19  fpi  tube  is  the  best  tube  tested 
(enhancement  of  2.6).  The  GEWA-K  19  fpi  tube  to  smooth  tube 
area  ratio  is  2.4,  while  the  GEWA-K  26  fpi  tube  to  smooth  tube 
area  ratio  is  3.2.  The  HIGH  FLUX  and  TURBO- B  tubes  are  the 
best  tubes  over  most  conditions  of  heat  flux  and  oil 
concentration,  with  both  having  fairly  equal  performance.  The 
exception  is  the  HIGH  FLUX  tube  at  high  heat  flux  and  oil 
concentration,  where  performance  drops  off  rapidly,  as  already 
mentioned.  Over  practical  heat  flux  ranges  (15-30  kW/m2) 
therefore,  the  porous  coated  and  modified  fin  tubes  are  the 
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best  for  all  oil  concentrations.  Although  the  TURBO- B  and  HIGH 
FLUX  tube  performance  is  degraded  at  10%  oil  concentration, 
in-service  oil  separators  should  maintain  a  3%  or  less 
concentration  in  an  actual  evaporator  system.  At  this  level, 
the  TURBO -B  and  HIGH  FLUX  tube  are  superior  at  normal  chiller 
heat  flux  between  10  and  30  KW/m2 . 

A  point  of  interest  worth  noting  was  that  the  Viton  and 
neoprene  0- rings  used  experienced  considerable  swelling  of 
approximately  30%  by  volume  when  subjected  to  R-124.  This  was 
enough  to  prevent  their  re-use  when  the  Teflon  inserts  were 
removed  to  change  the  tube.  Also  noted  was  that  the  0- rings 
only  swelled  when  taken  out  of  the  system.  When  undisturbed, 
they  maintained  a  good  seal  {this  might  not  be  true  if  the 
tube  were  sliding  past  the  0-ring) .  Also  significant,  after 
removed  from  use  within  the  apparatus  for  60  days,  both  types 
of  0- rings  did  shrink  back  to  their  original  size. 


Table  4.  ENHANCEMENT  SUMMARY  FOR  R-124 


Heat  Flux 

25  kW/m2 

90  kW/m2 

Oil  Concentration 

0% 

3% 

10% 

0% 

3% 

B 

Smooth 

1.0 

1 . 0 

1.0 

1.0 

1 . 0 

1 . 0 

19  fpi 

2.9 

2.7 

2.8 

2 . 1 

2.4 

2.6 

26  fpi 

| 

IB 

1 . 8 

IB 

IB 

1 . 5 

HIGH  FLUX 

6.2 

4 . 9 

4 . 8 

2.9 

2.3 

B 

TURBO - B 

6.0 

4 . 3 

B 

2.4 

2.2 

1 . 9 
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Flgur*  5.22  Effect  of  Oil  Concentration  on  Ileal  Transfer 
Coefficient  at  a  Heat  Flux  of  25  kW/m* 


OI'ZvU»/M)/»u»PUJ»03 


78 


Hc«t  Flu>/(W/mA2) 
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5.26  Coapariaon  of  All  Tub**  *t  10§  oil 
Concentration  in  R-124 


VI .  CONCLUSIONS 


1.  For  the  smooth  and  finned  tubes,  significant  improvements 
in  heat  transfer  (>50%)  were  obtained  in  pure  R-124 
compared  with  pure  R-114  at  the  same  saturation 
temperature.  This  was  attributed  to  the  higher  saturation 
pressure  of  R-124  activating  more  nucleation  sites  and  the 
higher  latent  heat  of  vaporization,  vapor  density,  and 
thermal  conductivity  of  R-124  over  R-114. 

2.  For  the  HIGH  FLUX  and  TURBO -B  surfaces,  performance  in  the 
two  refrigerants  was  similar,  since  the  number  of  active 
nucleation  sites  is  already  high. 

3.  For  pure  R-124,  the  HIGH  FLUX  tube  gave  the  largest 
enhancement  of  any  tube  over  the  whole  range  of  heat  flux. 

4.  For  the  smooth  and  finned  tubes,  significant  improvements 
in  the  heat  transfer  (>50%)  were  also  obtained  in  R- 
124/oil  mixtures  when  compared  with  R- 114/oil  mixtures  at 
the  same  saturation  temperature.  This  is  also  assumed  to 
be  due  to  the  higher  saturation  pressures  associated  with 
R-124. 

5.  Increasing  the  oil  concentration  in  R-124  by  up  to  6%  for 
the  smooth  and  finned  tubes  improved  heat  transfer 
performance  compared  to  the  pure  refrigerant.  For  example 
at  a  heat  flux  of  at  90  kW/nr ,  performance  for  the  smooth 
and  finned  tubes  increased  by  18%  and  35%  respectively.  At 
higher  oil  concentrations,  performance  deteriorated. 

6.  Any  increase  in  oil  concentration  in  R-124  for  the  HIGH 
FLUX  or  TURBO-B  tubes  led  to  a  decrease  in  performance 
such  that  at  high  heat  fluxes  and  high  oil  concentrations, 
the  HIGH  FLUX  tube  gave  only  35%  of  its  performance  in 
pure  R-124  and  the  19  fpi  tube  gave  the  best  overall 
performance . 

7.  At  chiller  design  heat  fluxes  (10-30  kW/m2)  and  oil 
concentrations  of  3%  maximum,  the  HIGH  FLUX  and  TURBO-B 
surfaces  still  provide  the  best  heat  transfer  which  is 
about  the  same  as  with  R-114.  Therefore,  similar  or 
slightly  better  performance  is  expected  with  R-124  when 
used  as  a  'drop-in'  replacement  for  R-114  in  a  refrigerant 
evaporator  operated  at  the  same  saturation  temperature. 

8.  Addition  of  oil  delays  the  point  of  incipience  (occurs  at 
a  higher  wall  superheat)  for  the  smooth  tube.  For  the  HIGH 
FLUX  tube,  oil  addition  has  no  measurable  effect  on  the 
wall  superheat  at  incipience. 
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VII.  RECOMMENDATIONS 


1.  The  remaining  enhanced  tubes  should  be  tested  in  R-124  at 
0%  ,  3% ,  and  10%  oil  concentrations  and  compared  to  the  five 
tubes  above  and  their  prior  performance  in  R-114  and  R- 
114/oil  mixtures. 

2.  Incipience  tests  performed  with  finned  tubes  and  TURBO-B 
tubes  should  be  conducted  to  compare  the  influence  of  oil 
on  wall  superheat. 

3.  The  influence  of  the  variation  in  time  between  test  runs 
on  boiling  nucleation  (ie.  onset  of  nucleate  boiling)  must 
be  investigated  more  thoroughly. 

4.  The  single  tube  apparatus  should  be  modified  to  test  an 
array  of  two  horizontal  boiling  tubes  and  compare  these 
results  with  the  R-114  tube  array  data  previously 
recorded. 

5.  A  prototype  chiller  presently  using  R-114  should  be  tested 
with  R-124  to  monitor  the  performance. 
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APPENDIX  A.  THERMOPHYSICAL  PROPERTIES  OF  R-124 


The  following  thermophysical  properties  of  saturated  R-124 
are  plotted  versus  temperature  in  degrees  Celsius.  They  were 
generated  from  REFPROP  [Ref.  12]. 


1 . 

Saturation  Pressure  /(kPa) 

2  . 

Liquid  Density  /(kg/m3) 

3  . 

Latent  Heat  of  Vaporization  /(kJ/kg) 

4  . 

Liquid  Specific  Heat  / (kJ/kg  K) 

5. 

Liquid  Viscosity  /(mP)  note:  1  mP=  lxlO-4  kg/m  s 

6. 

Liquid  Thermal  Conductivity  / (W/m  K) 
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Saturation  Pressure  for  R-124 


(BdiO/ssajd 


in 

+ 

CO 

VO 

« 

<N 

VO 

iH 

II 

4J 

10 

(0 

Q. 
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.94  6(T)+9.2081X10'2(T)2+6.9359X10  4(T) 


Liquid  Density  for  R-124 


(£W/3>i)/otiH 


n 

o. 
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1434.8-2.8619  (T)  -6 . 7267xlO‘3(T)2-7 . 2852xlO'5(T) 


Latent  Heat  of  Vaporization  for  R-124 
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159 .93-.  4609  (T)  -1 . 458xlO'5(T)2-l .  3715x10  5 (T) 


(M  8*/f>1)/<l3 


Liquid  Viscosity  for  R-124 


p»  »o  *r»  ■<*  m  CN 


(diuJ/asjA 


26xl0  2(T)+3.9458xl0  *(T)2-4. 193X10  6(T)3+2.07xl0'8(T) 


Thermal  Conductivity 


5191x10  -3 . 54  36x10  *  (T)  -1 . 9545xlO'7(T)  2+3 . 1835xlO‘9(T) 


APPENDIX  B.  REPRESENTATIVE  DATA  SET 


Oate  :  19  Nov  1992 

NOTE:  Program  name  :  DRPS8 
Olsk  number  »  01 

Oate  :  1 9  Nov  1992 

NOTE:  Program  name  :  ORPGB 
Disk  number  ■  00 
Neu  file  name:  01 119193 
TC  are  defective  at  locations  7  8 

Tube  Number:  13 


Data  Set  Number  •  I  Bulk  Oil  X  •  3.0 

TIME:  15:40:01 

TC  No:  I  2  3  4  S  E  7  8 

Temp  :  3.23  3.38  3.31  3.20  3.0S  3.32  -99.99  -99.99 

Twa  Tliod  Tliqd2  Tvapr  Psat  Tsump 

3.25  2.14  2.14  2.18  10.84  -13.3 

Thetab  Htube  Odp 

I.0G3  B.47SEv02  6.879E+02 


Oata  Sat 
TINE: 

TC  No: 

Number  •  2 

15=42:22 

1  2 

Bulk  Oil  X  * 

3  4 

3.0 

S 

6  7 

Temp  : 

3.58 

3.80 

3.72  3.61 

3.45 

3.64  -99.99 

T  w  a 

3.63 

Tliqd 

2.22 

Tliqd2 

2.24 

Tvapr  Psat 

2 . 25  10.91 

Tsump 

-13.2 

Thetab  Htube  Odp 

1.377  6 . 422E*02  8.84SE+02 


Oata  Set  Number  •  3  Bulk  Oil  X  »  3.0 

TINE:  15:45:20 

TC  No:  I  2  3  4  5  E  7  8 

Temp  :  3.85  4.35  4 . 35  4.19  3.87  3.95  -99.99  -99.99 

T«a  Tliqd  TliqdZ  Tvapr  Psat  Tsump 

4.09  2.09  2.21  2.17  10.83  -13. 1 

Thetab  Htube  Qdp 

1.921  6 . 523E»02  1.2S3E*03 


Oata  Set 

Number 

*  4  Bulk  Oil  X  - 

3.0 

TINE: 

15:47: 

:  12 

TC  No : 

1 

2  3  4 

6  6 

Te«p  : 

4.  10 

S.03  5.02  4.87 

4.46  4.28 

T  wa 

Tliqd 

Tliqd2  Tvapr  Psat 

T  sump 

4.62 

2.15 

2.27  2.27  10.92 

-13.0 

Thetab 

Htube 

Odp 

2.350 

6.986£< 

■02  1  . 642E+03 
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Oata  Set 

Number  “  5  Sulk  Oil  X* 

3.0 

TIN£: 

15:49:09 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

T  e*>p  ; 

4.30  5 . 75  5. 76 

5.44 

4.68 

4.56 

-39.99  -99.99 

Tua 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sump 

5 . 07 

2.05  2.20  2.23 

10.38 

-12.9 

Thetab 

Htube  Qdp 

2.345 

7. 8326*02  2.2236*03 

Cata  Set 

Number  *  6  Sulk  Oil  X  * 

3.0 

TIN£: 

15:51:30 

TC  No: 

1  2  3 

4 

5 

6 

7  3 

T  ere  : 

4.40  6.62  6.60 

6.21 

5.12 

4.73 

-99.93  -99.33 

T  ua 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sump 

5.61 

2.12  2. IS  2.28 

10.93 

-12.8 

Thetab 

Htube  Qdp 

3.331 

9.1546*02  3.0496*03 

Oata  Set 

Number  »  7  Sulk  Oil  X  ■ 

3.0 

TIME: 

15:54:01 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Temp  : 

4.75  7.39  7.27 

6.75 

5.38 

5.22 

-99.99  -99.99 

T  ua 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sump 

6.11 

2.15  2.21  2.28 

10.93 

-12.7 

Thetab 

Htube  Qdp 

3.831 

1.0316*03  3.9516*03 

Oata  Set 

Number  •  8  Bulk  Oil  X  * 

3.0 

TIME: 

15:56:43 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Temp  : 

5.14  8 . 28  8.15 

7.38 

5.55 

5.59 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sump 

6.66 

2.24  2.27  2.3S 

11.00 

-12.6 

Thetab 

Htube  Qdp 

4.302 

1.2596*03  5.4186*03 

Oata  Set 

Number  ■  9  SulkOilX* 

3.0 

TIME: 

15:59:50 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Temp  : 

4.78  6.89  8.17 

6.85 

4.88 

4.46 

-99.99  -39.99 

Tua 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sumo 

5.97 

2.29  2.24  2.32 

10.97 

-12.5 

Thetab  Htube  Qdp 

3 .550  2.1516*03  7.849E*03 

Data  Set  Number  «  10  Sulk  Oil  X  ■  3.0 


TIME: 

16:02:33 

TC  No: 

1  2 

3 

4 

5  6 

7 

8 

Temp  : 

4.83  5.18 

5.13 

4.84 

4.91  4.53 

-99.99 

-99.99 

Tua 

Tliqd  Tliqd2 

T  vapr 

Psat 

T  sump 

4.36 

2.30  2.25 

2.28 

10.93 

-12.4 

Thetab 

Htube  Qdp 

2.583 

4.1506*03  1.0726*04 

Oata  Set 

Number  ■  1  1 

Sulk  Oil  X  - 

3.0 

TIME: 

16:17:38 

TC  No: 

1  2 

3 

4 

5  6 

7 

8 

Temp  : 

5.14  5.66 

5.59 

5.22 

5.25  4.77 

-99.99 

-99.99 

T  ua 

Tliqd  TliqdC 

Tvapr 

Psat 

T  sump 

5.20 

2.20  2.12 

2.16 

10.82 

-11.9 

Thetab 

Htube  Qdp 

3.044  S . 093E  *03  1.5526*04 
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Oata  Sat 

Number  «  12  Sulk  Oil  X  • 

3.0 

TIME: 

16:20:31 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

T  erp  : 

S .79  6.34  6 . 28 

5.90 

6.01 

5.24 

-99.99  -99.99 

T  wa 

Tiiod  T 1  iqd2  Tvapr 

Psat 

T  sump 

S .  83 

2.26  2.19  2.20 

10.86 

-11.9 

Thetab 

Htube  Odp 

3.630 

6.3286*03  2.2976*04 

Data  Set 

Number  »  13  Bulk  Oil  X  * 

3.0 

TIME: 

16:23: IS 

TC  No : 

1  2  3 

4 

5 

6 

7  8 

Temp  : 

6.61  7.18  7.13 

6.69 

6.85 

5.99 

-99.99  -99.99 

Twa 

Tltqd  Tliqo2  Tvapr 

Psat 

T  sump 

6.60 

2.31  2.26  2.28 

10.93 

-11.8 

That ab 

Htube  Odp 

4.316 

7.794E+03  3. 364E*04 

Data  Sat 

Number  ■  14  8ulk  Oil  X  • 

3.0 

TIME: 

16:24:49 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Temp  : 

7.40  7.84  7 . 85 

7.36 

7.58 

6.49 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr 

Psat 

Tsump 

7.22 

2.35  2.27  2.32 

10.97 

-11.8 

Thetab 

Htube  Qdp 

4.89S 

9. 726E*03  4 . 76 1 E  *04 

Data  Sat 

NuPbar  •  IS  Sulk  Oil  X« 

3.0 

TIME: 

16:29:15 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Temp  : 

8.17  8.32  8.52 

8.1 1 

8.32 

7.09 

-99.99  -99.99 

Twa 

Tliqd  TUqd2  Tvapr 

Psat 

T  sump 

7.78 

2.29  2.23  2.28 

10.93 

-11.7 

Thetab 

Htube  Odp 

S.50I 

i.300E*04  7.1516*04 

Data  Sat 

NuPbar  t  16  Sulk  Oil  X* 

3.0 

TIME: 

16:33:55 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Terp  : 

8.79  8.96  9.33 

8.87 

9.12 

7.59 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sump 

8.37 

2.20  2.'4  2.20 

10.86 

-11.4 

Thetab 

Htuba  Qdp 

6.169 

1.5536*04  9.5776*04 

Oete  Set 

NuPbar  «  17  Sulk  Oil  X* 

3.0 

TIME: 

16:36:37 

TC  No: 

1  2  3 

4 

5 

6 

7  8 

Te«p  : 

7.44  7.74  7.92 

7.59 

7.77 

6.56 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr 

Psat 

T  sump 

7.24 

2.18  2.11  2.17 

10.83 

-11.4 

Thetab 

Htube  Odp 

S.070 

1.2306*04  6 . 235E*04 

Data  Sat 

NuPbar  •  18  Bulk  Oil  Xs 

3.0 

TIME: 

16:40:00 

TC  No: 

I  2  3 

4 

5 

6 

7  8 

Temp  : 

6.63  7.03  7.12 

6.77 

6.98 

5.96 

-99.99  -99.99 

Twa 

Tliqd  TliqdC  Tvapr 

Psat 

T  sump 

6.57 

2.24  2.18  2.24 

10.90 

-11.4 

Thetab 

Htube  Odp 

4.329 

9.657E*03  4.1806*04 
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Oata  Sat 

Number  “  '9  Bulk  Oil  X  • 

3.0 

T  :*£  : 

1 6 : *3 : 3* 

*C  Sc : 

;  3  * 

5 

5 

7  a 

T  e'i  : 

6.09  6.65  5.55  5.;" 

6.41 

5.54 

-99.99  -99.99 

T*a 

Tliqd  TliqoC  'vap<*  Psat 

7  S-»C 

6.  '* 

2.27  2.20  2.2*  10.90 

-i  i  .5 

Thetab 

HtuOe  Qdp 

3. 85-3 

7 . 8596*03  3.06'E-0* 

Data  Set 

Number  *  20  Sulk  Oil  X  » 

3.0 

Ting: 

16:45-59 

TC  No: 

1  2  3  * 

5 

6 

7  8 

Temp  : 

5.62  6.13  6.13  5.32 

5.92 

5.  19 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr  Psat 

T  sump 

5.70 

2.23  2.1*  2.19  10.35 

-11.5 

Thetab 

Htube  Qdp 

3.506 

6 . 7S4£*03  2 . 368E»0* 

Oata  Set 

Number  -  21  Sulk  Oil  X  » 

3.0 

TIME: 

16:48:34 

TC  No: 

1  2  3  * 

5 

6 

7  8 

Temp  : 

S.I7  5.73  5.71  5.36 

5.40 

4.85 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr  Psat 

T  Sump 

5.30 

2.31  2.22  2.25  10.91 

-11.7 

Thetab 

Htube  Qdp 

3 . 0*5 

S.S55E+03  1 . 692E+0* 

Oata  Set 

Number  »  22  Sulk  Oil  I  * 

3.0 

TIME: 

16:51 :59 

TC  No: 

12  3* 

5 

6 

7  8 

Temp  : 

*.69  5.20  S.IS  4.85 

4.92 

4.44 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr  Psat 

T  sump 

*.82 

2.27  2.17  2.21  10.87 

-ll.8 

Thatab 

Htube  Qdp 

2.613 

* . 52SE*03  1.I82E*04 

Oata  Set 

Number  •  23  Bulk  Oil  X  • 

3.0 

TIME: 

16:54:33 

TC  No: 

12  3* 

5 

6 

7  8 

T emp  : 

4.32  4.77  4.63  *.*5 

4.55 

4.2! 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr  Psat 

T  su"P 

«.*S 

2.30  2.21  2.27  10.92 

-11.9 

Thatab 

Htube  Qdp 

2.186 

3.6S3E-03  7. 9S5E*03 

Oata  Set 

Number  “  24  Bulk  vii  t  ' 

3.0 

TIME: 

16 :S7 :59 

TC  No: 

1  2  3  * 

5 

6 

7  8 

Temp  : 

3.94  4.31  4.23  *.I0 

4.11 

3.91 

-99.99  -99.99 

Twa 

TUqd  Tliqd2  Tvapr  Psat 

T  sump 

*.08 

2.2*  2.15  2.21  10.87 

-12.0 

Thetab 

Htube  Qdp 

1.863 

2.839E*03  5.299E»03 

Oata  Sat 

Number  «  25  Bulk  Oil  X  * 

3.0 

TIME: 

17:01:00 

TC  No: 

l  2  3  * 

5 

6 

7  8 

Temp  : 

3.56  3.91  3.34  3.75 

3.74 

3.60 

-99.99  -99.99 

Twa 

Tliqd  Tliqd2  Tvapr  Psat 

T  sump 

3.72 

2.22  2.08  2.1*  10.80 

-12.1 

Thetab 

Htube  Qdp 

1.577  2 . i SSe*03  3.41 9E*03 
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Oata  Sat  Nuxba*  -  2a  Suit  3*1  S  *  3  * 

TI«:  I7:«;:1S 


rc  No: 

12  3  4 

5  6  7  9 

T««o  : 

3.SS  3.83  3.77  j.7» 

369  3. Si  -99.99  -99.99 

r.« 

TllQd  t 

T  lw*p 

3.69 

:.3s  2.2i  i a. 92 

Thatab 

HtuD«  Qqq 

'  .428 

1  .8®9€»«3  2  .  S*?8€*03 

0*  1  a  Sat  Nu*bar  •  27  gait  Oil  t  • 

3.6 

TINE: 

1  7 : 86 : iS 

TC  No: 

'  2  3  4 

S  6  7  g 

Ta*p  : 

3.41  3. S3  3.S6  3.S2 

3.48  3. 4S  -99.99  -99.99 

Tua 

Tltqd  Tliqd2  Tvapr  Pjat 

T  aunp 

J.S* 

2.3S  2.2#  2.24  16.89 

-12.3 

Thatab 

Htuba  Odo 

1.266 

1 . S7IC*63  1 . 91 1 E*63 

Oata  Sal 

:  Nuabar  •  28  Sulk  Oil  t  - 

3.8 

TINE : 

17:81:84 

TC  N«: 

•  2  3  4 

S  6  7  8 

T.N>  : 

3.24  3.42  3.34  3.34 

3.38  3.32  -99.99  -99.99 

Tua 

Tliqd  Tliqd2  T»apr  P»at 

T  luae 

3.32 

2.2S  2. IS  2.14  18.88 

-12.4 

Thatab 

Htuba  Ode 

I.I6S 

1.1596*63  1.373£*63 

Oata  Sat 

Nu«bar  ■  29  Sulk  Oil  I  • 

3.8 

TINE: 

17: I8:SS 

TC  No: 

12  3  4 

S  6  7  8 

Ta«o  : 

3. IS  3.23  3.18  3. IS 

3.12  3.28  -99.99  -99.99 

Too 

Tliqd  Tliqd2  T*apr  Piat 

Tsune 

3.17 

2.28  2.2S  2.23  18.89 

-12. S 

Thatab 

Htuba  Odd 

.937 

1 . 85SE*83  9.8936*62 

Oata  Sat 

Nuabar  •  38  Sulk  Oil  l  * 

3.8 

TINE: 

17:13:34 

TC  No: 

<234 

S  S  7  g 

T,«*  : 

i.28  3.48  3. 38  3.33 

329  3.28  -99  99  -99.99 

Tub 

Tliqd  Tliqd2  T.aor  Ptat 

T  tune 

3  34 

2.21  2.88  2.18  18.77 

-12.6 

Thatab 

Htuba  Ode 

I.26S 

1.2468*83  i.S8t£*83 

Oata  Sat 

Nu«bar  *  31  Bulk  Oil  t  • 

3.8 

TINE: 

•7: 18=86 

TC  No: 

1  2  3  4 

S  S  7  8 

Taap  : 

3. IS  3.2S  3.28  3.IS 

3.13  3.23  -99.99  -99.99 

Tua 

Tliqd  Tliqd2  Tvapr  P»at 

T  »u«p 

3.H 

2.36  2.27  2.28  18.93 

-12.7 

Thatab 

Htuba  Ode 

.9*7 

9.922E*62  9.66i£*62 

Oata  Sat 

Nu«bar  •  32  lutk  Oil  I  * 

3.6 

TINE: 

17:17:39 

TC  No: 

12  3  4 

5  6  7  8 

Ta*p  : 

2.92  3.88  2.92  2.89 

2.91  3.63  -99.99  -99.99 

Tua 

TllQd  T|,ad?  T»apr  Plat 

T  aunp 

2.94 

2.15  2.17  2.15  18.81 

-12.7 

Thatab 

Htuba  Odp 

.79S 

8.888E*82  7.874E*82 

Oata  Sat 

Nu*bt'-  »  33  Suit  Oil  I  • 

3.6 

TINE: 

17:19:16 

TC  No: 

'  2  3  4 

5  6  7  8 

Ta»p  : 

2.92  2 . 95  2.88  2.86 

2.87  2.99  -33.59  -99.99 

Two 

TllQd  Tliqd2  T»apr  P»at 

T  lute 

2.91 

2.17  2.19  2. IS  18.82 

-12.8 

Thatab 

Htuba  Odp 

.7S3 

7 . 74SE  *82  S  8:96*62 
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APPENDIX  C.  SAMPLE  CALCULATIONS 


Data  file  number  D1221SM0  was  used  to  conduct  the  sample 
calculations.  The  saturation  temperature  was  2.28  °C  with  a 
single  smooth  tube  heat  flux  of  31290  W/m2 . 

A.  TEST-TUBE  DIMENSIONS 

0o= .015875  m 
Di=.0127  m 
D1=. 01245  m 
L= . 2032  m 
Lu=.0762  m 

B.  MEASURED  PARAMETERS 

V=141.95  volts 
1=2.34  amps 
Tl  =  16 . 84  °C 
T2-17.15  °C 
T3  =  16 . 42  °C 
T4=16 . 22  °C 
T5=15 . 66  °C 
T6=15 . 48  °C 
Tsat=2.28  °C 
kc=344  W/m-K 
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C.  OUTER  WALL  TEMPERATURE  OF  THE  BOILING  TUBE 


p=7r-D0=7T'  .015875  m  =.04987  m 

Ac  =  7r  (D02-D12)  /4  =  tt(  (.015875)2  -  (  .  '"'12  7 )  2  )  /4 
Ac=7 . 13x10 " 5  m2 

Qh=VI=141 . 95 • 2 . 34=  332.16  W 

Tavg=Ern/  n=E  6Tn/  6  =  1 6 . 3  °C 

Two=Tavg-!-  Qh-  (ln(D0/D1)/(2ir-L-kc)  )j 

r  ; 

Two=16 .3-t  332.16-  (In  (.  015875/ .  01245)  /  (2ir- L  •  kc)  )l 

Two=16 . 12  °C 


0=TWO-Tsat  =  16. 12-2, 28  =  13.84  °C 

D.  PROPERTIES  OF  R-124  AT  FILM  TEMPERATURE 

The  thermophysical  properties  of  R-124  were  obtained  from 
REFPROP  [Ref.  12]  and  are  shown  in  Appendix  A. 

Tf=(Two+Tsat)/2=(16.12+2.28)/2=  9.2  °C 

M  =  3.4  593-  (  .  0426) Tf+( 3. 9485x10 '4)Tf 2  -  ( 4 . 193x10  * 6 )  Tf 3  + 

(2 . 0709x10 '8)  Tf4 
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p=3.4593-  ( .0426)  • 9 . 2+ ( 3 . 9485x10 "4 )  ■  (9.2) 2- (4.193xl0‘6)  •  (9.2) 3  + 
(2 . 0709x10  "8 )  •  ( 9 . 2 ) 4 
p=3 .098xl0'4  N • s/m2 

p=1434 . 8  -  ( 2 . 8619 )  Tf  -  ( 6 . 72 67x10 ' 3 ) Tf 2  -  (7 . 2 852x10 '5) Tf 3 

p=1434 . 8  -  (2.8619)  -9.2-  (6 . 7267x10 '3)  •  ( 9 . 2 ) 2  -  ( 7 . 2852x10 ‘ 5 )  ■  (9.2) 3 

p=1407 . 84  kg/m3 

u=p/ p  =  3  . 098x10 "4/ 14 07. 84=  2.201xl0'7  m2/s 

k=7 . 5191x10 *2-  (3. 5436x10 '4)Tf-  (1 . 9545x10 ~7) Tf2+ (3 . 1835x10 ‘9)  Tf 3 
k=7 . 5191xlO'2  -  (3 . 5436xl0‘4)  •  (9.2)  -  (1 . 9545xl0'7)  ■  (9.2,2  + 

( 3 . 183  5x10" 9 )  ■  (9 .2) 3 
k=7. 192xl0'2  W/m-K 

Cp=l .  0542-t-  (2 . 1405x10  ‘3 )  Tf  +  ( 1 . 0709x10 '5 )  Tf  2  -  (6.4721xl0‘8)Tf3- 
(1. 4  324x10  "9)Tf4+  (4 . 13 6x10 " 11 )  Tf5 

Cp-1. 0542+ (2. 1405x10 '3)  ■  (9.2)  +  ( 1 . 0709x10 * 5 )  •  (9. 2)2-  (6.4721 
xlO"8)  •  ( 9 . 2 ) 3  -  (1. 4324x10 *9)  •  (9 . 2 ) 4+ (4 . 136x10’ 11 )  •  (9.2) 5 


Cp =1074. 74  J/kg-K 

a=k/ (p • Cp) =  7 . 192x10 '2/ (1407.84-1074.74) 
a=4.753xl0*8  m2/s 
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0=- (Ap/AT) / p-2  . 094xl0'3  (1/K) 


Pr=  u/a=4.63 


E.  HEAT - FLUX  CALCULATION 

The  average  natural  convection  heat • transfer  coefficient 
at  the  non-boiling  ends  of  the  test  tube  is  calculated  using 
the  Churchill  and  Chu  [Ref.  4]  correlation: 

For  simplification  let  &  represent  the  following  in  the 
correlation  : 

R=  [  (g-0-Do3  •  0  •  tanh  (m-  Lu)  )  /  ( u  ■  or  Lu  ■  m)  ] 1/6 
then 

h=k/DQ[  .6+.  3 87  [R/ (1+  (0 . 559/Pr)  9/16)  8/27]  ]2 
where 

m=  [  (h-p)  /  (kc- Ac)  ] 1/2 

and  an  iterative  process  is  used  in  program  DRPGB  to  compute 
the  value  of  h  by  starting  with  the  value  for  h  of  190 
(W/m2 • K) .  The  resulting  values  for  h  and  m  were: 

h=265 . 43  (W/m2 • K) 
m=23 .23  (1/m) 
therefore : 

Qf= (h- p- kc- Ac) 1/2 •  6 • tanh (m  -Lu)=7.44  W 
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F.  HEAT  FLUX  THROUGH  ACTIVE  BOILING  LENGTH 


Q=Qh-2 ■ Qf  =  3  32 . 16  -  (2-7. 44 ) =3 17. 2  8  W 
Ab=7r  •  D0  •  L=tt  ■  0 .015875  ■  0 .2032  =  1 . 0134x10 '2  m2 
q=Q/Ab=3 17 . 28/1 . 0134x10  "2  =  3 1 ,3  08  W/m2 
h=q/ 0  =  31,308/13.84=2262.1  W/m2  •  K 

The  following  results  were  produced  by  the  data 
acquisition  and  reduction  program  DRPGB: 
q=31290  W/m2 
0=13.84  °C 
h=2262  W/m2 • °K 


TABLE  5.  DIMENSIONS  OF  BOILING  TUBES  TESTED 


TUBE 

D1  mm 

D2  mm 

Di  mm 

Do  mm 

L  mm 

Lu  mm 

Kc 

W/moK 

smooth 

12.4 

15.9 

12.7 

15.9 

203.2 

76.2 

344 

GEWA-K 

19  fpi 

10.1 

12.7 

11 . 1 

12.7 

203.2 

76.2 

344 

GEWA-K 

26  fpi 

10.1 

12.7 

11.1 

12.7 

203.2 

76.2 

344 

TURBO  - B 

11.6 

13.8 

11 . 8 

15.8 

203.2 

76.2 

398 

HIGH 

FLUX 

12.9 

15 . 8 

13.2 

15 . 8 

203.2 

76.2 

45  * 

*  HIGH  FLUX  tube  porous  coating  is  on  a  copper/nickel  base 
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APPENDIX  D.  UNCERTAINTY  ANALYSIS 


Using  methods  of  uncertainty  analysis  outlined  by  Kline 
and  McClintlock  [Ref.  18]  and  performed  by  Sugiyama  [Ref.  2), 
four  data  points  are  analyzed.  The  data  points  chosen  were  for 
a  Smooth  tube  and  HIGH  FLUX  tube.  For  each  tube,  there  was  a 
data  point  at  a  high  heat  flux  and  one  at  a  low  heat  flux.  The 
following  is  a  sample  uncertainty  analysis  calculation  from 
data  file  number  D1220SM0  with  a  Smooth  tube,  0%  oil,  at  a 
heat  flux  of  86,000  W/m2 .  The  results  of  the  remaining  sample 
points  are  listed  in  Table  6.  Uncertainties  are  given  as  a 
percentage  of  the  calculated  parameter,  i.e.  a  relative 
uncertainty. 

A.  UNCERTAINTY  IN  HEAT -  TRANSFER  RATE 

Ig=1.99  volts  Vg=9 .325  volts 

Estimated  uncertainties  in  the  measured  quantities  I  and  V  : 
6I=+  0.025  amps 
6V=+  0.05  volts 


1=1.9182 • Is=3 . 82  amps 

V=25 • Vg=233 . 13  volts  Qh=VI 

where:  5=  uncertainty  in  measurement 
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<5Qh/Qh=  (  (6V/Vs)2+(5I/Is)2)1/2 
5Qh/Qh=  ( (.05/9. 325) 2+(. 025/1. 99) 2) 1/2 
6Qh/Qh=  1.37  % 

B.  UNCERTAINTY  IN  SURFACE  AREA 

Ab=  tt-Do-L 

Dq=  15.88  (mm)  6D0=  .1  (mm) 

L=  203.2  (mm)  <5 L=  .1  (mm) 

6Ab/Ab=  (  (  6D0/D0) 2  +  ( 5L/L) 2 ) 1/2 
6Ab/Ab=  ( ( . 1/15 . 88 ) 2  +  (. 1/203 .2) 2) 1/2 
6Ab/Ab=  .63  % 

C.  UNCERTAINTY  IN  WALL  SUPERHEAT 

Two-Tsat  Tgat=2 .23  °C  6Tgat=  .01  °C 

Two=  Tavg  -  Qh  C 1  n  ( VDj. )  /  ( 2  •  7T  •  L  •  kc )  ] 

Tn=  thermocouple  readings 

^avg=  ETn/n 

Tl-21.63  °C 

T2-21.41  °C 

T3  =  20 . 63  °C 

T4=20 . 94  °C 

T5=20 . 06  °C 

T6=2 0 . 40  °C 

Tavg=20.85  °C 
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standard  deviation=  (  (£  (Tn-Tavg)  2)  /n) 1(/2 
standard  deviation=  0.55  °C 

The  logarithmic  term  in  the  equation  for  Two  can  be 
neglected  for  the  uncertainty  analysis  as  it  is  very  small  in 
comparison  to  the  standard  deviation. 

Two=Tavg=20.8  5  °C 

6  Two= standard  deviation** .  55  °C 
AT=18 . 62  °C 

6 AT/ AT =  (  (<5Two/AT)  2+  ( 5Tgat/AT)  2) 1/2 
6AT/AT= ( (. 55/18. 62) 2+( .01/18. 62) 2) 1/2 
<5AT/AT=2 . 95  % 

D.  UNCERTAINTY  IN  HEAT  FLUX 

q=  (Qh-2-Qf)/Ab 

Qh=VI  Qh=890 . 56  W  <5Qh=12.2  W 

assuming  the  same  proportion  in  the  uncertainty  for  Qf: 

q=86000  W/m2  Ab=  0.010137  m2 

Qh- 2 • Qf =  871.81  W 

Qf =  9.38  W  5Qf=0 . 128  W 

5q/q= [ (6Qh/(Qh-2Qf) )2+(26Qf/(Qh-2Qf) )2+(6Ab/Ab)2]1/2 
5q/q=[ (12. 2/(890. 56-2-9. 38) )2+{2-(. 12 8)/{ 890. 56-2-9. 38) )2+ 

( 6 . 39x10 "5/.010137)23  1/2 
6q/q=  1.54  % 
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E.  UNCERTAINTY  IN  BOILING  HEAT- TRANSFER  COEFFICIENT 


h=q/AT 

6h/h=  [  (6q/q) 2+ (6AT/AT) 2] 1/2 
6h/h=  [(.0154)2+(.0295)2]1/2 
5h/h=  3.3  % 


TABLE  6.  UNCERTAINTY  ANALYSIS  OF  FOUR  DATA  POINTS 


Parameters 

Smooth  tube 

4019  W/m2 

Smooth  tube 

86000  W/m2 

High  Flux 

2244  W/m2 

High  Flux 

88860  W/m2 

6Qh/Qh  (%) 

3.52 

1.37 

8.11 

1.36 

6Ab/Ab  (%) 

.63 

.63 

.634 

.634 

5AT/AT  (%) 

5.20 

2.95 

17.4 

12.5 

6q/q  (%) 

1.29 

1.54 

8.41 

1.5 

6h/h  (%) 

5.4 

3.3 

19.3 

12.6 

All  uncertainty  analysis  calculations  were  conducted 
during  decreasing  heat  flux  runs.  The  wall  superheat 
uncertainty  (6AT/AT)  for  the  HIGH  FLUX  tube  was  significantly 
greater  than  that  of  the  smooth  tube  due  to  a  greater 
deviation  in  the  thermocouple  temperature  readings. 
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APPENDIX  E.  SETUP  PROGRAM 

The  following  setup  program  was  used  in  the  preparation  of 
the  apparatus  prior  to  commencing  testing.  The  program  was 
written  in  Hewlett-Packard  Basic  5.0  for  both  the  Hewlett- 
Packard  9300  and  9852A  series  data  acquisition/control  unit. 
The  setup  program  performs  the  following: 

1.  Monitor  the  Sump  temperature 

2.  Monitor  the  evaporator  liquid  temperature 

3.  Measurement  and  readout  of  all  thermocouple  channels. 

4.  Measurement  and  readout  of  the  power  supplied  to  the 

tube  cartridge  heater. 

5.  Measurement  and  readout  of  the  power  supplied  to  the 
auxiliary  heaters  if  used. 
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1  i  PROGRAM:  SETUP 

2  '  DATE:  AU6UST  3.1991 

4  i  PROGRAMMER;  LT  DEAN  SU6I YAMA 

5  1  MODIFIED  BY  LANNIE  LAKE  JAN  22.  1992 

10  COM  /Cc/  C(  7  > 

20  DATA  0. 1 0086091  .25727.94369  ,-767345.9295 .7802S595.9I  .-9247486589 .5.975=8 

1 1  .-2.661 92E»  l 3 

21  DATA  3. 94078E* I  4 

22  READ  Ct •  ) 

23  ON  KEY  1  .15  GOTO  27 

25  PRINTER  IS  I 

27  PRINT 

28  PRINT 

30  PRINT  USING  *4X. "SELECT  OPTION*" 

31  PRINT  USIN6  *6X  ."0-MONITOR  SUMP*" 

32  PRINT  USING  *6X.**l -MONITOR  LIQUID"* 

33  PRINT  USING  *6X .* ‘2-CHECK  THERMOCOUPLES* ” 

34  PRINT  USING  *6X .** 3-CHECK  MAIN  HEATER*** 

35  PRINT  USING  *6X." 4-CHECK  AUX  HEATERS"* 

36  PRINT  USING  *6X,* -5-EXIT  PR06RAH*  *  * 

37  PRINT  USIN6  *4X  .  *  “NOTE :  KEY  I  -  ESCAPE"* 

38  BEEP 

40  INPUT  Ido 

41  IF  I  do; 5  THEN  Ido-S 

42  IF  Ido-0  THEN  50 

43  IF  Ido- 1  THEN  155 

44  IF  Ido-2  THEN  173 

45  IF  Ido-3  THEN  195 

46  IF  Ido-4  THEN  195 

47  IF  Ido*5  THEN  231 

48  PRINT 
49' 

50  PRINT 

51  PRINT  ’SUMP  TEMPERATURE  DEG  C  * 

53  PRINT 

54  OUTPUT  709 1  *  AR  AF  11  ALII  VR5  * 

60  OUTPUT  709 | * AS  SA* 

70  So"-0 

80  FOR  J-l  TO  5 

90  ENTER  709 i E 

100  Su"-Sun»E 

1  10  NEXT  J 

120  Eave-Sun/S 

130  Te«p-FNTvjv( Eavi ) 

140  PRINT  USING  *4X ,MCO.OO*;T«np 

141  BEEP 

142  PRINT 

150  WAIT  5 
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'51  GOTO  50 
1 50 1 

'55  PRINT 

156  PRINT  -LIQUID  temperature  066  C 
153  PRINT 

159  OUTPUT  709i "AR  AP08  AL09  URS' 

1 60  Sum*0 

IS  I  FOR  I- 1  TO  2 

162  OUTPUT  709; "AS  SA' 

160  ENTER  709;E 

164  5un*SuwtE 

165  NEXT  I 

166  E*ve*Sun/2 

167  Tenp«FNTvs»( Eave i 

168  PRINT  USING  ’4X  ,MOC.OD*;Tenp 

169  BEEP 

170  WAIT  5 

171  60T0  155 
172' 

173  PRINT 

176  PRINT  -CHANNEL  TEMPERATURE  OE6  C* 

177  OUTPUT  709)  “AR  AF00  ALII  UR5* 

178  FOR  I«l  TO  12 

179  OUTPUT  709)  "AS  SA- 

180  Su*-0 

181  FOR  J« I  TO  5 

182  ENTER  709. E 

183  5u«“Su«»E 

184  NEXT  J 

185  E«v**Su*/5 

186  T«pip»FNTy*v(E  oil 

187  PRINT  TA8( 3  1) I iTABl IS  IjTanp 

188  NEXT  I 

189  BEEP 

190  WAIT  5 

191  60T0  173 
1941 

1 95  PRINT 

196  OUTPUT  709) ’ AR  AF20  AL22  uR5‘ 

197  FOR  I«1  TO  3 

198  OUTPUT  709. -AS  SA- 

1 99  Sum«0 

200  FOR  J-i  TO  5 

201  ENTER  709.E 

202  Suo'Sum’E 
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203  NEXT  J 

384  IF  1-1  THEN  Volt-Su«'5 
205  IF  1-2  AND  Iao*3  Then 

209  PRINT  'MAKE  SURE  V0LTA6E  BOX  IS  SET  TO  HA;N  HEATERS' 

2 1 0  Anp-Sun/S 

211  END  IF 

212  IF  1-3  AND  Ido-4  ThEN 

216  PRINT  'HAKE  SURE  vOUTASE  box  is  SET  to  *ux  HEATER S  ’ 

2 1 7  A«p-Su«/S 
2)3  END  IF 
219  NE«T  I 

223  Ahp»A8S< Anp * ) .9)92  > 

221  Voi  t *A8S( Vol t «2S  ) 

222  Poi-er-Uolt -Anp 

223  Re* i s t ence-Vo 1 t /Anp 

224  PRINT 

225  BEEP 

226  PRINT  *  VOLTAGE! V  )  CURRENT  (  A )  RESISTENCE!  oh*s >  POWER! U)* 

227  PRINT 

228  PRINT  USIN6  *  I  X  ,5<  flOOOO  .  00  ,4X  )*  I  Uol  t  ,A«p  .Resist*  nee  .Power 

229  WAIT  5 

230  GOTO  t  95 

231  BEEP 

232  PRINT 

233  PRINT  'THAT’S  ALL  FOLKS '* 

234  END 

235  OEF  FNTyJy(V) 

236  CON  /Co/  C<7> 

237  T-C<  0  ) 

238  FOR  I-l  TO  7 

239  T*T»C< I  > •  V"  I 

240  NEXT  I 

24)  t«T*8 . 626897E-2tT  *< 3.761 l 99E-3-T-5. 0689259E-5 ' 

242  RETURN  T 

250  FNENO 
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APPENDIX  F.  PROGRAM  DRPGB 


The  data  acquisition  and  reduction  program  DRPGB  is 
written  in  Hewlett-Packard  Basic  5.0  for  the  HP  9300  series 
computer  and  listed  on  the  following  pages. 


108 


10  i  FILE  NAME:  DRP6B 

20  l  DATE:  October  19.  1984 

30  I  REVISED :  July  M,  1991  BY  Oean  Sugiyana 

41  1  REVISED:  Aug  18.  1992  8V  Seorge  Bertteh  < including  R-124  oroperttesi 

50  COM  / Ido /  Ido 

80  PRINTER  IS  I 

70  CALL  Select 

80  INPUT  "WANT  TO  SELECT  ANOTHER  OPTION  (  1  -Y  ,0-N  )-> *  leal 

90  IF  Isal'l  THEN  SOTO  70 

1 00  BEEP 

M0  BEEP 

120  PRINTER  IS  1 

130  PRINT  ‘DATA  COLLECT  I ON/ REPROCESS  INS  COMPLETED" 
l 40  END 

150  SUB  Mam 

151  COM  / Iproo/  Ift 
160  COM  /I do/  Ido 

170  COM  /Ce/  C<  7  )  .leal 

180  COM  /Uil/  02, Oi  .Oo.L.Lu.Kcu 

1  90  DIM  Enf ( 12  1  .T< 12  )  .Ola< 13  >  ,02a<  i  3  )  .Dial  1  3  )  ,0oa(  13  )  ,La<  1 3  )  ,Lua<  1 3  >  ,Kcua<  1  3  > . 
Ett  19).Tn0(4)ClS] 

200  DAT*  «. I008S09I .25727.94389,-787345.8295,78025595.81 
210  OATA  -9247488589, 6. 97688E-1 I .-2 .66 l 92E» 1 3 .3 . 9407BE* 1 4 
220  READ  C(*> 

230  1  DATA  "Snooth" ."High  Flu* " . "Tharnoexel-E" , "Thernoexel-HE" 

240  DATA  Snooth, High  F lux  .Turbo-6 .High  Flux  Mod.Turdo-8  Mod 
250  READ  TnS( •  ) 

260  PRINTER  IS  701 

270  BEEP 

280  IF  Idd-4  THEN  PRINTER  IS  1 
290  IF  Ido-4  THEN  SOTO  2880 

300  I  INPUT  "ENTER  MONTH,  OATE  AnC  TIME  ( MM: DD : HH : MM :SS > "  ,0ata* 

310  ■  OUTPUT  709i "TO" (Date* 

320  i  OUTPUT  709 i "TO" 

330  1  ENTER  709iOataS 
340  PRINT 

350  1  PRINT  "  Month,  date  and  tine  :"iDate» 

351  PRINT  "  Data  :  ’  ,DATE0<  TImEOATE  ) 

380  PRINT 

370  PRINT  USIN6  "10X,""NOTE:  Progran  nax la  :  DRPSB’"" 

380  BEEP 

390  INPUT  "ENTER  DISK  NUMBER"  .On 

400  PRINT  USIN6  *  1 6X  , " ”0 1  at  nunbar  -  ••.22"1Dn 

410  BEEP 

420  INPUT  "ENTER  INPUT  MODE  < 0-3054A  .  I -FILE >"  . In 

430  BEEP 

431  Ift-0 

432  INPUT  "SELECT  FLUID  <0-1 14 .1-124)  -  .Ift 

433  BEEP 

440  INPUT  "SELECT  HEATIN6  NOOE  (0-ELEC:  I -WATER  1  *  , Ihn 
450  BEEP 

460  INPUT  "ENTER  THERMOCOUPLE  T7PE  ( 0-NEW  .  I -OLD >"  .leal 
470  IF  In-0  THEN 

480  SEEP 

490  INPUT  "SIVE  A  NAME  FOR  THE  RAW  DATA  F IlE"  ,02_f i let 
500  PRINT  USIN6  " 1 sx  ,""Nau  file  nana:  " * . 1 4A" i02_fi la« 

510  Size) -20 

520  CREATE  BOAT  02_f i let  .Size! 

530  ASSI6N  0FUe2  TO  02_filet 

540  > 

550  ’  OUMMY  FILE  UNTIL  Nrun  KNOWN 
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sea  oi_fii»«-‘Ounnr- 
S7a  CREATE  BOAT  OI_f  ila«  .Snel 

580  ASSI6N  0Filel  TO  Dt.filef 

590  OUTPUT  fFilel iOate< 

S00  IF  Ihn-0  THEN 
810  SEEP 

620  INPUT  'ENTER  NUMBER  OF  DEFECTIVE  TCS  (0-OEFAULT.  «  MAX.l'.Idtc 
630  IF  Idtc-0  THEN 
640  Ldtel-0 

650  Ldtc2-0 

651  Ldte3*0 

652  Ldtc4-0 

660  PRINT  USIN6  '  1 6X  ,  "No  defective  TCa  aviat'” 

670  ENO  IF 

680  IF  Idte«t  THEN 

690  BEEP 

700  INPUT  'ENTER  OEFECTIUE  TC  LOCATION  <1-8>'.Ldtc! 

710  PRINT  USIN6  '  1 SX  ,  "TC  la  defective  at  location  ".OO'iLdtcl 

720  Ldtc2*0 

721  Ldtc3-0 

722  Ldtc4-0 

730  ENO  IF 

740  IF  Idte«2  THEN 
750  BEEP 

760  INPUT  'ENTER  OEFECTIUE  TC  LOCATIONS  (  I  -a  >- .Ldtcl  ,Ldtc2 

770  PRINT  JSIN6  'I6X,"TC  ara  defective  at  location*  "  ,DD  ,4X  ,00'  iLdtc  I  ,Ldtc2 

771  Ldtc3“0 

772  Ldte4-0 

790  ENO  IF 

781  IF  Idte-3  THEN 

782  BEEP 

783  INPUT  'ENTER  OEFECTIUE  TC  LOCATION  (  I -8 )*  .Ldtcl  ,Ldtc2  ,Ldtc3 

784  PRINT  USING  '16X.*'TC  i*  defective  at  location  "  ,3<4X .00 >' iLdte I ,Ldte2  .Ld 

tc3 

785  Ldtc4-0 

786  ENO  IF 

787  IF  Idtc-4  THEN 

788  BEEP 

799  INPUT  'ENTER  OEFECTIUE  TC  LOCATION  (  I -8 ) ' .Ldt c I  ,Ldtc2 ,Ldtc3 ,Ldt c4 

790  PRINT  USIN6  ’16X,"TC  la  dafactive  at  location  "  ,4(  4X  ,DD  > '  iLdt  c  I  ,Ldt  c2  ,Ld 

t  c3  ,Ldtc4 

792  ENO  IF 

793  IF  Idtc>4  THEN 

800  BEEP 

810  PRINTER  IS  I 

820  BEEP 

830  PRINT  ' INUALIO  ENTRY' 

840  PRINTER  IS  701 

850  GOTO  610 

860  ENO  IF 

870  ENO  IF 

880  OUTPUT  0F  i  lei  iLdte  I  ,Ldtc2,Ldtc3.Ldtc4 

890  1  In* I  oot Ion 

900  ELSE 

9<0  BEEP 

920  INPUT  '6IUE  THE  NAME  OF  THE  EXISTING  DATA  PILE’  ,02_f i lal 

930  PRINT  USING  'ISX."01d  file  nan*:  "  , )  4A' ,  02_f  i  la* 

940  ASSI6N  0Fila2  TO  02_file» 

950  ENTER  0F l le2 iNrun 

960  ENTER  0Fila2iOold* 

970  PRINT  USIN6  '  1 6X  , "Thi*  data  aat  taxan  on  :  ' '  . 1 4A' lOolrt* 

980  ENTER  0FUa2iLdtcl  ,Ldte2  .Ldtc3  ,udtc4 

990  IF  Ldtcl >0  OR  Ldtc2>0  OR  Ldtc3>0  OR  Ldtc4>0  THEN 

1000  PRINT  USING  '  I6X  ."TMemocouple*  were  defect  iva  at  locat  tone  :  "  ,4<  30  ,4X  )’ i 
Ldtcl  ,Ldtc2  ,Ldtc3  ,Ldte4 
'010  ENO  IF 
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1020  ENTER  0F i i»i 1 1 1 1 
1030  END  IF 
1040  Idtc-0 

1050  IF  Ldtc1>0  THEN  Idte-Idtc+1 
10S0  IF  Ldtc2>0  THEN  Idte-tdtc*l 
(061  IF  Ldtc3>0  THEN  Idtc-Idte+1 
1 0S2  IF  Ldtc4>0  THEN  Idte-Idtc- t 
10701  IF  In-0  ANO  Ihn-I  THEN  1595 
1080  BEEP 

1090  input  'want  to  create  a  plot  file’  ib-n.i-v  >*  ,iBiot 

1100  IF  Iplot-1  THEN 
1110  8EEP 

1120  INPUT  •GIVE  NAME  FOR  PLOT  FILE*  .P_*  i  le* 

1130  CREATE  80AT  P_file«  .4 
1140  ASSI6N  #Plo»  TO  P_file* 

1 150  END  IF 
1180  IF  Ihn-I  THEN 
1 1 70  BEEP 

1180  INPUT  •UANT  TO  CREATE  Uo  FILE’  < 0-N . I -Y  )*  , Iuf 
1190  IF  Iuf-I  THEN 
1 200  BEEP 

1210  INPUT  “ENTER  Uo  FILE  NAHE‘  ,Uf lie* 

1220  CREATE  BOAT  Ufile*,4 
1230  ASSIGN  tUfile  TO  Ufile* 

1240  END  IF 
1 250  BEEP 

1250  INPUT  -want  TO  CREATE  Re  FILE’  (  0-N . 1 -V >•  . Ire 
1270  IF  Ire- I  THEN 
1280  BEEP 

1290  INPUT  ‘ENTER  Re  FILE  NAME' , Ref  lie* 

1300  CREATE  BOAT  Ref  lie*.  10 
1310  ASSIGN  tRef i le  TO  RefUe* 

1320  END  IF 
1330  END  IF 
1340  PRINTER  IS  1 
1350  IF  In-0  THEN 
1 3S0  BEEP 

1370  PRINT  USIN6  *4X  ,  *  "Select  tube  nunBef 
1380  IF  Ihn-0  THEN 

1390  PRINT  USIN6  -6X.--0  Snooth  4  inch  Ref- 

1400  PRINT  USIN6  -SX.**1  Snooth  4  inch  Cu  ( Pre*i/Sl ide  )* ’ " 

1410  PRINT  USIN6  "6X  . ”2  Soft  Solder  4  inch  Cu'"* 

1420  PRINT  USIN6  "6X.--3  Soft  Solder  4  inch  HI&H  FLUX*** 

1430  PRINT  USING  *SX.**4  Uielend  Herd  8  inch"* 

1440  PRINT  USING  *GX . *  *5  HIGH  FLUX  8  inch*** 

1450  PRINT  USING  *6X . *  *5  GEWA-K  40  Fine/in"’ 

1460  PRINT  USIN6  *5X , * *7  GEUA-K  26  Fin»/in*** 

1470  PRINT  USIN6  *SX  ,*  *8  GEWA-T  19  Fin»/in-** 

1480  PRINT  USIN6  *GX  . ■ *9  GEWA-T  OR  6EWA-TY  2G  Fins/in**- 

1490  PRINT  USING  *BX ,*  *  10  THERNOEXCEL-E" * ’ 

1500  PRINT  USIN6  *GX .**11  THERMOEXCEL-HE**- 
1510  PRINT  USIN6  "GX , *  *  1 2  TURBO-8--- 
1520  PRINT  USING  *6X , *  *  1 3  GEUA-X  19  Fin»/in-'- 
1530  ELSE 

1 540  PRINT  USING  *GX.-*0  Snooth  tuOe"" 

1550  PRINT  USING  *SX,-*1  High  Fluf 

l 5S8  PRINT  USING  *6X ,**2  Turoo-B"* 

1570  PRINT  USIN6  *GX.**3  High  Flu«  Hod" " * 

1580  PRINT  USING  -GX,**4  Turbo-8  Hod"- 

1590  ENO  IF 

>608  INPUT  Itt 

1610  OUTPUT  GFilel  iltt 

1620  ENO  IF 

1630  PRINTER  IS  701 

1640  IF  Itt<l0  THEN  PRINT  USING  MSX.’-Tube  Nunber:  **,0*iltt 

1 6S0  IF  1 1 1  >  9  THEN  PRINT  USING  •’GX.-'T-ioe  Nunoer:  '  .DC  lit*. 
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1660  IF  inn-1  T Hth  PRINT  uSihb  IB*,  uot  lyoe: 

1670  SEEP 

)S8S  INPUT  'ENTER  OUTPUT  VERSION  < 0-LON6 . 1 -SHORT  ,2-NONE >' ,Iov 
1690  SEEP 

1700  INPUT  'SELECT  1 0-LIQ  ,1 -UAP  ,2-< LIQ+UAP 1/2  1' . Ilqv 

17101 

1720'  DIMENSIONS  FO  TESTED  TUBES 

1730'  ELECTRIC  HEATED  MODE 

17401  D1-0ianeter  at  tharnoeouple  positions 

1750  OATA  .01 1 1 12S..01 1  1 125, .01 1 1125.. 01 29540,. 01 2446.. 01 29540.. 0100SSS 
1760  OATA  .0100965, .01 157  ..01 1S7  ,.01 1S7..01 1S7..01 157  ..0100965 
1770  READ  Dial*) 

1780  O)-Cla(Itt) 

1 790' 

I800i  02-Oianeter  of  test  section  to  the  base  of  fins 

1810  OATA  .015875 , .015875  ,  .015875 ,.015824,  .015875  .  .015824  ..01270 

1820  OATA  .0127, .0138. .0138, .0138, .0138, .0138, .0127 

1830  READ  D2a<»> 

18401 

18S0I  Di-Inside  dianeter  of  unenhanced  ends 

1860  OATA  .0127. .0127. .0127, .0132. .0127, .0132, .0111125, .0111125 
1870  DATA  .01 18,. 01 18  ..01 18  ,  .01 18 ..01 18 , .01 1 1 125 
1880  READ  Dial •  ) 

1890' 

1900'  Oo-Outside  dianeter  of  unanhanced  ends 

1910  DATA  . 01 5875.. 0! 5875,. 0IS87S,. 015824,. 01 5875.. 01 5824,. 01 270 ..01 270 
1920  DATA  .01331  ..01331  ,.01331  ,.0133) ,.0158, .0127 
1930  REAO  0oa<*> 

1940' 

1950'  L-Length  of  enhanced  surface 

i960  DATA  .1016..  1016..  1016,.  1016.. 2032,. 2032,. 2032,. 2032,. 2032,. 2032,. 2032,. 20 

32..  2032..  2032 
1970  READ  La< ♦ 1 
i960' 

19901  Lu-Length  of  unenhanced  surface  st  the  ends 

2000  OATA  .0254  .  .0254  ,  .0254  .  .0254 , .0782 , .0762  .  .0762  ,  .0782 . .0762 ..0762 . .0762 , .07 

62..  0762..  0762 
2010  REAO  Luel •  ) 

2020' 

2030'  Kcu-Thernel  Conductivity  of  tube 

2040  OATA  398 .344 ,344 ,45 .344 .45 .344 .344 ,398 ,398  .398  ,398  .398  .344 
2050  REAO  Kcua(*> 

2060  IF  Ihn-I  THEN 
20701 

2080'  Data  statenants  for  water  heating  node 
2090' 

2100  OATA  0.015875,0.015875,0.0169.0.0138.0.0169,0.0,0.0.0 
2110  REAO  02a< •  1 

2120  OATA  0.0127.0.0127,0.0145,0.0127,0.0145.0,0,0.0.0 
2130  REAO  Dial  <) 

2140  OATA  0.015875,0.015875,0.0169.0.015875.0.0)69.0,0,0.0.0 
2150  REAO  0oe< •  1 

2160  OATA  0.3048,0.3048,0.3048.0.3048.0.3048.0.0.0.0.0 
2170  REAO  La< •  1 

2180  OATA  0.0254,0.0254.0.0254,0.0254.9.0254.0.0.0.0,0 
2190  REAO  Lua( •  1 

2200  OATA  398.45,398,45,398,0.9.0,0.0 
2210  REAO  Kcual <  ) 

2220  END  IF 
2230  D2-02a( Itt  ) 

2240  Di-Oia(ltt) 

2250  Oo-Ooa(Itt) 

2269  L-La(Itt) 

2270  Lu-Lua(Itt) 

2289  Kcu-Kcua( Itt ) 

2290  Xn-.8 
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230(4  Ff-«.3 

2310  IF  Itt-0  THEN  Cf-1 .70E+9 

2320  IF  1 1 1 >0  THEN  CI-3.7037E+10 

2330  A-PI»(0o*2-0l'2 >/4 

2340  P-PI»Do 

2350  IF  Ihn-I  THEN 

2360  BEEP 

2370  INPUT  ‘TUBE  INITIATION  NODE.  (1-HOT  WATER , 2-STEAM  ,3-C0L0  WATER  >’  , 1 1  in 

2390  IF  Itin-1  THEN  PRINT  USIN6  '16X."TuBe  Initiate:  Hot  Water'" 

2390  IF  Itm-2  THEN  PRINT  USIN6  '16X,”Tube  Initiate:  Steen'" 

2400  IF  Itin«3  THEN  PRINT  USIN6  '16X."TuBe  Initiate:  Cold  Water'" 

2410  INPUT  'TEMP/VEL  HOO£:  ( 0-T-CONST ,U-DEC 1 1 -T-DEC  .U-CONST i  2-T-INC  .U-CONST  )' 
I  tv 

2420  IF  I  tv-0  THEN  PRINT  USIN6  '  1 6X  .*  'Teno/Vel  Mode:  Conatant/OeereasinQ'  *  * 

2430  IF  Itv-1  THEN  PRINT  USIN6  '  1 6X  ."Teno/Vel  Mode:  OeereasmQ/Constant  ■ ' ' 

2440  IF  Itv-2  THEN  PRINT  USIN6  “  1SX  .“Tena/Vel  Mode:  Inereaa mQ/Constant  * ' ' 

2450  INPUT  'WANT  TO  RUN  WILSON  PLOT’  ( 1 -Y  ,0-N  ) “  .lull 

2460  IF  Ihn-1  ANO  Iuil-0  THEN 

2470  IF  Itt-0  THEN  Cl-. 032 

2480  IF  Itt-t  OR  Itt-3  THEN  Ci-.0S9 

2490  IP  Itt-2  OR  Itt-4  THEN  Ci-. 062 

2500  BEEP 

2510  INPUT  'ENTER  CI  ( OEF :  WH- .032 ,HF- .059 ,TB- . 062  )' ,Ci 
2520  PRINT  USIN6  '  1  ex  . "Sieder-Tata  "' 

2530  PRINT  USINS  '16X  Constant  -",Z.4D'iCi 
2540  ENO  IF 

2550  END  IF 

2560  IF  Ihn-I  ANO  In-1  ANO  luil-t  THEN 

2570  IF  Itt-0  THEN  Cl-. 032 
2580  IF  Itt-I  OR  Itt-3  THEN  Cl-. 059 
2590  IF  Itt-2  OR  Itt-4  THEN  Cl-. 062 
2600  ASSI6N  *FUe2  TO  • 

2610  CALL  WilsonlCf  ,Cl  ) 

2620  ASSI6N  0FiIa2  TO  02_fila» 

2630  ENTER  0Fila2iNrun.Oold0.Ldtel.Ldte2.Ldtc3.Ldtc4.Itt 
2640  ENO  IF 
2650  Nsub-0 

2660  IF  Ido-4  THEN  Ihn-1 
2670  IF  Ihn-1  THEN  NsuB-8 
2680  Ntc-6 

2690  IF  Ihn-0  THEN  Ntc-12 

2700  J-1 

2710  S«*0 

2720  Sy-0 

2730  Sx»-0 

2740  Sxy-0 

2 ’50  Repeat:  i 

2760  IF  In-0  THEN 

2780  Ido-2 

2790  ON  KEY  1 ,15  RECOVER  2750 
2800  PRINTER  IS  I 

2810  PRINT  USINS  '4X. "SELECT  OPTION"' 

2820  PRINT  USINS  'SX. "0-TAKE  DATA'" 

2830  IF  Ihn-0  THEN  PRINT  USIN6  *6X. "1-SET  HEAT  FLUX'" 

2840  IF  Ihn-1  THEN  PRINT  USIN6  'SX,"l-SET  WATER  FLOW  RATE'" 

2850  PRINT  USINS  '6X. "2-SET  T»et'" 

2860  PRINT  USIN6  *4X."N0TE:  KEY  1  -  ESCAPE'" 

2870  8EEP 

2880  INPUT  Ido 

2890  IF  Ido>2  THEN  Ido-2 

2900  IF  IdO-0  THEN  4440 

29t0< 

2920'  LOOP  TO  SET  HEAT  FLUX  OR  FLOWMETER  SETTIN6 
2930  IF  I  do- 1  THEN 

2940  IF  Ihn-0  THEN 
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2950  OUTPUT  709|'AR  AF20  AL21  UR5‘ 

2960  BEEP 

2970  INPUT  'ENTER  OESIRED  Ode'  .Oaop 

2980  PRINT  USIN6  *4X  ,*  'DESIRED  Qdo  ACTUAL  Qdp'" 

2990  Err- 1000 

3000  FOR  I- I  TO  2 

3010  OUTPUT  709i 'AS  SA' 

3020  Sun-0 
3030  FOR  J i-l  TO  5 
3040  ENTER  709 i E 
3050  Sun-Su«fE 
3060  NEXT  Ji 

3070  IF  I-l  THEN  Uolt-Sun/S 
3080  IF  1-2  THEN  Anp-Sun/S 

3090  NEXT  I 

3091  PRINT  'U* : '  ,Uol t 

3092  PRINT  ' Is : ' ,Awp 
3100  Anp-ABS< Anp*1 .9182 > 

3110  Uolt-ABS<Uolt«2S 1 

3120  Aqdp-Uolt »A«p/< PI »02»L ) 

3130  IF  ABS( Aqdp~0qdp  >>Err  THEN 

3140  IF  Aqdp>0qdp  THEN 

31 S0  BEEP  4000  ..2 

3160  BEEP  4000,. 2 

3170  BEEP  4000,. 2 

3180  ELSE 

3190  BEEP  250.. 2 

3200  BEEP  250.. 2 

3210  BEEP  250.  .2 

3220  END  IF 

3230  PRINT  USIN6  "4X  ,HZ .30E  ,2X  ,MZ . 30£* (Oqdp . Aqdp 

3240  WAIT  2 

3250  SOTO  3000 

3260  ELSE 

3270  BEEP 

3280  PRINT  USIN6  '4X  ,MZ .30E  ,2X  ,HZ . 30E* lOodo ,Aqdp 

3290  Err-500 

3300  WAIT  2 

3310  SOTO  3000 

3320  END  IF 

3330  ELSE 

3340  BEEP 

3350  INPUT  'ENTER  FLOWMETER  SETTIN6'  ,F«o 

3360  60T0  2810 

3370  ENO  IF 

3380  ENO  IF 

3390< 

3400i  LOOP  TO  SET  Ts«t 
3410  IF  Ida-2  THEN 
3420  IF  Ikdt*1  THEN  3470 
3430  BEEP 

3440  INPUT  'ENTER  DESIRED  Tut', Olid 


3450' 

PRINT  USING  '4X ,"  OToot  AT  sot 

Roto 

Tv 

Roto'" 

3460 

Ikdt-t 

3470 

OIdl-0 

3460 

Old2-0 

3490 

Nn-1 

3S00 

Nri-Nn  MOO  15 

3510 

Nn-Nn»l 

3520 

IF  Nrs-1  THEN 

3530 

IF  Ihn-0  THEN  PRINT  USIN6  ’4X." 

Tsot 

Tldl 

Tld2  Tv 

Tsu«p"* 

3540 

IF  IhB-1  THEN  PRINT  USIN6  ’4X," 

Tsot 

Tldl 

Tld2  Tv 

Tsunp  Tjnlo 

t  Tpilo  Tout"' 

3550  ENO  IF 

3560  IF  IHm-0  THEN  OUTPUT  709i'AR  AF08  AL I  1  URS' 
3570  IF  Ihn-1  THEN  OUTPUT  709; 'AR  AF0  ALS  UR5" 
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3593  FOR  1-1  TO  5 

3590  IF  Ihn-0  AND  1>4  THEN  3S50 

3500  Sun-0 

3610  OUTPUT  709 i "AS  SA- 
3620  FOR  Ji-1  TO  20 
3630  ENTER  709tEliq 
36*0  Sun-Sun-El iq 
3650  NEXT  Ji 
3660  Eliq-Sun/20 
3670  Tld-FNT v»v( E 1  iq  > 

3680  IF  I-l  THEN  Tldl-Tld 

3690  IF  1-2  THEN  Tld2-Tld 

3700  IF  1-3  THEN  Tv-Tld 

3710  IF  I-*  THEN  Tsunp-Tld 

3720  IF  1-5  THEN  Tinlet-Tld 

3730  IF  1-6  THEN  Tout-Tld 

3740  NEXT  I 

3750  IF  Ihn-t  THEN 

3760  OUTPUT  ?09l *AR  AF00  AL00  UR5‘ 

3770  OUTPUT  709 i "AS  SA’ 

3780  Sun-0 

3790  FOR  Kk-1  TO  20 

3800  ENTER  709 i E 

3810  Sun-Sun+E 

3820  NEXT  Kk 

3830  Enff 7 >-ABS< Sun/20 > 

3840  Tp  i  1 e-En  f ( 7  )  / 3 . 96E-4 
3850  END  IF 

3860  Atld«<TId1+Tld2>».5 

3870  IF  ABS<  At  1 d-Ot 1 d  ) >  . 2  THEN 

3880  IF  At ld>Ot Id  THEN 

3890  BEEP  4000  ..2 

3900  BEEP  4000.. 2 

3910  BEEP  4000.. 2 

3920  ELSE 

3930  BEEP  250.. 2 

39*0  BEEP  250.. 2 

3950  BEEP  250,. 2 

3960  ENO  IF 

3970  Err  I -At  Id-Old  I 

3980  Oldl-Atld 

3990  Err2-Tv-01d2 

4000  01d2-Tv 

40 1 0  IF  T 1 d I > 1 00 .  THEN  4060 

4020  IF  Ihn-0  THEN  PRINT  USING  '4X ,Sl HOOD . OD ,2X ) *  I Dt Id .Tldl  ,T ld2 .Tv .Tsuno 
4030  IF  Ihn-1  AND  Idp-0  THEN  PRINT  USIN6  "*X ,7< MOD .OD ,2X 1 ’ tOt Id .Tldl .Tld2 .Tv .Te 
unp  ,T inlet  ,Tpi le 

4040  IF  Ihn-i  and  Idp»*  THEN  PRINT  USIN6  ‘ 4X ,5( HDD .00  ,2X  >  ,3< M30 .DO ,2X )*  tOt Id  ,Tl 

dl  .Tld2  ,Tv  ,T»uno  .Tmlet.Tpile.Tout 

4050  WAIT  2 

4060  SOTO  3500 

4070  ELSE 

4080  IF  A8S( At ld-Ot Id )> . 1  THEN 

4090  IF  Atld>Dtld  THEN 

4100  BEEP  3000,. 2 

4110  BEEP  3000,. 2 

4120  ELSE 

4130  BEEP  800,. 2 

4140  BEEP  800,-2 

4150  ENO  IF 

4160  Errl-Atld-Oldl 

*170  Oldl-Atld 

4180  Err2-Tv-01d2 

4190  01d2-Tv 

*200  IF  Ihn-0  THEN  PRINT  USIN6  *4X ,S< MOOD .DQ ,2X  >' iDt Id  .Tldl  ,Tld2 .Tv  .Tiunp 
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*210  IF  Ihn-I  THEN  PRINT  US  I N6  **X .S< HOO .DD  ,2X  i  ,3< M30 . DO  . ' *  : ' ;Ot Id .-ldt  .T id2  ,T v 
,T»unp  ,T  imat .  tpua  ,iout 
*220  WAIT  2 
*230  60T0  3500 
*240  ELSE 
4250  BEEP 
*260  Erri-Atld-Oldl 
*270  Oldl-Atld 
*280  Err2-Tv-01d2 
*290  01d2-Tv 

*300  IF  Ihn-0  THEN  PRINT  USIN6  **X  ,S<  MDDD  .DO  ,2X  )•  tOt  Id  .Tldl  ,Tld2  ,Tv  ,T  sunp 

*310  IF  Ihn-I  THEN  PRINT  US IN6  **X  ,B(  HDD.  DD  ,2X  ) •  iDt  Id  .TIdl  ,Tld2  ,Tv  .Tsunp  .T i Met 

.Toils  .Tout 

*320  WAIT  2 

4330  SOTO  3500 

*340  END  IF 

*350  END  IF 

*360  ENO  IF 

*370'  ERROR  TRAP  FOR  Ido  OUT  OF  BOUNDS 

*380  IF  Ido>2  THEN 

4390  BEEP 

*400  SOTO  2810 

4*10  ENO  IF 

4*201 

4*30'  tAKE  DATA  IF  In-0  LOOP 
4*40  IF  Ikol-1  THEN  4*80 
4450  BEEP 

4460  INPUT  -ENTER  BULK  OIL  X*  .Bop 
*470  Ikol-I 

4480  IF  Ihn-0  THEN  OUTPUT  709t‘AR  AF00  ALII  URS’ 

*490  IF  Ihn«l  THEN  OUTPUT  709t'AR  AF0  AL5  UR5' 

4S00  IF  Ihn-0  THEN  Ntc-12 

4510  FOR  I-l  TO  Nte 
4520  OUTPUT  709 i *AS  SA‘ 

4530  Sun-0 

4540  FOR  Ji-I  TO  20' 

4550  ENTER  709 i E 
4560  Sun-Sun+E 

4570  IF  I-< 9-NsuB  )  OR  I-<10-Nsud)  THEN  Et ( J 1 - 1  )-E 
4580  NEXT  Ji 
*590  Kdl-0 

4600  IF  I-<  9-Nsub  )  OR  I-<  10-Nsub  )  THEN 
4610  Eava-Sun/20 
4620  Sun-0. 

4630  FOR  Jk-0  TO  19 

4640  IF  ABS<  Et<  Jk  )-Ea»a  X5.0E-6  THEN 

4650  Sun-Sun+Et< Jk > 

4660  ELSE 
4670  Kdl-Kdl+t 
4680  ENO  IF 
4690  NEXT  Jk 

*700  IF  I-(  9-Naub  )  OR  I-<I0-Naud>  THEN  PRINT  USIN6  ■*X,"KdI  -  *  *  .DD"  |Kdl 
4710  IF  Kdl>l0  THEN 
*720  BEEP 
*730  BEEP 

*7*0  PRINT  USIN6  **X,'*Too  nuch  icattarino  in  data  -  rapaat  data  aat"*" 

4750  SOTO  2800 
*760  ENO  IF 
4770  ENO  IF 

4780  Enf ( I >-Sun/<  20-Kdl  1 

4790  next  I 

*800  IF  Ihn-I  THEN 

*810  OUTPUT  709 i -AR  AF00  AL00  URS’ 

*820  OUTPUT  709 i "AS  SA* 

4830  Sun-0 

4840  FOR  Kk-I  TO  20 
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4352 

4362 

4878 

4888 

4898 

4988 

4318 

4928 

4938 

4948 

49S8 

4988 

4978 

4988 

4998 

5888 

5818 

5828 

5830 

5848 

5858 

5868 

5070' 

5888' 

5898 

5108 

5110 

5128 

5138 

5148 

5150 

5168 

5170 

5188 

5190 

5288 

5210 

5220 

5230 

5240 

5258 

5251  ' 

5260 

5270 

5288 

5298 

5380 

5310 

5328 

5338 

5348 

5358 

5360 

5370 

5388 

5390 

5408 

5418 

5420 

5438 

5448 

5458 

5460 

S470 

5488 

5498 


’0315 

E 

NEXT  Kk 

E«f<7  >-ABS(Sun>/20 
END  IF 

IF  Ihn-8  THEN 

OUTPUT  709 i "AR  AF28  AL21  UR5' 

FOR  1-1  TO  2 
OUTPUT  709i  "45  SA‘ 

Su«-0 

FOR  Ji-I  TO  5 
ENTER  709 i E 
Sum-Su—E 
NEXT  Ji 

IF  1-1  THEN  Ur-Sun/S 
IF  I«2  THEN  Ir- Sum/5 
NEXT  I 
END  IF 
ELSE 

IF  Ihm-0  THEN  ENTER  0Fila2  ;Boo  .Told*  ,Enf  (  •  )  .Ur  ,Ir 
IF  Ihn-1  THEN  ENTER  8F i Ie2 1 Bop , Told* ,Enf( •)  ,Fm» 

END  IF 

CONUERT  wf'S  TO  TEMP  ,UOLT  .CURRENT 

Tya-0 

FOR  >1  TO  Ntc 
IF  Idtc>0  THEN 

IF  I-Ldtcl  OR  I-Ldtc2  OR  I-Ldtc3  OR  1-Ldte4  THEN 

T ( I  1—99.99 

60T0  5248 

END  IF 

END  IF 

IF  Itt<4  AND  Ihn-0  THEN 

IF  1)4  ANO  I <9  THEN 

T<  I  1—99.99 

60T0  5240 

END  IF 

ENO  IF 

T (  I  >-FNTv»v<Enf<  I  1  1 
NEXT  I 

IF  Itt<4  THEN 

UARNIN6:  4  BAD  THERMOCOUPLES  ARE  NOT  POSSIBLE  FOR  ITT<4 
FOR  1-1  TO  4 

IF  I-Ldtcl  OR  I-Ldtc2  THEN 

Tua-Tw« 

ELSE 

Tya-Tya+Tl I  ) 

ENO  IF 
NEXT  I 

Twa-Tua/1 4-Idtc  1 
ELSE 

IF  Ihn-I  THEN  5450 
FOR  I- I  TO  8 

IF  I-Ldtcl  OR  I-Ldte2  OR  I-Ldte3  OR  I-Ldtc4  THEN 

Twa-Tya 

ELSE 

Tya-Twa»T( I  ) 

END  IF 
NEXT  I 

Ty-Tya/( 8-Idtc 1 
END  IF 

Tld*T<  9-Naub ) 

Tld2-T( I0-N«u6  1 
Tlda-<  Tld+Tld2  >».S 
Tv-T ( 1 l-Naub) 

IF  It t<3  ANO  Ihn-0  THEN 
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5500  Tld2— 99.99 
5510  Ty<  T  (  |0  )+T(  in)/: 

5520  END  IF 

5530  Tsu«P-T< 12-NauB > 

5540  IF  IH«-0  THEN  5570 
5550  Tinl«t-TI !3-N»ub> 

5560  Taut«T< 14-Naub  ) 

5570  IF  lh»-0  THEN 
5580  A*p-A8S< Ir«I .9182 > 

5590  Ueit»A8S< Ur >*2S 
5600  Q*Vo! t •Amp 
5610  END  IF 

5620  IF  Itt"0  AND  lhn«0  THEN 
5630  Kcu“FNKeu( Tu ) 

5640  ELSE 

5650  Kcu*Kcua( Itt ) 

5660  END  IF 

5680'  FOURIER  CONDUCTION  EQUATION  UITH  CONTACT  RESISTANCE  NE6LECTED 
5690  IF  Ihn-0  THEN  Tw«Tw-0»LQ61 02/DI  )/< 2»PI*Keu»L  ) 

5700  IF  IIqv«0  THEN  Tsat*Tlda 

5710  IF  IIqv*1  THEN  T*at«< Tlda*Tv  )• .5 

5720  IF  Ilqv-2  THEN  T*at-Tv 

5730  IF  Ihn-1  THEN 

S740  Tavg«T  mitt 

5750  6rad»37.98S3+. 104388»Tavg 

5760  Tdroo“ABS( E«f ( 7  1  )• 1 .E«-6/< I0»6rad> 

5770  T»vQC«Tinl'at-Tdrop*  .5 

5780  IF  A8S( Tavg-Tavgc  )> .01  THEN 

5790  T«vq«<  Tavg+Tavgc  )•  .  5 

5800  60T0  5750 

5810  END  IF 

59201 

58301  COMPUTE  UATER  PROPERTIES 
5840  IF  Ih««l  THEN 
5850  Ku«FNKw<  Tavg ) 

5860  Huw«-FNHuu<  Tavg  > 

5870  Cou-FNCoul TavQ ) 

5880  Pru»FNPrw( Tavg ) 

5890  Rhow«FNRhoM<Tavo> 

5900  Twj«Tavg 
5910' 


5930  Hdot-3*96S7E-3+Fna*(3.6l955E-3-F«*«(8.82006E-6-Fn»»(  1  .23688E-7-Fn»»4 .3 1  897 


E-10  >  >  > 

5940'  Hdot-Hdot' 


( 1 


.0365-T inlet  *< 1 . 96644E-3-T inlat »5.252E-6  )  )/l .0037 


5950  Kdt-0 

5960  Q-Hdot*CpM»TdroB  , 

5970  Lntd»Tdrop/L06< ( Tinlat-Taat  >/< T tnlat-Tdrop-T»at )  > 
5980  Uo*Q/(PI,Oo«L«LHtd ) 

5990  Rw«0o»L06(Do/0i >/4  2.  •Kcu  > 

6000  Tw-Taat+Fr«L«td 

6010  Uw«Hdot/<Rho«*PI,0l‘2/4 ) 

6020  Ra«*Rha«»Uw»Di/Huua 

6030  Hi«C t»K«/Ol»Ra«* .8»Prw"( 1/3.  )•<  Huwa/FNMun< T«i >  >* . 14 

6040  T«ne»Tavg-Q/<PI»Oo«L*Hi  ) 

6050  IF  A8S<  Tmi-Tuic)>.0'  THEN 

6060  T«»“( T»i+T»»ic  )».5 

6070  SOTO  6030 

6080  END  IF 

6090  Tui-(Twl+T«ic )».S 

6100  He»l /< 1 /Ue-Oe/(Oi»Hi )-Rw  ) 

61 10  END  IF 

6120  ENO  IF 

6130  IF  Ihn-1  THEN 

6140  Tfi«taO«Q/(Ho*PI»Do*L  > 
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6158  T«“T  satvThetao 
5188  ELSE 

6178  Thatab-T«-T»at 
6188  ENO  IE 
6198  IF  Thatab<0  THEN 
6288  BEEP 

6218  INPUT  “TyALL<TSAT  (0-CONTINUE.  1 -ENO)*. lav 

6228  IP  lav-8  THEN  SOTO  2760 

6238  IF  lav- I  THEN  7220 

6240  ENO  IF 

6258> 

6260'  COMPUTE  VARIOUS  PROPERTIES 
6270  Tf ii«-(Tu*Taat  )«.S 
6280  Rho-FNRho(Tfiln) 

6290  Mu-FNMui Tf i In ) 

6308  K-FNK(Tfiln) 

6310  Cp-FNCp(Tf iln > 

6320  Beta-FN8eta(Tf ll 
6330  Hf g«FNHf g( T»at > 

6340  Ni-Mu/Rho 
6350  Alpha«K/< Rho«Cp > 

6360  Pr-Ni /Alpha 
6370  Psat«FNP*at< Tjpt > 

63801 

6381  PRINT  "Rho: “ ,Rho 

6382  PRINT  “Mui'.Mu 

6383  PRINT  “K:“,K 

6384  PRINT  “Cp:“.Cp 

6385  PRINT  ■Tfilni;-  .Tfiln 

6386  PRINT  “Hfg: "  ,Hfg 

6387  PRINT  "Psat : "  ,Psat 

6388  PRINT  "Tsat ! "  ,T aat 

6390i  COMPUTE  NATURAL-CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 
6400'  FOR  UNENHANCED  ENO(S) 

6410  Hbar-190 

6420  Fa-< Hbar*P/<Keu«A  )  )* ,5»Lu 
6430  Tanh-FNTanh( Fa  ) 

64«0  Tneta-ThataP«Tanh/Fa 

6450  X*«<  9.81 •Bata»ThataO»Do*3»Tann/(Fa»Nj» Alpha  )  >* . 166667 

6460  Yy-<  !♦(  . 559/Pr  )*(  9/16  > )‘(  8/27) 

6470  Hbarc«K/Do*< .6+.3B7*Xx/Yy  )*2 
6480  IF  ABSI ( Hbar-Hbarc >/Hbarc >>.001  THEN 
6490  Hbar-< Hbar+Hbarc  )• .5 
6500  GOTO  6420 

6510  ENO  IF 

6511  PRINT  “NAT  CONV  h:“,Hbar 
65201 

65301  COMPUTE  HEAT  LOSS  RATE  THR0U6H  UNENHANCED  ENO(S> 

6540  01-( Hbar»P»Keu»A )" .5«Thatab*Tanh 

6550  Qc-Q-2»Q1 

6560  As-PI«02«L 
6570' 

6580'  COMPUTE  ACTUAL  HEAT  FLUX  ANO  80ILIN6  COEFFICIENT 

6590  Odp-Oe/Aa 

6600  Htuba-Qdp/Thatad 

6610  C»f-(Cp»Thatab/Hfg  )/<  Qdp/< Mu*Hfg  )•( .014/(9.81 •Rho  )“ .5  )*( 1 /3.  >«Pr" 1 .7 ) 
6620> 

6630'  RECORD  TIME  OF  DATA  TAKIN6 
6640  IF  1—0  THEN 
6650  OUTPUT  709i “TO* 

6660  ENTER  709iTold* 

6678  ENO  IF 
66801 

6690'  OUTPUT  DATA  TO  PRINTER 
6700  PRINTER  IS  701 
6710  IF  Iov-0  THEN 


672®  PRINT 

673®'  PRINT  U51N6  *10X,**Oata  Set  Nunner  •  *  *  ,000  ,2x  .  •  'Bulk  Oil  i  •  *",Q0.0,5X,i 
4A*  i  J  ,Bep  .Told® 

6731  PRINT  USIN6  *10X,**Oata  Set  Nunber  -  ••  ,000  ,2X  .  •  ‘Bulk  Oil  t  -  ", DD.D'iJ.B 


op 

6732  PRINT*  TIME : *  ,T1ME*( TIMEDATE > 

674®  IF  Ihn-0  THEN 

6750  PRINT  OSINS  * I0X . *  *TC  No;  1  2  3  4  5  S  7 

8*  •  • 


6760  PRINT  US I NS 
7 ) ,T( 8  ) 

677®  PRINT  US I  NS 
6780  PRINT  USINS 
Ld2 .Tv ,Psat .Tsunp 
6790  PRINT  USINS 
6800  PRINT  USINS 
6810  ELSE 
SS20  PRINT  USINS 
Ho*** 


*  10X  ,**Tenp  :**  ,8(  IX  .HDD. 00  )*  ;T<  I  )  ,T(2  )  ,T(3 >  ,T( 4  )  ,T( 5  )  .T( 6  >  ,T( 

“10X,*“  T-a  Tiiod  T 1 iad2  Tvepr  Pset  T>unp*** 

*10X  , 2 ( MOD . 00 , 1 X  >  , 1 X .HDD. DO .IX  ,2<  IX .MOD. 00) ,2X .MOD.D* ;TW .Tld ,T 

*10X.”  Thetab  Htube  Odp"** 

*  t  OX .MD0.3D.1X .M2. 3DE, 1 X .MZ.30E* ; Thetab  .Htube  ,Qdo 

*I0X.”  Fns  Uw  T sat  Tinl  Torop  Thetab  a  Uo 


6830  PRINT  USINS  *10X  .4(20. 00, IX  )  .2. 30. IX  ,00. 00. IX  ,3(MZ.30E.1X  >*iFns  ,U-.Teat  ,Ti 

nlet  , Tdrop .Thetab ,Qdp , Uo .Ho 

6840  END  IF 

685®  ENO  IP 

6860  IF  lov-l  THEN 

6870  if  j-i  then 

6880  PRINT 

6890  IF  Ihn-0  THEN 

6900  PRINT  USINS  *I0X,**  RUN  No  OilX  T»at  Htube  Qdo  Thetab*** 

6910  ELSE 

6920  PRINT  USINS  *I0X,**  FMS  OIL*  TSAT  HTUBE  QDP  THETAB'  *  * 

6930  END  IF 

6940  ENO  IF 

6950  IF  Ihn-0  THEN 

6960  PRINT  USIN6  *  1 2X  .30 ,4X  ,00  ,2X  ,M00 .DO  ,3(  1 X  ,MZ . 3DE  >*  i  J  .Bop  .Tsat  .Htube  .Qdp  .The 
tab 


6970  ELSE 

6980  PRINT  USING  *  1  2X .30 ,4X ,00 ,2X , MOO .00 ,3(1 X ,MZ . 3DE  ) * iFn*  .Bop  ,T»at  .Htube  .Odp  ,T 
netao 

6990  END  IF 
7000  ENO  IF 
7010  IF  In-0  THEN 
7020  SEEP 

7030  INPUT  *0K  TO  STORE  THIS  DATA  SET  (1 -Y ,0-N >7*  ,0k 
7040  ENO  IF 

7050  IF  Ok-I  OR  In- 1  THEN  J-J  +  I 
70S0  IF  Ok-I  AND  In-0  THEN 

7070  IF  Ihn-0  THEN  OUTPUT  *Fi lei (Bop  , To  Id*  ,Enf< •  1  . Ur . Ir 
7080  IF  Ihn-I  THEN  OUTPUT  OFilel  <Bop  .Told*  ,Enf<  o  ,Fm 
7090  END  IF 

7100  IF  Iuf-I  THEN  OUTPUT  0Uf llaiU- ,Uo 
7110  IF  Ire-1  THEN  OUTPUT  *RefileiFns  , Re- 

7120  IF  (In-1  OR  Ok-I)  ANO  Iplot-1  THEN  OUTPUT  0Plot iQdp  .Thetab 
7130  IF  In-0  THEN 
7140  BEEP 

715®  INPUT  'WILL  THERE  BE  ANOTHER  RUN  ( i -y  ,0-N >?■  ,6o_on 
7160  Nrun-J 

7170  IF  6o_on-0  THEN  7220 
7180  IF  6o_on<>0  THEN  Repeat 
7190  ELSE 

7200  IF  J<Nrun*1  THEN  Repeat 
7210  ENO  IF 
722®  IF  In-0  THEN 
723®  BEEP 

7240  PRINT  USIN6  *  1  OX  ,* "NOTE :  **,22.**  data  run*  -ere  stored  in  file  **.10A*iJ- 
1  ,02_f ile* 

725®  ASSI6N  OF  lie  I  TO  • 
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7260  OUTPUT  #F i  ie2 iNrun- 1 

7270  ASSIGN  CF  1 1  a  1  TO  01_filei 

7280  ENTER  0FilBliOatsi.Ldtc1.Ldtc2.Itt 

'OSC  OUTPUT  0FilB2iOatai.Ldtet.Ldtc2.Itt 

7300  FOR  I- 1  TO  Nrun-I 

7310  IF  Ihn-0  THEN 

7320  ENTER  0Fi lal iBoe  .Toldi  .Enf ( •  1  .Ur  ,Ir 
7330  OUTPUT  #Fila2iSop  , Toldi  .Eiif  <  •  )  .Ur  ,Ir 
73*0  ELSE 

7350  ENTER  #F  i  la  1  1 8op  .Toldi  .Enf  <•)  .Fns 
7360  OUTPUT  #Fila2i8op. Toldi  ,E«f<*>,F«s 
7370  ENO  IF 
7380  NEXT  1 

7390  ASSI6N  #F 1 1 a  I  TO  ♦ 

7*00  PURSE  * DUMMY* 

7*10  ENO  IF 
7*20  BEEP 
7*30  PRINT 

7**0  IF  Ip  lot* 1  THEN  PRINT  USIN6  *10X.**NOTE:  "*.22.**  X-Y  pairs  uere  stored 
plot  data  file  *  * . I0A* i J-l  .P_f ilei 
7*50  ASSI6N  0Fila2  TO  • 

7*60  ASSIGN  iPlot  TO  • 

7*70  IF  Iuf-I  THEN  ASSI6N  iUfile  TO  • 

7*80  IF  Ira-1  THEN  ASSIGN  iRefile  TO  » 

7490  CALL  Stats 
7500  BEEP 

7510  INPUT  ‘LIKE  TO  PLOT  DATA  ( 1 -Y ,0-N  )’ *  ,0k 
7520  IF  Ok- I  THEN  CALL  Plot 
7530  SUBENO 
75*01 

7550'  CURUE  FITS  OF  PROPERTY  FUNCTIONS 
7560  OEF  FNKcu(T) 

7570'  OFHC  COPPER  250  TO  300  X 
7580  Tk-T+273. IS  'C  TO  K 

7590  K-43*-. ) 12»Tk  . 

7600  RETURN  K 
7610  FNENO 

7620  OEF  FNMu(T) 

7621  COM  /Iprop/  I  ft 

7622  IF  Ift-0  THEN 

7630'  170  TO  360  K  CURUE  FIT  OF  UISCQUSITY 

7640  Tk-T+273. 15  1 C  TO  K 

7650  Hu-EXP<  -* .  *636*<  101  1  ,*7/Tk  1  Ml  .0E-3 

765)  END  IF 

7652  IF  Ift-1  THEN 

7653'  223  TO  373  K 

7657  Mu-3. *593-*. 2SE-2»T»3.9*85E-*»T-2-*. I 93E-6»T*3+2.0709E-8*T** 

7658  Mu-Mu* 1 . 0E-4 

7659  ENO  IF 

7660  RETURN  Mu 
7670  FNENO 

7600  OEF  FNCp(T) 

7681  COM  /Iprop/  I  ft 

7682  IF  Ift-0  THEN 

7690'  180  TO  *00  K  CURUE  FIT  OF  Cp 
7700  Tk-T*273. 15  'C  TO  K 

7710  Cp-. *0188*1 ,6S007E-3*Tk*1 .51 *9*E-o*Tk *2-6 .67B53E-10*Tk*3 

7720  Co-Co* 1 000 

7721  ENO  IF 

7722  IF  Ift-1  THEN 

7723'  223  TO  373  K  CURUE  FIT  Co<Rl2*> 

7725  Co- I .0542*2. I405E-3*T*I . 0709E-5*T*2-6 . *72 1 E-8*T*3-1 . *32*E-9*T*«** . i 36E-1 
T*5 

7726  Co-Co* 1000 

7727  ENO  IF 
7730  RETURN  Cp 
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77*0  FNENO 

7750  OEF  FNRhol T ) 

7751  COM  /loroe/  I  ft 

7752  IF  tft-0  THEN 

77S0  Tk-T*273. 15  !C  TO  K 
7770  X-l-< 1 ,8«Tk/7S3.95>  'K  TO  R 

7780  Ro-36. 32+61  . I464I**X*<  1  /3 )+!6 .*1 801 5+X  +  l 7 .476838+X* .S+1 .  1 1 9828+X ' 2 

7790  Ro-Ro/. 062428 

7791  END  IF 

7792  IF  Ift-I  THEN 

7793 1  223  TO  373  K  (R-124) 

7796  Ro«l*3*.8-2.8619»T-S.72S7E-3«T-2-7.28S2E-5»T‘3 

7797  END  IF 
7800  RETURN  Ro 
78 1 0  FNENO 

7820  DEF  FNPr(T) 

7830  Pr-FNCp( T >+FNMu( T )/FNK( T ) 

7840  RETURN  Pr 
7850  FNENO 

7860  OEF  FNK(T) 

7861  COM  /Iprop/  Ift 

7862  IF  Ift-0  THEN 
78701  T<360  K  WITH  T  IN  C 

7880  k-. 071 -.000261 *T 

7881  END  IF 

7882  IF  Ift-t  THEN 

7883!  223  TO  373  K  (R-124) 

7885  k*7.S191E-2-3.S436E-4*T-1 . 9S4SE-7.T‘2+3 . 1 835E-9+T'3 

7886  END  IF 
7890  RETURN  K 
7900  FNEND 

7910  OEF  FNTanh(X) 

7920  P-EXP(X) 

7930  0- I /P 

7940  Tanh-( P-Q  )/( P+Q  1 
7950  RETURN  Tanh 
7960  FNENO 
7970  OEF  FNTv»v(U) 

7980  COM  /Cc/  C( 7 )  ,leai 
7990  T-C( 0  ) 

8000  FOR  I-l  TO  7 
8010  T-T-Ct I >«U* I 
8020  NEXT  I 
8030  IF  lcal-1  THEN 

8040  T-T-6 .7422934£-2+T»< 9. 0277043E-3-T+( -9.325991 7E-5 > > 

8050  ELSE 

8060  T-T *8 . 626897E-2+T  +  1 3 . 76 1 99E-3-T+S.0689259E-5  ) 

8070  ENO  IF 
8080  RETURN  T 
8090  FNENO 
8100  OEF  FNBata(T) 

8110  Rop-FNRho<T+. I  > 

8120  Ron-FNRho<  T-.  1  ) 

81 30  Bata— 2/ ( Rop+Row  )•(  Rop-Ro* )/ .  2 
8140  RETURN  Seta 
8 1 50  FNEND 

8160  OEF  FNHfg(T ) 

8161  COM  /Iprop/  Ift 

8162  IF  Ift-0  THEN 

8170  Hfg-1 . 3741 344E+5-T+1 3.309436  I E+2+T  + 1 .2165143) 

8171  END  IF 

8172  IF  Ift-1  THEN 

8173'  223  TO  373  K  CURVE  FIT  <RI24) 

8176  Hfg-159  93- . 4609+T-1 .4S8E-3+T-2- 1 .3715E-S+T-3 

8177  Hf g-Hf g+ 1 000 

8178  ENO  IF 
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8188  RETURN  Hfg 

8 l 88  FNENO 

8280  OEF  FNP»at(Te  ) 

8201  COM  /Iprop/  Ift 

8202  IF  Ift-0  THEN 

8210'  0  TO  80  dag  F  CURVE  FIT  OF  Paat 
8220  Tf-i.8«TC*32 

8230  Pa-5.94SS2SfTf*< . 153S2082*Tf»< I .4840963E-3-Tf.9.6lS067IE-6 ) ) 

8240  Pa-Pa-14.7 

8241  ENO  IF 

8242  IF  Ift-I  THEN 

8243'  223  TO  373  K  1R-1241 

8245  Pa-( 1 62. 68-5. 946«Tc*9. 2081 E-2*Te '2+6. 9359E-4*Tc'3 )• . 1 4504 

8247  Pg-Pa- 14.7 

8248  ENO  IF 

82S0  IF  Pg>0  THEN  '  ♦-PSI6,—  in  H0 

8260  Psat-Pg 
8270  ELSE 

8280  Psat-Pg* 2 9. 92/ 14.7 
8290  ENO  IF 
8300  RETURN  P.at 
8310  FNENO 

8320  OEF  FNHsnoothl X  ,Bop  ,l*at  ) 

8330  DIM  A<5>  ,B(S)  ,C(5>  ,0(5) 

8340  OATA  .20526 ..25322  ..319048  ..S5322,. 79909,1 .00258 
8350  DATA  .74515. .72992. .73189. .71225, .68472. .64197 
8360  OATA  . 41092, .17726. .25142. .54806,. 81916. 1 .084S 
8370  DATA  .71403, .72913. .72565, .696691  ,.665867. .61889 
8380  READ  A(  •  )  ,B(  •  )  ,C(  •  >  ,0<  •  > 

8390  IF  Boo<6  THEN  I-Bop 
8400  IF  Bop-6  THEN  1-4 
8410  IF  Bod-10  THEN  I-S 
8420  IF  Isat-I  THEN 
8430  Hs"EXP<  A( I  )*8< I  )-L06<  X  >  I 
8440  ELSE 

8450  Ha-EXP<C< I  >*0< I  )*L0G(X  ) ) 

8460  END  IF 
8470  RETURN  Ha 
8480  FNENO 
8490  OEFFNPoly(X) 

8500  COM  /Coly/  A( 1 0 . 10  )  ,C< 10  )  ,8( 5  1  ,Npp  .Iornt  ,0po  ,  I  log  ,1 f n  ,  I  join .Njoin 
8510  Xl-X 
8520  Poly-B(0) 

8530  FOR  1-1  TO  Nop 

8540  IF  Ilog-I  THEN  XI-L061X) 

8550  Poly-Poly*B( I >*Xl  *  I 
8560  NEXT  I 

8570  IF  l log- I  THEN  Poly-EXP( Po ly ) 

8580  RETURN  Poly 
8590  FNENO 
8600  SU8  Poly 

8610  DIM  R<  10), S<  10),Sy(  12),S*(  !2).X*(  1001. Yy<  100) 

8620  COM  /Cply/  A(  10 , 10  )  ,C<  1 0  )  ,8<  5  )  ,N  ,lprnt  ,0po  ,  I  log  .1  f  n  .  i  join  .Njom 
8630  COM  /Xxyy/  Xp( 25  )  ,Yp< 25  ) 

8640  FOR  1-0  TO  4 
8650  B< I >-0 
8660  NEXT  I 
8670  BEEP 

8680  INPUT  'SELECT  <  0-FIlE , 1 -KEYBOARD , 2-PR06RAM  > '  ,  In 
8690  In-In*l 
8700  BEEP 

8710  INPUT  'ENTER  NUMBER  OF  X-Y  PAIRS'. Np 
8720  IF  In-l  THEN 
8730  BEEP 

8740  INPUT  'ENTER  OATA  FILE  NAME' ,0_f i la» 

8750  BEEP 
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8760  INPUT  "LIKE  to  EXCLUDE  OATA  PAIRS  ( I -y ,0-N  )’ '  ,  leg 
8770  IF  lad- I  THEN 

8780  SEEP 

8788  INPUT  'ENTER  NUMBER  OF  PAIRS  TO  BE  EXCLUDED*  ,Ipa* 
8800  END  IF 

8810  AS5I6N  #F  i  la  TO  0_fila8 
8820  ELSE 

8830  BEEP 

8840  INPUT  'UANT  TO  CREATE  A  DATA  FILE  <1 -Y  ,0-N  >7'  ,Yas 
8850  IF  Yaa-l  THEN 

8860  BEEP 

8870  INPUT  '6IVE  A  NAME  FOR  DATA  FILE* ,0_f i lat 
8880  CREATE  BOAT  0_fila«,5 
8890  ASSIGN  IF  i  la  TO  D_fila« 

8900  ENO  IF 
8910  ENO  IF 
8920  BEEP 

8930  INPUT  'ENTER  THE  ORDER  OF  POLYNOMIAL'  .N 

89*0  FOR  1-0  TO  N>2 

8950  Sy<l>-0 

8960  Sx<  I  >-0 

8970  NEXT  I 

8980  IF  lad- I  ANO  In- I  THEN 

8990  FOR  1-1  TO  Ipax 

9000  ENTER  0F i la i X  ,Y 

9010  NEXT  I 

9020  ENO  IF 

9030  FOR  1-1  TO  Np 

9040  IF  In-1  THEN 

9050  IF  Opo-2  THEN  ENTER  0FilaiX,Y 

9060  IF  Opo<2  THEN  ENTER  0FilaiY,X 

9070  IF  Opo-1  THEN  Y-Y/X 

9080  IF  I  log- 1  THEN 

9090  IF  Opo-2  THEN  Xt-X/Y 

9100  X-L061X) 

9110  IF  Opo-2  THEN  Y-LOG(Xt) 

9120  IF  Opo<2  THEN  Y-l06<Y) 

9130  ENO  IF 
91*0  ENO  IF 
9150  IF  In-2  THEN 
9160  BEEP 

9170  INPUT  'ENTER  NEXT  X-Y  PAIR'.X.Y 

9180  IF  Yaa-I  THEN  OUTPUT  0FileiX,Y 

9190  ENO  IF 

9200  IF  In<3  THEN 

9210  Xx<I>-X 

9220  Yy(  I  )-Y 

9230  ELSE 

92*0  X-Xp< I-i  > 

9250  Y-YptI-1) 

9260  ENO  IF 
9270  Rt  0  l-Y 
9280  $y(  0  )-Sy(  0  )*Y 

9290  S(  1  )-X 
9300  S*< I  >-Sx<  I  >*X 

9310  FOR  J-1  TO  N 
9320  RC J >-RC J-1  )*X 
9330  Sy<  J  >-Sy<  J  1*RU  ) 

93*0  NEXT  J 

9350  FOR  J-2  TO  N«2 

9380  St  J  »-S( J-1  >•* 

9370  Sx(J)-S*(J)*S(J) 

9380  NEXT  j 
9390  NEXT  I 

9*00  IF  Yaa-1  AND  In-2  THEN 
9410  SEEP 
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9420  eR  INT  USIN6  -12X.UL,”  *-»  same  -ere  storeo  i-  'ne  .  -  «jr  ;No,u_me* 
5430  END  IF 

9440  Sx<0>-Np 
9450  FOR  1-0  TO  N 
9450  C(I)-Syll) 

9470  FOR  J-0  TO  N 
9460  A<  I  ,  J  )-Sx< I»J  > 

9490  NEXT  J 
9500  NEXT  I 
9510  FOR  1-0  TO  N-l 
9520  CALL  Divide! I  ) 

9530  CALL  Subtract! I- 1  ) 

9540  NEXT  I 

9550  0<N  >-C(  N  )/At  N  ,N  ) 

9560  FOR  1-0  TO  N-l 
9570  e<N-l-I>«C<N-!-I> 

9580  FOR  J-0  TO  I 

9590  8  <  N- 1  - 1  )-8(  N- 1  -  I  >-A<  N-l-I  ,N-J>»8tN-J  ) 

9600  NEXT  J 

9610  BIN- I -I l-8( N- I -I )/A(N-l-I .N-l -I > 

9620  NEXT  I 

9630  'PRINTER  IS  701 

9640  'PRINT  8( • > 

9650  'PRINTER  IS  705 
9660  IF  Iprnt-0  THEN 

9670  PRINT  USING  *  I  2X  “EXPONENT  COEFFICIENT’” 

9660  FOR  1-0  TO  N 

9690  PRINT  USIN6  *  I  5X  ,00  ,SX  ,HO.  70E*  1 1  ,B<  I  > 

9700  NEXT  I 
9710  PRINT  ’  ’ 

9720  PRINT  USIN6  ’ I 2X  ’DATA  POINT  X  Y  Y( CALCULATED  )  OISCR 

EPANCY”’ 

9730  FOR  I- I  TO  Np 
9" 40  Yc»e<0) 

9750  FOR  J- I  TO  N 
9760  Ye«Yc*8< J  )»X«( I  )*J 
9770  NEXT  J 
9780  0-Yy<  I  >-Yc 

9790  PRINT  USIN6  ’  1  5X  ,30  ,4X  ,  4  (  HD .  5DE  ,IX)'il  ,Xx<  I  )  ,Yy(I  ),Yc  ,0 

9800  NEXT  I 

9810  END  IF 

9620  ASSI6N  0FiIe  TO  • 

9830  SU8EN0 
9840  SUB  Divide' M  ) 

9950  COM  / Coly/  At  1 0 , 1 0  >  ,C<  1 0  >  ,B<  5  >  ,N  .Iprnt  ,0po  ,1  log  .  I  fn  ,1  join  .Njoin 
9860  FOR  I-M  TO  N 
9070  Ao-AU  ,H> 

9880  FOR  J-M  TO  N 

9890  A(  I  ,  J  >— Ai  I  ,  J  )/Ao 

9900  NEXT  J 

9910  C( I > — C ( I  >/Ao 

9920  NEXT  I 

9930  SUBENO 

9940  SUB  Subtree t ( K  ) 

9950  COM  /Cply/  At  1 0 , 10  )  ,C<  10  >  ,B<  5  '  ,N  .Iprnt  ,0po  ,  I  log  .1  f  n  ,1  join  .Njoin 

9960  FOR  I-K  TO  N 

9970  FOR  J-K-l  TO  N 

9960  A(  I  ,J  )-A(K-l  ,J  >  —  A(  I  ,J  ) 

9990  NEXT  J 

10000  Cl  I  >-C( K - 1  ) — C <  I  ) 

10010  NEXT  I 
10020  SUBENO 
10030  SUB  PI  in 

10040  COM  /Cply/  At  10 , 1 0  )  ,L< 10 >  ,6(5 )  ,N , Iprnt  ,0po  ,  I  log  .Ifn  .  I  join .Njoin 
10050  COM  /Xxyy /  X«<  25  )  ,Yy( 25  ) 

10060  PRINTER  IS  70S 
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10070  SEEP 

10080  INPUT  'WANT  TO  PLOT  Uo  v»  U«Y  < 1 -Y ,0-N ) '  , luo 
10030  IF  Iuo-0  THEN 
10100  BEEP 

10110  INPUT  'SELECT  (0-h/hBI  sana  tub« , I -h(HF  )/h( an  )*  .Irt 
10120  SEEP 

10130  INPUT  'SELECT  h/h  RATIO  11 -FILE ,0-COHPUTED >' .Ihrat 
10140  IF  Ihrat-0  THEN 
101 S0  8EEP 

101S0  INPUT  "WHICH  Tsat  (  1 -6 .7  ,0— 2 .2  >*  .laat 

10170  ENO  IF 

10160  Xnin-0 

10190  Xnax- 10 

10200  Xstap-2 

10210  IF  Irt-0  THEN 

10220  Ynin-0 

10230  Ynax-1 .4 

10240  Yatao-.2 

10250  ELSE 

10260  Ynin-0 

10270  Ynax-1 5 

10260  Ystap-S 

10230  ENO  IF 

10300  ELSE 

10310  Odo-2 

10320  Ynin-0 

10330  Ynax-1 2 

10340  Ystep-3 

10350  Xnin-0 

10360  Xnax-4 

10370  Xstap-1 

10380  ENO  IF 

10390  IF  Ihrat- 1  THEN 

10400  Ynin-0 

10410  Ynax-1 8 

10420  Ystap-3 

10430  Xnm-0 

10440  Xnax-9 

10450  Xstap-2 

10460  ENO  IF 

10470  BEEP 

10480  PRINT  * INiSP l (IP  2300.2200.8300,66001' 

10490  PRINT  'SC  0 . 100 ,0 , 1 00iTL  2.0i* 

10500  Sfx-1 00/ ( Xnax-Xnin  1 
10S10  Sfy-l00/(  Ynax-Ymn  > 

10520  PRINT  'PU  0.0  PO' 

10530  FOR  Xa-Xnin  TO  Xnax  STEP  Xatao 

10540  X-( Xa-Xnin  )«Sfx 

10550  PRINT  'PA'iX,',0i  XT i ' 

10560  NEXT  Xa 

10570  PRINT  'PA  100.0iPUi' 

10580  PRINT  'PU  PA  0.0  PO” 

10590  FOR  Ya-Yntn  TO  Ynax  STEP  Ystap 
10600  Y-t Ya-Yntn )»Sfy 
10610  PRINT  'PA  0  , ' t Y  ,  'YT* 

10620  NEXT  Ya 

10630  PRINT  *PA  0.100  TL  0  2' 

10640  FOR  Xa-Xnin  TO  Xnax  STEP  Xitap 

10650  X-( Xa-Xnin >»Sf x 

10660  PRINT  'PA' iX  ,'  .1001  XT' 

10670  NEXT  Xa 

10680  PRINT  'PA  100,100  PU  PA  100,0  PO’ 

10690  FOR  Ya-Ymn  TO  Ynax  STEP  Ystap 

10700  Y-(  Ya-Ymn  »•  S P y 

10710  PRINT  'PO  PA  100,*  ,Y  , ' YT  * 

10720  NEXT  Ya 
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0730  PRINT  "PA  100,100  PU" 

07*0  PRINT  "PR  0,-;  ift  1.5.2' 

0750  FOR  Xa-Xnm  TO  Xnax  STEP  Xatao 
0760  X-<  Xa-Xmn  )«Sfx 
0770  PRINT  "PR" i X , ’  ,01  " 

0780  IF  luo*0  THEN  PRINT  "CP  -2  .-1 |LB" iXai * " 

0790  IF  Iuo-I  THEN  PRINT  "CP  -1 .5  ,*>  iLB"  -.Xai  "  " 

0800  NEXT  Xa 

0810  PRINT  "PU  PR  0.0" 

0820  FOR  Ya-Ynm  TO  Ynax  STEP  Yateo 
0830  IF  RBSIYaXI  .E-S  THEN  Ya-0 
08*0  Y-<  Ya-Ymn  )»Sfy 
0850  PRINT  "PR  0  , " i Y  , " " 

0860  IF  Iuo-0  THEN  PRINT  "CP  -*  . 2S-.LB"  I  Ya I  * " 

0870  IF  Iuo-I  THEN  PRINT  "CP  -3  . 25 1 LB" ; Ya  I " " 

0880  NEXT  Ya 

0890  Xlabal*-*Oil  Percent' 

0900  IF  Iuo-0  THEN 

0910  IF  Irt-0  THEN 

0920  Ylabel*-"h/h0X" 

0930  ELSE 

09*0  Ylabelf-'h/hanooth" 

0950  END  IF 

0960  PRINT  "SR  1 .5 ,2 iPU  PR  50. -10  CP* i -LEN( X Labels >/2i "0iLB* iXlabelti * ' 

0970  PRINT  "PR  -11,50  CP  0  ,*  t-LEN(Y label* >/2«S/6i "01  0 . 1 iLB* i YiabelSi ' * 

0980  PRINT  "CP  0,0" 

0990  ELSE 

1000  PRINT  "SP0iSP2" 

1010  PRINT  "SR  1 .2 ,2.*|PU  PR  -8 ,35|0I  0,I|LBU|PR  1  ,0.5iLBoiPR  -1  ,0.5(LB  lkU/« 

1020  PRINT  "PR  -1  ,0.5iSR  l.1.5iLB2iSR  I.S.2|PR  .5,.5iLB.iPR  ,S.0iLBK>" 

1030  PRINT  "PR  *2  ,  —  1@i 01  1  ,0:LBU|PR  .*,-1iLBwiPR  1  .  .SiLBl «/s >* 

10*0  PRINT  "SP0iSP1" 

1050  END  IF 
1060  Ion-0 
1070  BEEP 

1080  INPUT  "URNT  TO  PLOT  ORTR  FROH  R  PILE  ( 1 -Y ,0-N  >7"  ,Okp 

1090  Icn-0 

1100  IF  Oko-1  THEN 

1110  BEEP 

1120  INPUT  "ENTER  THE  NRME  OF  THE  ORTR  FILE" ,D_file* 

1130'  IF  Iuo-0  THEN 
I  1 40  BEEP 

1150  INPUT  "SELECT  (0-LINERR,  1 -L06<  X  ,Y  >'  ,  I  log 
1160'  END  IF 

1170  RSSI6N  0Fi la  TO  0_fila» 

1 1 80  BEEP 

1190  INPUT  "ENTER  THE  BE6INNINS  RUN  NUMBER"  ,Hd 
1200  BEEP 

1210  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PRIRS  STORED" .Npaira 
1220  IF  Iuo-0  RNO  Ibrat-0  THEN 
1230  BEEP 

12*0  INPUT  "ENTER  OESIREO  HERT  FLUX',0 

1250  END  IF 

1260  BEEP 

1270  PRINTER  IS  1 

1280  PRINT  USIN6  "*X  ." "Select  a  symbol:""" 

1290  PRINT  US  I  NS  "*X,""1  star  2  Plus  HQn""" 

1300  PRINT  USIN6  "*X,""3  Circle  *  Souere""" 

'310  PRINT  USIN6  "*X , " "5  Ronoua " " * 

1320  PRINT  USIN6  "*X.""6  Rignt-side-uo  triengle"*" 

1330  PRINT  USIN6  "*X,""7  Ue-»ide-aoun  triengle""" 

13*0  INPUT  Sy* 

1350  PRINTER  IS  70S 
1380  PRINT  "PU  01" 

1370  IF  Syn-I  THEN  PRINT  "SM«" 
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M382  IF  Syn-2  THEN  PRINT  "SM*" 

11390  If  Syn-3  Th£N  rRiNi  bMo 

1 1 400  Nn-4 

11410  IF  I loo”'  THEN  Nn-1 
1)420  IF  Md> 1  THEN 
11430  FOR  I-t  TO  (Md-1  ) 

11440  ENTER  BFileiXa.Ya 
1 14S0  NEXT  I 
1  use  ENO  IF 

11470  IF  lhrat»0  THEN 
11480  01-0 

11490  IF  I  log* 1  THEN  0-L061Q) 

11500  ENO  IF 

11510  FOR  1-1  TO  Npaica 

11520  IF  Iuo-0  ANO  Ihret-0  THEN 

11530  ENTER  tFilaiXa  ,B( •  1 

11540  Ya-8(0> 

IIS50  FOR  K«l  TO  Nn 
1 1560  Ya-Ya+B< K >*0*K 
1 1570  NEXT  K 
1 1580  ENO  IF 

11590  IF  luo-1  OR  Ihrat-1  THEN 
1 1S00  ENTER  0FilaiXa,Ya 
11610  IF  Iuo*>  THEN  Ya-Ya/1000 
11620  END  IF 

11630  IF  Iuo-0  AND  lhrat-0  THEN 
11640  IF  I  log-1  THEN  Ya-EXP(Ya) 

11650  IF  I log-0  THEN  Ya-Ol/Ya 

11660  IF  Irt-0  THEN 

11670  IF  Xa-0  THEN 

1 1680  Yo-Ya 

11690  Ya-1 

1 1 700  ELSE 

11710  Ya-Ya/Yo 

1  1 720  ENO  IF 

11730  ELSE 

11740  Han-FNHanootmO.Xa.Isat  ) 

11750  Ya-Ya/Han 
11760  ENO  IF 
11770  ENO  IF 
11780  X*(I-1  )-Xa 
11790  Yy( 1-1  )-Ya 
11800  X*(  Xa-Xnm  )»Sf* 

1 181  Y-< Ya-Ynxn )»STy 

l 182 ^  IF  Syn>3  THEN  PRINT  "SM" 

11830  IF  Syn<4  THEN  PRINT  "SR  1. 4,2.4" 

11840  PRINT  "PA"  ,X  ,Y  ,*" 

11850  IF  Syn>3  THEN  PRINT  "SR  1.2. 1.6" 

11860  IF  Syn-4  THEN  PRINT  "UC2 .4 .99 ,0 .-8 ,-4 ,0  ,0 .8  ,4  .0 .  i " 
11870  IF  Syn-5  THEN  PRINT  "UC3  ,0 .99  ,-3  .-S  .-3  .6  .3  .6  .3  .-6i  " 
11880  IF  Syn-6  THEN  PRINT  "UC0 .5. 3 .99 ,3 .-8 ,-6 ,0 .3 ,Bi " 
11890  IF  Syn«7  THEN  PRINT  "UC0  .-5 .3 ,99 ,-3 ,8 .6 ,0 ,-3  ,-8i " 
11900  NEXT  I 
11910  8EEP 

11920  ASSI&N  0Fila  TO  • 

11930  ENO  IF 
11940  PRINT  "PU  SH" 

1 1 950  BEEP 

11960  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  1 1 -Y  ,0-N >’ "  ,OXp 
11970  IF  OXp-l  THEN 
I  I  990  BEEP 
11990  PRINTER  IS  1 

12000  PRINT  USIN6  *4X  ," "Salaet  Una  type:""" 

12010  PRINT  USIN6  "6X  . " "0  Solid  Una*" 

12020  PRINT  USIN6  *6X,""l  Oaahad"" 

12030  PRINT  USIN6  *6X ,**2 ,  ,  ,5  Longer  line  -  dean"** 
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12040  INPUT  Ipn 
ixdad  FAINTER  la  70S 
1 2060  SEEP 

12070  INPUT  “SELECT  ( 0-LlNEAR  ,  1 -L06C  X  ,Y  >  1“  ,Uoa 

12080  Iprnt«1 

12090  CALL  Poly 

12100  IF  luo-l  THEN 

12110  SEEP 

12120  INPUT  “DESIRE  TO  SET  X  Lower  and  Uppar  Limit  (  l-Y  ,0-N  >7“  ,I«lm 

12130  IF  Ixlim-0  THEN 

12140  Xll-0 

121S0  Xul-7 

12160  END  IF 

12170  IF  Ixlim-I  THEN 

12180  SEEP 

12190  INPUT  “ENTER  X  Lower  Limit'. XU 
12200  SEEP 

12210  INPUT  “ENTER  X  Upper  Limit*  .Xui 
12220  ENO  IF 
12230  END  IF 

12240  FOR  Xa-XIl  TO  Xul  STEP  Xstep/2S 
12250  Icn-Icn+1 
12260  Ya-FNPoly( Xa ) 

12270  IF  Iuo-l  THEN  Ya-Ya/1000 

12280  Y-<Ya-Ynm)«Sfy 

12290  X-( Xe-Xnin  >»Sfx 

12300  IF  Y<0  THEN  Y-0 

12310  IF  Y> 1 00  THEN  60T0  12410 

12320  Pu-0 

12330  IF  Ipn-1  THEN  Idf-Icn  MOO  2 
12340  IF  Ipn-2  THEN  Idf-Icn  NOD  4 
12350  IF  Ipn-3  THEN  Idf-Icn  MOO  8 
12360  IF  Ipn-4  THEN  Idf-Icn  MOO  16 
12370  IF  Ipn-5  THEN  Idf-Icn  MOD  32 
12380  IF  Idf-1  THEN  Pu-1 
12390  IF  Pu-0  THEN  PRINT  “PA“  .X  ,Y  ,*PO* 

12400  IF  Pu-I  THEN  PRINT  “PA“ ,X ,Y . “PU* 

12410  NEXT  Xa 
12420  PRINT  “PU“ 

12430  SOTO  11070 
12440  END  IF 
12450  BEEP 

12460  INPUT  -WANT  TO  QUIT  ( 1 -Y  ,0-N  )7“  ,  Iqui t 
12470  IF  Iouit-1  THEN  12490 
12480  SOTO  11070 
12490  PRINT  “PU  SP0“ 

12500  SUBENO 
12510  SUB  Stata 
12520  PRINTER  IS  701 
12530  J-0 
12540  K-0 
12550  BEEP 

12560  IF  Ip  lot  —  1  THEN  ASSIGN  0Fila  TO  P  file* 

12570  BEEP 

12580  INPUT  “LAST  RUN  No7(0-QUIT  1“  ,Nn 

1 2590  IF  Nn-0  THEN  12950 

12600  Nn-Nn-J 

12610  Sx-0 

12620  Sy-0 

12630  Si-0 

12640  Sxa-0 

1 2650  Sys-0 

12660  Sz*-0 

12670  FOR  I- I  TO  Nn 

12680  J-J+l 

12690  ENTER  IFileiQ.T 
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2780  h-Q/T 
2718  Sx-Sx+G 
2720  Sxa-Sx*+0*2 
2730  Sy-Sy»T 
2740  5y»*Sy»+T*2 
27S0  Sz-Sz+H 
2760  Szs-Szj*h*2 
2770  NEXT  I 
2780  Qave-Sx/Nn 
2790  Tava-Sy/Nn 
2800  Have-Sz/Nn 

2810  Sdavq-SQR< A8S( <Nn»Sxs-Sx*2 >/<Nn«<Nn-l  )))) 

2820  Sdavt-SQR< A8S< <Nn»Sys-Sy*2 >/<Nn»<Nn-1 ) ) ) ) 

2830  Sdavh-SQR(A8S< <Nn»Sz»-Sz*2 >/(Nn*(Nn-1 )>)> 

2840  Sh-l00»Sdavh/Hava 
28S0  So-i00«Sdavq/Oave 
2860  St-l00»Sdevt/Tava 
2870  IF  K-1  THEN  12930 
2880  PRINT 

2890  PRINT  USIN6  '1  IX  .'•"ORTA  FILE:  "  ,I4A*  tFile* 

2900  PRINT 

2910  PRINT  USIN6  MIX, "RUN  Htuba  SdavH  Qdp  SdevQ  Thetab  SdavT" 

2920  K-1 

2930  PRINT  USIN6  *  I  I X ,00 ,2< 2X ,D. 3DE , lx ,30 .2D ) ,2X  ,DD .30 , 1 X  ,3D . 2D” i J  .Hava  ,Sh  ,Qave 

Sq  ,Tava  ,St 

2940  60T0  12570 

2950  ASSISN  0Fi lai  TO  - 

2960  PRINTER  IS  1 

2970  SUBENO 

2980  SUB  Coaf 

2990  COM  /Col//  A(  10 , 1 0 ) ,C< 10 > .8(5 ) ,N .Iqpnt ,0oo ,I1oq .1 fn .1  join .Njoin 
3000  BEEP 

3010  INPUT  '6IVE  A  NAME  FOR  CROSS-PLOT  FILE*  ,Cpf» 

3020  CREATE  BOAT  CpfS.2 
3030  ASSI6N  BFile  TO  Cpf« 

3040  BEEP 

3050  INPUT  ’SELECT  <  0-LlNEAR , 1 -L06<  X  ,Y  >  >*  ,1  log 
3060  BEEP 

3070  INPUT  -ENTER  OIL  PERCENT  (- 1 -STOP >*  .Boo 

3080  IF  Bop<0  THEN  13120 

3090  CALL  Poly 

3100  OUTPUT  BFilaiSop  ,B< •  ) 

3110  60T0  13060 
3120  ASSI6N  tFila  TO  • 

3130  SUBENO 

3140  SUB  Ullsonl Cf ,Ci > 

3150  COM  /Uil/  02  ,0i  ,0o  ,L  .Lu  ,Keu 
3160  OIM  Enf  < 1 2  > 

31701  ULISON  PLOT  SUBROUTINE  DETERMINE  CF  AND  Cl 
3-80  BEEP 

3190  INPUT  ‘ENTER  DATA  FILE  NAME"  ,Fila0 

3200  BEEP 

3210  PRINTER  IS  1 

3220  PRINT  USIN6  *4X. "Salad  option:*" 

3230  PRINT  USIN6  *4X,"  0  Uary  Cf  and  Ci*" 

3240  PRINT  USIN6  *4X,**  1  Fix  Cf  Vary  Cl*** 

3250  PRINT  USIN6  *4X,"  2  Uary  Cf  Fix  Ci*" 

3260  INPUT  ‘ENTER  OPTION*  .Iefix 

3270  PRINTER  IS  701 

3280  IF  Iefix-0  THEN  13320 

3290  IF  Icf ix>0  THEN  BEEP 

3300  IF  Icfix-1  THEN  INPUT  ’ENTER  Cf.Caf 

3310  IF  Iefix-2  THEN  INPUT  ‘ENTER  CI’.Ci 

3320  PRINTER  IS  I 

3330  INPUT  ’Want  To  Uary  Coaf f “»< 1-7 ,0-N )* .Iceoaf 
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13340  IF  Iccoef«l  THEN  INPUT  "ENTER  COEFF"  R 

13350  PRINTER  IS  701 

13360  IF  Icfix-0  OR  Icfix-2  THEN  Cfa-.004 

13370  IF  tefix«l  THEN  Cfa-Csf 

13380  Cia-Ci 

13330  Xn-.8 

13400  Fr.,3 

i34i0  :j-e 

13420  Rr-3. 

13430  IF  Iccoaf-1  THEN  Rr-R 

13440  PRINTER  IS  I 

134S0  PRINT  Oo.Ol.Kcu 

13460  ASSI6N  CFtli  TO  Filel 

13470  ENTER  BFileiNrun ,Oate*  .Ldtcl .Ldte2  .Itt 

13480  Ru“0o*L06( Do/Oi )/(2»Kcu  1 

13430  Sx-0 

13500  Sy-0 

13510  Sxy«0 

13520  Sx2»0 

13530  Sy2-0 

13540  FOR  I-l  TO  Nrun 

13550  ENTER  0FileiBoo  .Told*  .Enf ( •  1  ,F«a 

135601  CONVERT  EMF'S  TO  TEMPERATURE 

13570  FOR  J-t  TO  5 

13580  T( J  )«FNTv»v<  Enf  ( J  ) ) 

13530  NEXT  J 

13600  T»at-(T<  I  >+T<2)  >».5 
13610  Tavg“T(5) 

13620  Grad*37 . 9853+ . 104388*Tavg 

1 3630  Tdrop«Enf (7  )• 1  . E+G/<  1 0.«6rad  1 

13640  Tavgc«T( 5 >-Tdrop* .5 

13650  IF  ABSITavg-Tavge )>.0!  THEN 

13660  Tavg“(Tavg+Tavgc  )».S 

13670  SOTO  (3620 

13680  END  IF 

136301 

137001  Conpute  properties  of  water 
13710  KwFNKw(Tavg) 

13720  Muua*FNMuw(Tavg> 

13730  Cpw-FNCpw< Tavg  ) 

13740  Prw*FNPrw< Tawg ) 

13750  Rhow-FNRhowfTavg  1 
1 3760 ' 

13770'  Compute  properties  of  Freon-114 
1 3780  Lntd"Tdrop/L06( ( T( 5  >-T»at )/( T( S  )-Tdrop-T aat  ) ) 

13790  IF  Jj-0  THEN 
13800  Tw«Tsat+Fr»Cntd 
13810  Theteb"Ty-Tsat 
13820  4j-l 
13830  END  IF 
13840  Tf«( Tw+T»at  )• . 5 
13850  Rho-FNRho<Tf ) 

13860  Mu«FNMu<Tf  1 
13870  X-FNK(Tf) 

13880  Cp“FNCp( T  f  ) 

13890  8eta»FNBeta<  T  f  ) 

13300  Hfg-FNHfg<T»at  ) 

13310  Nl«Mu/Rho 
13320  Aloha»K/( Rho»Cp  ) 

13930  Pr«Nt/Alpha 
139401 

139501  Analyiti 
139601  COMPUTE  MDOT 
13970  A«PI»<  Do‘2-0i *2  )/* 

13980  P-PI*0o 

13990  Hdot«3.9657E-3+Fne*( 3.6 1 955E-3-Fns»l 8. 82006E-6-Fns»<  1 .23688E-7-Fna»4 . 3 1 89? 
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E- 1  ®  >  >  > 

4006  OMdot •Cp««TaroB 

46101  COMPOTE  NATURAL-CONUECT I UE  HEAT-TRANSFER  COEFFICIENT 
4020)  FOR  UNENHANCEO  ENO(S) 

4030  Hbaral90 

4040  Fea(Hbar»P/(Kcu*A  )  )* .5*Lu 
4050  TanhaFNTanh<Fa > 

43S5  ThetaaThatab*Tann/'Fe 

4070  Xxa<  9.81 «0at4» That ab»0o'3»Tanb/( F#«Nj  »A]pna  /  >* ,  ISB6E7 
4080  ¥y“( !♦<  .S59/Pr )*( 9/ 1 6 >  )* ( 8/27  ) 

4090  HbarcaK/Do»< .S+.387»Xx/Yy >*2 

4100  IF  ABSI  (Hbar-Hbarc >/Hbar >> .001  THEN 

4110  Hbara< Hbar*Hbare >• .S 

4120  SOTO  14040 

4130  END  IF 

41401 

4 1 50 1  COMPUTE  HEAT  LOSS  RATE  THROUSH  UNENHANCEO  ENOS 
4160  Ql •( Hbar«P»Kcu»A  )" .S*ThetaoaTanh 
4170  Qc-Q-2«01 
4180  AsaPI«02»L 

41901  COMPUTE  ACTUAL  HEAT  FLUX 
4200  QdoaQc/Aa 

4210  IF  Icfix-0  OR  Ief ix> I  THEN  Caf»l /Cf * ( 1 . /Hr  ) 

4220  ThatabaCaf/C0*Hfg«<Qdp/<Mu»Hfg>«< .01 4/ <9. 81 *Rho>  >* . 5  > “ ( I /Rp >»Pr*1  .7 

4230  HoaQdp/Thatab 

4240  0««gaaHo/Cf 

4250  Uoa0/ ( PI *Oo*L*Lnt d  > 

4260  UwaMdot/(Rhow»PI*0i"2/4 1 
4270  RawaRhow*Uw«Oi/Muwa 
4280  TwiaTu+0»Ri»/(PI»0o*L  ) 

4290  GanaaKw/0i*Rew*  .8»Pru*(  1  /3 .  )•<  Mut*a/FNMu«i(  Tui )  )*.14 

4300'  PRINTER  IS  1 

4310  Yw(  1  ./Uo-Rw>a0naga 

4320  Xwa0nagaa0o/(6ana*0i  ) 

♦330  Sx-Sx+Xu 
4340  Sy-Sy+Yu 
4350  5xyaSxy+Yw»Xu 
4360  Sx2-Sx2+Xu»Xu 
♦370  Sy2-Sy2*Y*»Y« 

4380  NEXT  I 

4390  ASSI6N  #Fila  TO  • 

4400  Ma<  Sx»Sy-Nrun*Sxy  )/( Sx»Sx-NrunaS»2  > 

4410  Ca<  Sy-Sx»H)/Nrun 

4420  IF  Ieflx-0  OR  lcfix-3  OR  Iefix-4  THEN 

4430  Cica1/n 

4440  Cfc-I/C 

4450  END  IF 

4460  IF  Iefix-t  THEN 

4470  Cic-I/M 

4480  Cfe«Cf 

4490  ENO  IF 

4500  IF  Icfix-2  THEN 

4510  Cic-Cl 

4520  Cfc-l/C 

♦530  ENO  IF 

4540  IF  A8S< ( Ci-Cje  )/Cle  )> .001  OR  A8SI ( Cf-Cf c  )/Cf c >> .001  THEN 
4550  Cia<Ct*Cte  >».5 
4560  Cfa(Cf«Cfe  )».5 
4570  PRINTER  IS  1 

4580  PRINT  USIN6  "2X  ‘  Caf  a  *  *  ,MZ .30E  .2x  .* •  Ci  -  “  .M2.30E' iCaf .Ci 

4590  PRINTER  IS  701 

4600  60TO  13460 

4610  ENO  IF 

4620  PRINT 

4630  PRINTER  IS  701 

4640  PRINT  USIN6  *23X.”  Cf  Ci*** 
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14653  PRINT  USING  "8*  , " "ASSUMED  ' " ,M2 .30E  ,3X  .i:  .5DE" ; :ra  ,Zia 

14663  PRINT  USIN6  "8X  "CALCULATED  ' " ,MZ . 3DE .3* .M2 .30E" ; C ■> f .2; 

14670  PRINT 

1*680  Su«2"Sy2-2*N*Sxy-2«C*Sy*,H‘,2*Sx2'*2*H»C*5x+Nrun*C*2 
14690  PRINT  USIN6  '10X,”Suix  of  Squ«r«»  -  " ‘ ,Z .3DE" i Sum2 
1  *700  PRINT  USIN6  •  10X  ,”Co«ff  tei«nt  -  ”,O.ODO"iRp 
14710  SU8ENO 
1*720  OEF  FNMuu(  T  ) 

1*730  A«2*7.8/(Tf133.!S> 

1*7*0  Mu“2 . 4E-5* i 0* A 
14750  RETURN  Mu 
1*760  FNENO 
1*770  DEF  FNCpu(T) 

1*780  Cpu-*.21 1208S8-T«(2.26826E-3-T*(*.*236!E-S*2.7l*28E-7*T  ) ) 

1*790  RETURN  Cpu* 1 000 

1*800  FNENO 

1*810  OEF  FNRhoul T  1 

1*820  Ro«999.529*6*T*< .01 269-T*< 5.48251 3E-3-T*1 .23*i*7E-S)1 

14830  RETURN  Ro 

1*840  FNENO 

1*850  OEF  FNPru(T) 

I *860  Prw-FNCput T )*FNMuui T  )/FNKu( T  ) 

1*870  RETURN  Pru 
1*880  FNENO 
14890  OEF  FNKu(T) 

14900  X*<T*273. 151/273. IS 

1*910  Ku—  .92247*X*(2.839S-X*(  1 .8007-X*(  .S257?-.073***X  ) ) ) 

14920  RETURN  Kw 
1*930  FNEND 
1*940  SUB  Plot 

14950  COM  /Cply/  A<  10 .101  ,C<  10  1  ,8(5)  .Nop.Ip.-nt  .Opo  .Iloo.Ifn  .  Ijoin  .Njoin 
14960  DIM  B s  <  3  > 

1*970  INTE6ER  tl 
1*980  PRINTER  IS  I 
1*990  Idv-0 
15000  BEEP 

15010  INPUT  ’LIKE  DEFAULT  VALUES  FOR  PLOT  ( 1 -V .0-N )?• ,Idv 
1 5020  Ooo«0 
15030  BEEP 

150*0  PRINT  USIN6  "*X.”S«l«ct  Option:*” 

15050  PRINT  USIN6  *6X.**0  q  v«rsu»  d«ltt-T""" 

15060  PRINT  USIN6  *6X,”1  h  vtrsuJ  deltq-T”* 

15070  PRINT  USIN6  *6X.”2  h  v«r»u»  q*” 

15080  INPUT  Opo 
15090  BEEP 

15100  INPUT  "SELECT  UNITS  < 0-SI  . f -EN5LISH 1* . lun 
151 10  PRINTER  IS  705 
15120  IF  Idv0 1  THEN 
15130  BEEP 

151*0  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  X-AXIS"  ,Cx 
15150  BEEP 

15160  INPUT  "ENTER  NUMBER  OF  CYCLES  FOR  Y-AXIS", Cy 
15170  SEEP 

15180  INPUT  "ENTER  MIN  X-VALUE  (MULTIPLE  OF  l0)".Xnin 
15190  BEEP 

15200  INPUT  "ENTER  MIN  Y-UALUE  (MULTIPLE  OF  10)",Ynin 

15210  ELSE 

15220  IF  Opo-0  THEN 

15230  Cy-3 

152*0  Cx-3 

1 5250  Xnm«  .  1 

15260  Y*in»l00 

15270  ENO  IF 

15280  IF  Opo-1  THEN 

15290  Cy-2 

15300  Cx*2 
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1  S3 ' 8  Xnin-. t 
1 5320  rmn-100 
15330  END  IF 
15340  IF  Ooo-2  THEN 
15350  IF  Iun-0  THEN 
15360  Cy-3 
15370  Cx-2 
1S360  Xn in* 1 000 
15390  Ynin-100 
1S408  ELSE 
15410  Cy-3 
15420  Cx*3 
15430  Xmn-100 
15*40  Ynin-10 
1 5450  END  IF 
15460  ENO  IF 
15470  ENO  IF 
1 5480  BEEP 

15*90  PRINT  " INiSPI i IP  2300  ,2200 ,8300 ,68001 ’ 

15500  PRINT  ‘SC  0 , 1 00  ,0 , 1 00 1  T[_  2  ,0i  ' 

15510  Sfx-100/Cx 
15520  Sfy-100/Cy 
15530  BEEP 

155*0  INPUT  "WANT  TO  8Y-PASS  CASE’  <  I -Y  ,0-N  >•  .  Ibyts 
15550  IF  Ibyp-I  THEN  16790 
15560  PRINT  -PU  0,0  PO* 

1 5570  Nn*9 

15580  FOR  I- I  TO  Cx+I 
15590  X«t-Xmn«l0*<I-l  > 

15600  IF  I -Cx+I  THEN  Nn-1 

15610  FOR  J-1  TO  Nn 

15620  IF  J-1  THEN  PRINT  "TL  2  0' 

15630  IF  J-2  THEN  PRINT  ‘TL  t  0' 

15640  X*»X«t+J 

15650  X-L6T (Xa/Xxiinl+STx 

15660  PRINT  'PA'|X.',0|  XT  i ' 

15670  NEXT  J 
15680  NEXT  I 

15690  PRINT  -PR  100,0iPur 
15700  PRINT  “PU  PA  0,0  PO" 

15710  Nn-9 

15720  FOR  1-1  TO  Cy+1 

15730  Y»t-Ynin+ 1 0“ ( 1-1  1 

157*0  IF  I -Cy+1  THEN  Nn-1 

15750  FOR  J-1  TO  Nn 

15760  IF  J-1  THEN  PRINT  "TL  2  0" 

1S770  IF  J-2  THEN  PRINT  'TL  1  0" 

1S780  Y«-Y«t+J 

15790  Y-L6T (Y«/Yntn)«Sfy 

15800  PRINT  'PA  0  , ' i Y  , ' YT ' 

15810  NEXT  J 
15820  NEXT  I 

15830  PRINT  'PA  0,100  TL  0  2* 

<58*0  Nn-9 

15850  FOR  1-1  TO  Cx+1 

15860  X«t-X«iin+!0-<  1-1  1 

15870  IF  I-Cx+t  THEN  Nn-1 

15880  FOR  J-1  TO  Nn 

15890  IF  J-1  THEN  PRINT  'TL  0  2' 

15900  IF  J>l  THEN  PRINT  'TL  0  I* 

15910  X«-X«t+J 

15920  X-L6T1 X«/Xnin  )»Sfx 

15930  PRINT  *PA*  «X , * ,100i  XT' 

15940  NEXT  J 
15950  NEXT  I 

15960  PRINT  ‘PA  100,100  PU  PA  100,0  PO* 
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15970  Nn«S 

' 53S8  FOR  I-i  TO  Cy-1 

15990  Yat-Yntn.i0*<  I-l  1 

16000  IF  I-Cy+1  THEN  Nn- I 

15010  FOR  J-l  TO  Nn 

16020  IF  J-1  THEN  PRINT  *TL  0  2* 

16030  IF  J>l  THEN  PRINT  *TL  0  I* 

16040  Ya-Yat  »J 

16050  y-l6T< Ya/Ymin  >«Sfy 

16060  PRINT  -PD  PA  100,*,Y,*YT* 

16070  NEXT  J 
16080  NEXT  I 

16090  PRINT  'PA  100,100  PU* 

16100  PRINT  *PA  0,-2  5R  1.5,2* 

16110  Ii-LST(Xnin) 

16120  FOR  1-1  TO  C*-M 
16130  Xa-Xnin* 1 0' < I- I  ) 

16140  X-L6T<  Xa/Xmn  >»Sfx 
16150  PRINT  "PA* iX  .*  ,0i * 

16160  IF  ii>-0  THEN  PRINT  *CP  -2  ,-2 iLB 1 0 1 PR  -2  ,2 iLB* < 1 1 i *  * 

16170  IF  1 1 <0  THEN  PRINT  *CP  -2  ,-2 1 LB  1 0 iPR  0 ,2 iLB* i Ii  i  *  * 

16180  Ii-Ii*t 

16190  NEXT  I 

1 6200  PRINT  *PU  PA  0.0* 

16210  I i-LST ( Ymn ) 

16220  Y10-10 
16230  FOR  1-1  TO  Cy*1 
16240  Ya-Ymn- 1 0" ( I- 1  1 
16250  Y-L6T<Ya/Yrun)«Sfy 
16260  PRINT  *PA  0,*|Y,‘* 

16270  PRINT  *CP  -4  . 25 iLB 1 0 iPR  -2  .2 118* 1 1 u *  * 

16280  Ii-Ii*1 
16290  NEXT  I 
16300  BEEP 

16310  INPUT  -WANT  USE  OEFAUIT  LABELS  (  1 -Y ,0-N )7*  ,IdI 
16320  IF  IdlOl  THEN 
16330  BEEP 

16340  INPUT  "ENTER  x -LABEL*  .XIabel* 

16350  BEEP 

16360  INPUT  "ENTER  Y-LABEL*  ,YlaOel» 

16370  ENO  IF 

16380  IF  0oo<2  THEN 

16390  PRINT  *SR  I  ,2iPU  PA  40,-14t* 

16400  PRINT  *LB(T|PR  -1.6,3  PO  PR  1.2,0  PUtPR  .5 ,-4 iLBuoiPR  .5,1 i* 

15410  PRINT  *LB-T iPR  .5  ,-t iLBsat iPR  ,5,1»* 

16420  IF  Iun-0  THEN 
16430  PRINT  *LB )  (XI* 

16440  ELSE 

16450  PRINT  *LB  1  <  F  )  * 

16460  END  IF 
16470  END  IF 
16480  IF  Ooo-2  THEN 
16490  IF  Iun-0  THEN 

16500  PRINT  *SR  1  .5 ,2  iPU  PA  40,-14iLBq  (U/ntSR  l.I.SiPR  0.5,1  tLB2iSR  1.5,2iPR 
0.S.-1  iLB  1* 

16510  ELSE 

16520  PRINT  *SR  1 .5 ,2iPU  PA  34,-14|LBo  <Btu/hriPR  ,5..5tLB.iPR  .5. -.Si- 
16530  PRINT  *  LB  ft  i  PR  .S.IiSR  1,l.5iLB2iSR  1.5,2iPR  .S.-liLB)i* 

16540  ENO  IF 
16550  ENO  IF 
16560  IF  Ooo-0  THEN 
16570  IF  Iun-0  THEN 

16580  PRINT  *SR  1 . 5 ,2 iPU  PA  -!2,40iOI  0,1iLBq  (U/mPR  -l  ,0.5iSR  l.1.SiLB2iSR  I 
.5 ,2 iPR  l , . 5 1 LB  1  * 

16590  ELSE 

16600  PRINT  *SR  I  .5  ,2iPU  PA  - 1 2  ,32  <01  0,liL9o  iStu/hnPR  -  .5  .  .5iLB .  iPR  .5, .Si* 
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16610  PRINT  ”L0Tt i SR  l  , I . 5 1  PR  -  1  , ,5i^82 iPR  1  .  .Si  SR  1 .5 ,2:1.3  " 

16620  END  IF 
16630  END  IF 
18640  IF  Ooo>0  THEN 
16650  IF  lun-0  THEN 

16660'  PRINT  "SR  I .5 ,2iPU  PR-l2.38tOI  0.1iLBh  (W/mPR  -l..StSR  1,!.SilB2iSR  I 
.5  ,2  i  PR  .5. .Si* 

16670  PRINT  *SR  1  . 2  ,2 . 4  iPU  PR  -12  .3*7 1 01  B.liLBhiPR  1  ,0.5iL8oiPR  -1  ,0.5iL8  <U/n 

16680  PRINT  "PR  -1  ,  ,5iSR  1.1.  S  tl_B2 1  SR  1.5,2iPR  .S,.5tLB.«PR  .S,0tLBK>* 

16690  ELSE 

16700  PRINT  *SR  I .5  ,2iPU  PR  -1 2  .28 iDI  0,liL8h  <  Btu/Hr i PR  - .5  . . SiLB. iPR  .5. .Si* 

16710  PRINT  "LSf t iPR  -l.-SiSR  1.I.SiL82iSR  I . 5 . 2 i PR  .S..5iLB.iPR  .S..5iLBF>* 

16720  ENO  IF 

16730  END  IF 

16740  IF  Idl-0  THEN 

16750  PRINT  *SR  I.S.2iPU  PR  50,-16  CP* i -LEN( Xlabel* >/2 ; *0iLB* iXlabalSi *  * 

16760  PRINT  *PR  -14,50  CP  0 , * i -LEN< YlaOe 1« 1/2-S/6 1 *DI  0 , 1 iLB* i YladelCi *  * 

16770  PRINT  *CP  0,0  01* 

16780  END  IF 
16790  Ion-0 
16800  Xll-1  .E+6 
16810  Xul  — 1  .6+6 
16820  Icn-0 
16830  lfn-0 
16840  Ijoin-I 
16850  SEEP 

16860  INPUT  -URNT  TO  PLOT  DRTR  FROM  R  FILE  < 1 -Y  ,0-N >7*  ,0k 
16870  IF  OX- 1  THEN 
16880  SEEP 

16890  INPUT  -ENTER  THE  NAME  OF  THE  ORTR  FILE*  .D_file« 

16900  RSSIGN  0Fil«  TO  0_file» 

16910  SEEP 
16920  BEEP 

16930  INPUT  "ENTER  THE  BE3INNIN6  RUN  NUMBER* ,Md 
16940  BEEP 

16950  INPUT  'ENTER  THE  NUMBER  OF  X-Y  PRIRS  STORED* .Noairs 
169601  8EEP 

169701  INPUT  'CONNECT  ORTR  UITH  LINE  t 1 -Y , 0-N  >7 * , Ic I 

16900  BEEP 

16990  PRINTER  IS  1 

17000  PRINT  USIN6  * 4X  "Salad  ■  aynool:*** 

17010  PRINT  USIN6  *6X,*-|  Star  2  Plua  HQn*** 

17020  PRINT  USIN6  *6X  ,**3  Circle  4  Souere*** 

17030  PRINT  USIN6  *6X.**5  RonDus"* 

17040  PRINT  USIN6  *6X . " *6  Right-side-uo  triangle*"* 

17050  PRINT  USIN6  *SX  , * *7  Uo-»ide-do«n  triangle*** 

17060  INPUT  Syn 
17070  PRINTER  IS  705 
17080  PRINT  *PU  01* 

17090  IF  Syn-I  THEN  PRINT  *SH»* 

17)00  IF  Sy«-2  THEN  PRINT  *SH«* 

17110  IF  Syn-3  THEN  PRINT  *SMo* 

17120  IF  Md>1  THEN 
17130  FOR  1-1  TO  (Md-1  ) 

17140  ENTER  0Fi lei  Ye  .Xa 

17150  NEXT  I 

17160  END  IF 

17170  FOR  1-1  TO  Npeiri 

17180  ENTER  VF 1 1 e i y e  ,  X  a 

17190  IF  1-1  THEN  Ql-Ya 

17200  IF  1-Noaire  THEN  Q2-Ya 

17210  IF  Ooo-l  THEN  Ya-Ye/Xe 

17220  IF  Ooo-2  THEN 

17230  O-Ya 

17240  Ya*Ya/Xa 
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r:s®  <«*Q 

1726®  eno  if 

1727®  if  mom  then  xn-xa 

1728®  IF  Xa>Xul  THEN  Xul-Xa 
1729®  IF  Iun-1  THEN 

'  17300  IF  0oo<2  THEN  Xa-Xa-1.8 

1731®  IF  Ooo>0  THEN  Ya-Ya*.t76l 
17320  IF  Obo-0  THEN  Y«-Y«*.317 
17330  IF  Obo-2  THEN  Xa-Xa».3l? 

"  173*0  ENO  IF 

17350  X-LST1 Xa/X»in  )»Sf* 

1 7380  Y-lST(  Ya/Ymn  >»Sfy 
17370  Kj-0 

17380  CALL  Sy«b(X  ,Y  .Sy«  ,Icl  ,K j  ) 

17390  SOTO  17520 

17*00  IF  Sy«>3  THEN  PRINT  *S«* 

17*10  IF  Sy*<*  THEN  PRINT  'SR  1.*.2.*' 

17*20  IF  Ic 1*0  THEN 
17*30  PRINT  *PA*  .X  ,Y 
17**0  ELSE 

17*50  PRINT  'PA'  .X  ,Y  ,  "PO " 

17*60  ENO  IF 

17*70  IF  Syn>3  THEN  PRINT  "SR  1.2. '.6* 

17*80  IF  Sy««*  THEN  PRINT  *UC2  ,* . 99  ,0 . -8 ,-* ,® ,0 . 8 .*  .01  * 
17*90  IF  Sy«-S  THEN  PRINT  'UC3  ,0 .99  , -3  . -S  ,-3 , 6  ,3  ,6  ,3  ,-6 1 ' 
17500  IF  Sy«-6  THEN  PRINT  'UC0 ,5. 3 , 99 , 3 , -8  .-6 ,8  ,3  ,8 1 * 

1 75 1 0  IF  Sy««7  THEN  PRINT  'UC0  .-S.3  .99  ,-3  ,8  ,6  .0  .-3  ,-8 1 ' 
17520  NEXT  I 
17530  PRINT  -PIT 
'75*0  BEEP 

17550  INPUT  "WANT  TO  LA8EL7  ( 1 -Y  ,0-N  )*  .11*0 
17560  IF  I  lab-1  THEN 
17570  PRINT  'SP0|SP2‘ 

17580  BEEP 

17590  IF  K lab-0  THEN 
17600  X  lab-5 
17610  Y lab-85 

17620  INPUT  'ENTER  INITIAL  X  ,Y  LOCATIONS' ,XlaO .Ylab 
17630  XI t-Xlab-5 
176*0  Yt t-Ylab+8 
17650  PRINT  'SR  I .1 .S' 

17660  PRINT  'SHiPA* ,Xtt ,  Ytt ,'LB  X  Haat  Fila' 

17670  Yt*-Ytt-3 

17660  PRINT  'PA' ,Xtt  ,Ytt  ,'LB  Oil  Flux  Nana' 

17690  IF  Syw-1  THEN  PRINT  *SH»* 

17700  IF  Sy«-2  THEN  PRINT  'SH-' 

17710  IF  Syn-3  THEN  PRINT  'SHo' 

17720  Klab-I 
1773®  ENO  IF 
'77*0  K j  —  1 

17750  CALL  SynblXlab .Ylab  .Syn.Ici  ,Kj  ) 

17760  PRINT  'SR  1.1 . 5  <  SH' 

17770  IF  Syn<*  THEN  PRINT  'PR  2.0' 

17700  0EEP 

17790  INPUT  'ENTER  BOP' .Bob 

17800  IF  Boa< 10  THEN  PRINT  *PR  2  ,0tLB' i8ob  i  ” 

17810  IF  8oo>9  THEN  PRINT  'PR  . 5  ,0 iLB'  i Bob  i  '* 

«  '7820  Ihf-0 

17030  IF  01 >02  THEN  IMf- 1 

178*0  IF  Ihf-0  THEN  PRINT  'PR  *,®iLBInc‘ 

17850  IF  IHf-1  THEN  PRINT  'PR  *  ,0iLBObc' 

7  17860  PRINT  'PR  2 ,8iLB' tD_f ilati " 

1787®  PRINT  'SP®iSPl i SR  1.5.2' 

17080  Yiab-Ylab-5 
1709®  ENO  IF 
1790®  BEEP 
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17910  ASSI6N  IFilt  TO  • 

17929  xu-Xll/  I  .2 

17930  Xul-Xul»! .2 
179101  SOTO  8040 
17950  END  IF 

17960  PRINT  'PU  SM* 

17970  SEEP 

17980  INPUT  ‘WANT  TO  PLOT  A  POLYNOMIAL  ( 1 -Y , 0-N  >7 ■ ,50_en 

17990  IF  6o_on-1  THEN 

13000  BEEP 

18010  PRINTER  IS  1 

18020  PRINT  USIN6  '4X.'*Salact  lina  tyoa:'*' 

18030  PRINT  USIN6  '6X,"0  Solid  lma"' 

18040  PRINT  USIN6  '6X  ,*  *  1  Oaanad''' 

18050  PRINT  USIN6  ■6X1**2,..5  Longer  line  -  daeh"' 

18060  INPUT  Ipn 
18070  PRINTER  IS  705 
18080  BEEP 

18090  INPUT  'SELECT  < 0-LIN  .  I -L06< X  ,  Y  ) ) '  ,  I  log 

18100  Iprnt-1 

181 10  CALL  Poly 

18120  IF  Ifn-0  THEN 

18130  BEEP 

18140  INPUT  'ENTER  NUMBER  OF  FILES  TO  JOINT' ,Njom 
ieiS0  END  IF 
18160  Ijoin-0 

18170  IF  Ifn<Njoln  THEN  Ijoin-I 
18180  IF  Ifn-0  OR  Ijom-I  THEN 
18190  FOR  Ij-0  TO  3 
18200  BK  Ij  >-B»<  Ij  )+B<  I j  ) 

18210  NEXT  Ij 

18220  Ifn-Ifn+I 

18230  END  IF 

18240  IF  Njoin-Ifn  THEN 

18250  FOR  Ij-0  TO  3 

18260  B(  Ij  )*8»(  Ij  )/Njom 

18270  Bs< Ij  >-0 

18280  NEXT  Ij 

18290  Ifn-0 

18300  ELSE 

18310  60T0  16850 

18320  END  IF 

18330  BEEP 

18340  input  'ENTER  Y  LOWER  AND  UPPER  LIMITS' ,Y1 1  ,Yul 

18350  FOR  Xx-0  TO  Cx  STEP  Cx/200 

18360  Xa-Xmn-IO'Xx 

>8370  IF  Xa<Xll  OR  Xa>Xul  THEN  18640 

18380  Ien-Ien*l 

18390  Pu-0 

18400  IF  lon-1  THEN  Idf-Icn  MOD  2 
18410  if  Ion-2  THEN  Idf-Icn  MOD  4 
18420  IF  Ion-3  THEN  Idf-Icn  MOD  8 
18430  IF  Ipn-4  THEN  Idf-Icn  MOD  16 
18440  IF  !pn«5  THEN  Idf-Icn  MOO  28 
18450  IF  Idf-1  THEN  Pu-I 
18460  IF  Opo-0  THEN  Ya-FNPoly< Xa > 

18470  IF  Ooo-2  AND  I log-0  THEN  Ya-Xa/FNPoly( Xa > 

18480  IF  Opo-2  ANO  I  log- 1  THEN  Ya-FNPoly( Xa  ) 

18490  IF  Opo-1  THEN  Ya-FNPoly(Xa  l 

18500  IF  YaXYntn  THEN  18640 

18510  IF  Ya<Y 1 1  OR  Ya>Yul  THEN  18640 

18520  IF  Iun-I  THEN 

18530  IF  Opo<2  THEN  Xa-Xe»1.8 

18540  IF  Opo>0  THEN  Ya-Ya» .175! 

18550  IF  Opo-0  THEN  Ya-Ya» .317 

18560  IF  Opo-2  THEN  Xa-Xa«.317 
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:g='a  enc  :p 

1858®  v ««_6T t.  Ya/ thin  / * S r y 

1859®  X-i_6T<Xa/Xmn>-Sfx 

1 880®  IF  Y<3  THEN  Y-0 

1861®  IF  Y> 1 00  THEN  SOTO  18640 

18620  IF  Pu-0  THEN  PRINT  *PA- ,X  ,Y  .'PO* 

18630  IF  Pu-I  THEN  PRINT  'PA*  ,X  ,Y  ,"PU" 

1864®  NEXT  Xx 
18650  PRINT  *P U- 
18660  SOTO  1 S850 
18670  END  IF 
18680  SEEP 

18690  INPUT  “WANT  TO  PLOT  REILLY'S  DATA?  I l-Y  ,0-N >■  , Irly 

1870®  IF  Ooo-0  OR  Opo-1  THEN 

1871®  XU-3 

18720  Xul-20 

18730  ENO  IF 

1874®  IF  Opo-2  THEN 

18750  XU- 10000 

18760  Xul-100000 

18770  ENO  IF 

18780  IF  Irly-I  THEN 

18790  YU-20 

18800  Yul-70 

18810  BEEP 

18820  INPUT  -ENTER  LOHER  ANO  UPPER  Y-lIHITS  FOR  PL0TTIN6'  ,Yil  ,Yul 

18830  FOR  Xx-0  TO  Cx  STEP  Cx/200 

18840  Xa-Xn m*  1 0"Xx 

18850  IF  Xa<XU  OR  XaXXui  THEN  18980 

1 8860  X 1 -l06( Xa  ) 

18870  IF  Opo-0  THEN  Yl— 2.5402837E-I*X1»(4.9720151-XI»2.5134787E-I  ) 
18880  IF  Opo-I  THEN  Y 1— 2  .S402837E-1  *X1  •(  3 . 97201 51  -X 1  -2 .51 34787E- 1  1 
18890  IF  Opo-2  THEN  Yl— 3 . 7073801E- 1  ♦*  1  •<  8 . 72591 90E- 1  -X I  -6 . 08Z6842E-3  > 
I89ea  Ya-EXP<  Yl  ) 

18910  Y-L6T(Ya/Ynin.>-Sfy 
18920  X-l6T<  Xa/Xnin  )«Sf x 
18930  Ipu-0 

18940  IF  Y<YU  THEN  Iou-1 
18950  IF  Y>Yul  THEN  60TO  18980 
18960  IF  Ipu-0  THEN  PRINT  *PA’  ,X  ,Y . 'PO- 
18970  IF  Ipu-I  THEN  PRINT  ’PA’ .X ,Y . -PU* 

It  1  NEXT  Xx 
16. -0  PRINT  -PU- 
19000  ENO  IF 
19010  BEEP 

19020  INPUT  -WANT  TO  PLOT  ROHSENOw  CORRELATION?  ( 1 -Y  ,0-N ) *  , Iroha 

19030  IF  Iroha- 1  THEN 

19040  YU-15 

19050  Yul-00 

<9060  SEEP 

1 9070  INPUT  -ENTER  Taat  <  Oag  Cl-.Taat 

19080  Ca'-.0040 
19090  SEEP 

19100  INPUT  -ENTER  Caf  (DEF-0.004 >*  ,Caf 
19110  Tf -Taat+2 

19120  FOR  Xx-0  TO  Cx  STEP  Cx/200 
19130  Xa-Xmh«10"Xx 
19140  IF  xa<xu  OR  XaXXuI  THEN  >9420 
19150  xi-LOS(Xa) 

1916®  IF  Opo<2  THEN  Tf»Taat*Xa/2 
19170  Rho-FNRho( T f  ) 

19180  K-FHX<  Tf  ) 

19190  Hu-FNHulTf  ) 

19200  Cp-FNCo(Tf) 

19210  Hf Q-FNHf q<  T  aat  ) 

19220  Ni-Hu/Rho 
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' 3;33  --•Cd»Hu - k 

>9240  Onega-Cs I •*fg/ Co* v '  . 0’ ■ 9 • 3 ' *Rnc  .  .3. 

19250  IF  Qoo-0  THEN  Ya«(Xa/0«aga >"3 

192S0  IF  0oo-l  THEN  Ya-( Xa/0«ega  )  "3/Xa 

19270  IF  Ooo-2  THEN  Ya-Xa' ( 2 . /3 1/Ontga 

19290  IF  Obo-2  THEN 

19290  Tfe*Tsat*Xa/Ya*.5 

19300  IF  ABS<Tf-Tfc  >>.01  THEN 

19310  Tf»(Tf*Tfc  )•  .5 

19320  SOTO  19170 

19330  ENO  IF 

19340  ENO  IF 

19350  Y-L6T(  Ya/Ynm)»Sfy 
19360  X>L6T(Xa/Xnin)*Sfx 
19370  Ipu-0 

19390  IF  Y<Y11  THEN  Ipu-1 
19390  IF  Y>Yul  THEN  19420 
19400  IF  Iou-0  THEN  PRINT  "PA- .X  .Y.'PD* 

19410  IF  Ipu-t  THEN  PRINT  ‘PA- ,X . Y  .  ■ PU* 

19420  NEXT  Xx 
19430  PRINT  ‘f>U‘ 

19440  ENO  IF 
19450  BEEP 

19460  INPUT  "WANT  TO  QUIT  ( I -Y  ,0-N )• , Iqt 

19470  IF  lot- 1  THEN  19490 

19480  60T0  16800 

19490  PRINT  ’PU  PA  0,0  SP0" 

19500  SUBENO 

19510  SUB  Synb<X  .Y.Sym.Icl  ,Kj  ) 

19520  IF  5y«>3  THEN  PRINT  -SH' 

19530  IF  Sy«<4  THEN  PRINT  *SR  1.4.2.4" 

I  9540  Y«d-0 

19550  IF  Kj-1  THEN  Yad-.8 
19560  IF  Icl-0  THEN 
I  9570  PRINT  "PA*  ,X  ,Y*Yad ,  *  * 

19580  ELSE 

19590  PRINT  -PA- ,X .Y+Yad.'PO" 

I  9600  ENO  IF 

19610  IF  Syn>3  THEN  PRINT  "SR  1.2. 1.6* 

19620  IF  Syn-4  THEN  PRINT  'UC2 ,4 ,99 ,0 ,-8 .-4 ,0 ,0 .8 .4 ,0i * 
19630  IF  Syn-5  THEN  PRINT  *UC3 .0 , 99 .-3 ,-6 ,-3 ,6  ,3  .6  .3  ,-S ; * 
19640  IF  Syn-6  THEN  PRINT  'UC0 .5.3 .99 ,3 ,-B ,-6 .0  .3  ,8 1 * 

19650  IF  Sy«-7  THEN  PRINT  'UC0 . -5. 3 , 99 , -3 .8 ,6 .0 .-3  . -8 1 * 
19660  IF  Kj-1  THEN  PRINT  "SfliPR  0.-.8- 
19670  SUBENO 
1 9680  SUB  Purg 
19690  BEEP 

19700  input  -ENTER  file  NAME  TO  BE  DELETED' ,F i la* 

19710  PUR6E  Filat 
19720  SOTO  19690 
19730  SUBENO 
19740  SUB  Tdcn 
19750  COM  /Ce/  C(7).leal 
19760  DIM  Enf ( 1  > 

19770  DATA  0.10086091.25727.94369.-767345.8295.78025595.81 
19780  OATA  -9247486589.6.976886*11  , -2. 66I92E-I3. 3. 94078E-M4 
19790  READ  CIO 
19800  BEEP 

19810  INPUT  -6IVE  A  NAME  FOR  FILE  TO  8E  CREATED- ,Fi iat 
19820  BEEP 

19830  INPUT  'SELECT  TUBE  < 0-HH , t -HF ,2-WT >• , I t t 
19840  SEEP 

19850  INPUT  -SELECT  THERMOCOUPLE  TYPE  ( 0-NEW , 1 -OLD )' .leal 
19860  IF  1 1 1 <2  THEN  Oi-.0127 
19870  CREATE  BOAT  Filat, 4 
19880  ASS16N  BFila  TO  Fila« 
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3890  jU’PU"  :  ie  i  1 1  •• 

19900  J-0 
19910  BEEP 

19920  INPUT  "ENTER  MONTH,  DATE  AND  TIME  <MM:DO:HH:MM:SS )* .Oate* 

19930  OUTPUT  709i*TD’iOate0 
19940  OUTPUT  709i ’TD" 

19950  ENTER  709. Date* 

19960  PRINTER  IS  I 
19970  PRINT 

19980  PRINT  ‘  Month,  date  and  tine:  'lOateC 

19990  PRINT 

20000  PRINT  USINS  ’10X,"  Fns  Tin  Tav  U«* 2  Tdrop’” 

20010  IP  K-0  THEN 
20020  PRINTER  IS  701 
20030  PRINT 

20040  PRINT  ’  Month,  date  and  tine:  ’tOataC 

20050  IF  Itt-0  THEN  PRINT  USINS  ’t0X.’’Tube  Type:  Uieland  Snooth’’’ 

20060  IF  Itt-1  THEN  PRINT  USIN6  ’10X,’*Tube  Type:  Hi0h  Flux’" 

20070  IF  Itt-2  THEN  PRINT  USIN6  '  1 0X  . "Tupe  Type:  Turbo-B”’ 

20080  PRINT 

20090  PRINT  USINS  ’10X,"  Fns  Tin  Tev  Vu*2  Tdrop’" 

20100  PRINTER  IS  1 
20110  K-1 
20120  END  IF 
20130  BEEP 

20140  INPUT  'ENTER  FLOWMETER  READING’  ,Fna 
20150  OUTPUT  709 1 ’ AR  AF0  AL4  UR  1 ’ 

20160  FOR  L-0  TO  4 
20170  OUTPUT  709 i ’AS  SA’ 

20180  IF  L>0  AND  L<4  THEN  20260 
20190  S-0 

20200  FOR  1-0  TO  9 
20210  ENTER  709 i E 
20220  S-S+E 
20230  NEXT  I 

20240  IF  L-0  THEN  Enf i 0 )-ABS( S/ 1 0 > 

20250  IF  L-4  THEN  Enf ( 1  >-ABS< S/ 1 0 > 

20260  NEXT  L 

20270  OUTPUT  709 1  * AR  AF00  AL00  UR1* 

20280  OUTPUT  709 i ’AS  SA’ 

20290  Etp-0 

20300  FOR  1-0  TO  9 

20310  ENTER  709iEt 

20320  Etp-Etp-Et 

20330  NEXT  I 

20340  Eto-Etp/10 

20350  Tin-FNTvev(Enf( 1  ) ) 

20360  Tew-FNTvev( Enf < 0  )  ) 

20370  6rad*37 . 9853+ .  1 04308»Tm 

20380  Mdot«3 . 9657E-3nFn»«<  3.S 1 955E*3-Fns»<  8. B200SE-6-Fns»(  I .23688E-7-Fns»4.3!897 
E- 10  > ) ) 

20390  Uu-Mdot/< I000«PI*Oi  ’2  >«4 
20400  Tdrop-Etp»1 .E»S/( 10»6rad ! 

20410  PRINT  USINS  ’ I0X ,3< 00.00  ,4X  ) , I X ,2 .00  ,4X ,MZ . 40’ iFna .Tin  ,Tev  ,V«*2 .Tdrop 

20420  8EEP 

20430  INPUT  ’WANT  TO  ACCEPT  THIS  DATA  SET’  ( I -V ,0-N )* ,0X 

20440  J.J  +  1 

20450  IF  Ok-0  THEN 

20460  J-J-1 

20470  SOTO  20130 

20480  ELSE 

20490  OUTPUT  0F lie i Fns  ,Enf<  •  )  ,Etp 
20500  PRINTER  IS  701 

20510  PRINT  USINS  ’10X .3(00. 00, 4X ), IX  ,2.00, 4X  , M2. 40’iFna  .Tin, Tev, Uu*2, Tdrop 
20520  PRINTER  IS  1 
20530  BEEP 
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20543  INP'j*  'wlui-  7m£k£  BE  AnC~hE5  OATH  £E“"  ■  *  ,0«n  5- 

:055a  IF  6o_on-l  THEN  Z0120 

70560  END  IF 

20570  ASSISN  IFiit  TO  • 

205S0  PRINTER  IS  701 

20590  PRINT  » 

20600  PRINT  US  INS  *  1 0X  ,  ‘"NOTE:  **.ZZ.**  data  tats  art  stored  m  file  **,1SA*iJ,F 

1  let 

20610  PRINTER  IS  1 
20620  SU8END 
20630  SUB  Uoprt 
20640  PRINTER  IS  1 
20650  BEEP 

20660  INPUT  'Entar  Uo  File  Nana*  ,File» 

20670  BEEP 

20680  INPUT  "NunOer  of  Data  Runs* , Nrun 

20690  INPUT  *0o  You  Want  a  Plot  File’d-Y ,0-N>* .Iolot 

20700  BEEP 

20710  IF  Iplot-1  THEN 

20720  INPUT  *6tve  Plot  File  Nana*  ,P_f lla« 

20730  CREATE  BOAT  P_file«,4 
20740  ASSI6N  BPlot  TO  P_fila« 

20750  END  IF 
20760  PRINTER  IS  701 
20770  PRINT 
20780  PRINT 

20790  PRINT  US  INS  *  1 0X  ."*  Water  Ual  Uo*** 

20800  ASSI6N  0File  TO  File* 

20810  IF  Iplot-1  THEN  ASSISN  IFilet  TO  P.file* 

20820  FOR  1-1  TO  Nrun 
20830  ENTER  •FilaiUu.Uo 

20840  IF  Iplot-1  THEN  OUTPUT  BFiletiUw.Uo 
20850  PRINT  USIN6  *  1 SX .0 . DO  ,6X  ,MZ .306* iVu .Uo 
20860  NEXT  I 
20870  ASSISN  tFila  TO  • 

20880  ASSI6N  #F 1 1 a  I  TO  • 

20890  PRINT  USIN6  *I0X,**NOTE:  **.ZZ.**  data  aata  era  atorad  in  file  **.l5A*iNru 
n ,Fi let 

20900  IF  Iplot-1  THEN 

20910  PRINT  USIN6  " 10X "NOTE:  **,ZZ."  X-Y  Pairs  are  stored  in  file  **.15A*!Nru 
n,P_file* 

20920  ENO  IF 
20930  PRINTER  IS  I 
20940  SUBENO 
20950  SUB  Select 
20960  COM  /I  dp/  I  dp 
20970  BEEP 
20980  PRINTER  IS  1 

20990  PRINT  US  I  NS  *4X.  "Select  option:*** 

21000  PRINT  US  I  NS  *6X,**  0  Tallinn  data  or  ra-groctiiinj  previous  date*** 

21010  PRINT  US  INS  *6X  ,  * ’  I  Plotting  data  on  Log-Log  *** 

21020  PRINT  USINS  *6X  ,  *  *  2  Plotting  data  on  Linear*** 

21030  PRINT  USINS  *6X  ,  *  *  3  HaKe  cross-plot  eoefft  fils*** 

21040  PRINT  USIN6  *6X  ,**  4  Re-eireulete  uatar*** 

21050  PRINT  USINS  *6X  , *  *  5  Purge*** 

21060  PRINT  USINS  *6X  ,**  6  T-Orop  correct  ion* *  * 

21070  PRINT  USINS  *6X  ,  *  *  7  Print  Uo  File*’* 

2 1080  PRINT  USIN6  *6X,**  B  Modify  x-Y  file***  ( 

21090  PRINT  USINS  *6X  ,  *  *  9  Hove*** 

21100  PRINT  USINS  *6X .*  *  10  Cono/F  ixup *  *  * 

21110  INPUT  Ido 

2H20  IF  Idp-0  THEN  CALL  Mein 

21130  IF  I  dp- 1  THEN  CALL  Plo* 

21140  IF  I  dp-2  THEN  CALL  Pirn 
21150  IF  I dp-3  THEN  CALL  Coef 
21160  IF  I dp-4  THEN  CALL  Main 
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Z'‘~?  IF  lcomS  'n£N  Cnu-  -urg 
2 '  ;  50  IF  I ao “6  'nEN  CAuw  Tacn 

21190  IF  l dp-7  THEN  CALL  Uoprt 
21200  IF  I  dp-8  THEN  CALL  Xynod 
21210  IF  ldp-9  THEN  CALL  Nova 
*  21220  IF  I dp- 10  THEN  CALL  Conb 

21230  SUBEND 
21240  SUB  Xynod 
2I2S0  PRINTER  IS  1 
-  21260  BEEP 

21270  INPUT  ‘ENTER  FILE  NAME‘,Fila» 

21280  ASSISN  0F 1 1  a  I  TO  FUal 
21290  BEEP 

21300  INPUT  ‘ENTER  NUMBER  OF  X-Y  PAIRS'  ,Np 
21310  BEEP 

21320  INPUT  ‘ENTER  NEW  FILE  NAME‘,File2* 

21330  CREATE  BOAT  Fila20.S 
21340  ASSI6N  0FiIa2  TO  Fila2* 

21350  BEEP 

2 1360  INPUT  ‘ENTER  NUMBER  OF  X-Y  PAIRS  TO  BE  DELETED"  .Ndal 
21370  IF  Ndal-0  THEN  21410 
21380  FOR  1-1  TO  Ndal 
21390  BEEP 

21400  INPUT  ‘ENTER  DATA  SET  NUMBER  TO  BE  DELETED*  ,Nd( I > 

21410  NEXT  I 

21420  FOR  J-1  TO  Np 

21430  ENTER  0Filal iX  ,Y 

21440  FOR  1-1  TO  Ndal 

21450  IF  Nd< 1 1-1  THEN  21490 

21460  NEXT  1 

21470  OUTPUT  0Fila2lX,Y 

21400  PRINT  J.X.Y 

21490  NEXT  J 

21500  PRINTER  IS  701  . 

21510  ASSI6N  0F 1 1 a  I  TO  • 

21520  ASSISN  0Fila2  TO  . 

21530  SUBENO 
21540  SUB  Mova 
215501  FILE  NAME:  MOUE 
215601 

2  1 570  OIM  Bap <  66  )  ,A(66  >  ,8<66>.C(66  >  .0(66)  ,E(66  )  ,F<66  >  ,6(66  >  ,H<  66  ) ,  J<  66  )  ,K<  65  )  ,L( 
66 1 ,M< 66  ) 

21580  OIM  Told*<66  It  1 4 1  ,N<66  >  ,Vr<  66  >  ,Ir<  66  > 

21590  BEEP 

21600  INPUT  ‘OLO  FILE  TO  MOUE‘ ,02_fUa» 

21610  ASSISN  0Fila2  TO  02_Mla* 

21620  ENTER  0Fi la2 iNrun  .Data*  .Ldtci .Ldte2 ,Itt 

21630  FOR  I-l  TO  Nrun 

21640  ENTER  0Fila2iBop< I  )  ,Told*< I  ) 

21650  ENTER  0Fila2iA(  I )  ,B(  I ) ,C(  1  >  ,D(  I  )  ,E<  I )  ,F(  I )  ,6(  I )  .H(  I )  ,J(  I )  ,K(  I )  ,L<  I  >  ,M<  I )  ,N 

( I  1 

21660  ENTER  0Fila2(Vr< I >  ,Ir( I  ) 

21670  NEXT  I 

21680  ASSISN  *Fila2  TO  • 

21690  BEEP 

21700  INPUT  ‘SHIFT  DISK  ANO  HIT  CONTINUE’ . OX 
21710  BEEP 

«  21720  INPUT  ‘INPUT  BOAT  SIZE’  .Siza 

21730  CREATE  BOAT  02_ftla*  .Sue 
21740  ASSI6N  *Filal  TO  02_f i la* 

21750  OUTPUT  0Fi lal iNrun  .Datal.Ldtcl  ,Ldtc2  .1 tt 
*>  21760  FOR  1-1  TO  Nrun 

21770  OUTPUT  BFilal (Bop<  Z  > ,  To  I d* ( I  > 

21780  OUTPUT  0F 1 1  a  1  iA<  I  1  ,d(  1 )  ,C(  I  1  ,D(  I  )  ,E(  I  )  ,F(  I  )  ,6(  I )  ,H(  I )  ,J(  I  )  ,K(  I )  ,L(  I  >  ,M(  1 1  , 
N<  I  ) 

21790  OUTPUT  0Filal iUr( I  )  ,Ir< I  ) 
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2! see  NEXT  ! 

21818  ASSIGN  CFilel  TO  • 

2 1820 1  RENAME  -TEST*  TO  D2_fila* 

2)838  BEEP  2088.-2 
21848  SEEP  4088,-2 
21858  SEEP  4880.-2 
2 1860  PRINT  "DATA  FILE  MOUEO' 

21878  SU8ENQ 
21680  SU8  Conb 
21898'  FILE  NAME:  COMB 
21988* 

21910  DIM  Enf< 12) 

21928  SEEP 

21930  INPUT  'OLO  FILE  TO  FIXUP' ,D2_fxia* 

21940  ASSIGN  0Fxla2  TO  02_fxla» 

21950  Oi_fila*-'TEST' 

21968  CREATE  BOAT  01  fila*,20 
21970  ASSIGN  fFilal  TO  01_fila» 

21980  ENTER  #Fila2iNrun .Date* .Ldtel ,Ldtc2 , 1 1 1 
21990  Nrunn«20 

22000  IF  K-0  THEN  OUTPUT  0F i La  1 iNrunn  .Data*  ,Ldte 1  .Ldtc2,Itt 
22010  FOR  1-1  TO  Nrun 

22020  ENTER  0Fi la2 iBop .Told*  .Enf ( •  )  ,Ur .Ir 
22030  OUTPUT  8FilaliBoo. Tola*, Enf<»>. Ur. Ir 
22040  NEXT  I 

22050  ASSIGN  0Fila2  TO  • 

22060'  RENAME  'TEST'  TO  02_fxla* 

22070  BEEP  2000.-2 
22080  BEEP  4008.-2 
22090  BEEP  400B..2 
22100  BEEP 

22110  INPUT  'UANT  TO  ADO  ANOTHER  FILE  ( I -V  ,0»N  )7 *  .Oka 
22120  IF  Oka-1  THEN 
22130  K-1 
22140  BEEP 

22150  INPUT  *6IUE  NEW  FILE  NAME ' ,Nf i let 
22160  ASSI6N  0Fila2  TO  Nfilet 
22170  GOTO  21980 
22180  ENO  IF 

22190  ASSI6N  0Fiie2  TO  • 

22200  SU8END 


Oate  :  19  Jan  1993 

NOTE:  ProQran  nana  :  0RP6B 
Oiak  number  •  00 
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